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PREFACE 


The  Engineering  Design  Handbook  Series  of  the  Army  Materiel  Command 
is  a  coordinated  series  of  handbooks  containing  basic  information  and  fundamental 
data  useful  in  the  design  and  development  of  Army  materiel  and  systems.  The 
handbooks  are  authoritative  reference  books  cf  practical  information  and  quanti¬ 
tative  facts  helpful  in  the  design  and  development  of  Army  materiel  so  that  it 
will  meet  the  tactical  and  the  technical  needs  of  the  Armed  Forces. 

This  handbook.  Aerodynamics,  is  one  of  a  series  on  ballistic  missiles.  It  presents 
design  information  peculiar  to  ballistic  missile  flight,  i.e.,  information  outside 
the  scope  of  basic  subsonic  and  supersonic  aerodynamics  which  may  be  termed 
“conventional  aerodynamics.”  Accordingly,  this  handbook  treats  the  problems 
unique  to  ballistic  flight,  namely,  high-speed,  high  temperature  and  high  altitude 
effects.  Empirical  as  well  as  theoretical  data  are  presented  in  the  tables  and  graphs 
which  accompany  the  text  material.  As  of  the  cut-off  date  for  publication,  the 
data  and  theory  are  believed  to  be  current,  representing  the  latest  thinking  and 
observations  in  the  areas  cf  atmospheric  and  geophysical  parameters,  typical 
ballistic  trajectories,  physical  properties  of  air  particularly  at  elevated  temper¬ 
atures,  pressure  distributions  over  bodies  at  high  speeds  as  required  for  determining 
aerodynamic  loads  and  heating  effects,  aerodynamic  heating  and  cooling  methods, 
hypersonic  aerodynamic  forces  and  stability  considerations,  and  techniques  and 
problems  of  aerodynamic  experimentation.  Appropriate  engineering  approximations 
have  been  introduced,  where  applicable,  to  simplify  the  complex  theory  and  make 
the  text  adaptable  to  the  solution  of  practical  problems. 

The  original  manuscript  for  this  handbook  was  supplied  by  Vitro  Laboratories 
for  the  Engineering  Handbook  Office  of  Duke  University,  prime  contractor  to 
the  Army  Research  Office-Durham. 

Elements  of  the  U.S.  Army  Materiel  Command  having  need  for  handbooks 
may  submit  requisitions  or  official  requests  directly  to  Publications  and  Reproduc¬ 
tion  Agency,  Letterltenny  Army  Depot,  Chambersburg,  Pennsylvania  17201. 
Contractors  should  submit  such  requisitions  or  requests  to  their  contracting 
officers. 

Comments  and  suggestions  on  this  handbook  are  welcome  and  should  be 
addressed  to  Army  Research  Office-Durham,  Box  CM,  Duke  Station.  Durham, 
North  Carolina  27706. 
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CHAPTER  1* 

AERODYNAMICS  IN  BALLISTIC  MISSILE  DESIGN 


1-1  AERODYNAMIC  PROBLEM  AREAS 

The  design  and  specification  of  a  ballistic  mis¬ 
sile  of  one  or  more  stages  require  areodynamic  in¬ 
formation  somewhat  wider  in  scope  than  that  usually 
required  for  totally-guided  missiles.  This  is  due 
mainly  to  enlarged  operational  conditions.  Of  great 
importance  for  ballistic  missiles,  for  instance,  are 
the  thermal  effects  engendered  by  high  velocities, 
particularly  during  re-entry  into  the  higher  density' 
regime  cf  the  atmosphere.  Furthermore,  ballistic 
missiles  encounter  significant  variations  in  atmos¬ 
pheric  conditions  as  they  traverse  wide  altitude 
ranges.  In  addition,  aerodynamic  factors  which  in¬ 
fluence  the  precision  of  aiming  and  launching  be¬ 
come  important  since  trajectories  cannot  readily  be 
altered  after  ballistic  flight  has  begun. 

Aerodynamic  phenomena  bear  on  numerous 
facets  cf  missile  design  and  operation : 

(1)  Overall  performance  and  flight  mechanics  are 
influenced  by  resultant  aerodynamic  forces. 

(2)  Stability  and  control  depend  upon  aero¬ 
dynamic  moments  and  stability  derivatives. 

(3)  Structural  design  in  the  low  temperature  range 
depends  upon  steady  and  time-dependent 
aerodynamic  loads. 

(4)  Material  selection  and  structural  design  in  the 
high  temperature  range  depend  upon  a  know¬ 
ledge  cf  aerody  namic  heat  inputs  and  losses, 
as  well  as  cooling  techniques.  Factors  of  this 
type  dominate  the  re-entry  capabilities  cf  high 
performance  missiles.  In  fact,  trajectories  may 
be  programmed  to  alleviate  high  temperature 
effects. 

(5)  The  effectiveness  cf  detection,  discrimination, 
tracking  and  interception  of  high  speed  mis¬ 
siles,  as  well  as  communication  to  or  from 
the  missile,  bear  strongly  on  the  physical  state 
of  high  temperature  gases  in  the  immediate 
vicinity  of  the  vehicle  and  in  its  wake.  In 
particular,  electron  concentrations  and  radia- 

*  This  volume  was  written  by  M.  H.  Bloom.  S.  V. 
Nardo  aided  with  the  chapter  on  aerodynamic  testing,  and 
P.  M.  Sforza  assisted  in  arrangement  and  editing.  All  are 
from  Poh  technic  Institute  cf  Brooklyn.  C.  D.  Fitz.  then  of 
Vitro  Laboratories,  presented  constructive  suggestions. 


tion  properties  are  significant  in  this  regard. 

(6)  Aerodynamic  experimentation  and  testing, 
and  the  establishment  cf  model-scaling  proce¬ 
dures  in  connection  with  a  given  design  are 
important  problem  areas  in  themselves. 

(7)  The  magnitude  cf  the  dispersion  cf  the  missile 
at  impact  may  be  influenced  by  various  aero¬ 
dynamic  disturbances.  These  disturbances 
may  be  introduced  by  meteorological  condi¬ 
tions,  as  well  as  by  the  aerodynamic  con¬ 
sequences  cf  structural  inaccuracies. 

(8)  Aeroelastic  phenomena  such  as  flutter  depend 
upon  steady  and  time-dependent  aerodynamic 
loads,  and  upon  aerodynamic  derivatives 
which  act  as  forcing  and  damping  functions. 

The  process  of  designing,  building  and  using 
a  missile  is  a  “closed-loop”  system  in  which  in¬ 
formation  is  constantly  fed  back  and  forth  between 
various  phases  of  the  system,  engendering  correc¬ 
tions.  modifications,  and  compromises  in  each  phase. 
Thus  requirements  set  by  other  than  aerodynamic 
factors,  such  as  structural  design  or  handling,  may 
alter  the  ideal  aerodynamic  configuration.  An  im¬ 
portant  role  of  aerodynamics  is  to  aid  in  the  evalua¬ 
tion  of  the  effects  of  these  changes. 

Furthermore,  meteorological  information  is  re¬ 
quired  for  the  determination  of  aerodynamic  phe¬ 
nomena.  Aside  from  the  usual  atmospheric  charac¬ 
teristics,  features  such  as  winds,,  gusts  and  high 
altitude  properties  may  be  cf  importance.  Moreover, 
phenomena  such  as  bombardment  by  meteoric  or 
electrically  charged  particles  may  influence  the  con¬ 
dition  cf  aerodynamic  surfaces  and  consequently 
the  aerodynamic  ballistics. 

1-2  OBJECTIVES  CF  AERODYNAMICS 
HANDBOOK 

Significant  amounts  cf  information  in  the  areas 
previously  cited  lie  within  the  scope  of  basic  sub¬ 
sonic  and  supersonic  aerodynamics,  and  are  not 
peculiar  to  ballistic  missile  flight.  With  some  liberty, 
these  aspects  may  be  termed  “conventional  aero¬ 
dynamics.”  On  the  other  hand,  high-speed,  high- 
temperature  and  high-altitude  effects  are  essentially 
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unique  to  ballistic  flight.  As  a  result  they  receive 
major  emphasis  in  this  chapter. 

Atmospheric  and  geophysical  parameters  are 
discussed  in  Chapter  2;  typical  ballistic  trajectories 
in  Chapter  3;  physical  properties  of  air.  particularly 
at  elevated  temperatures,  in  Chapter  4;  pressure 
distributions  over  bodies  at  high  speeds,  as  required 
for  determining  aerodynamic  loads  and  aerodynamic 


heating  effects,  in  Chapter  5 ;  aerodynamic  heating 
and  some  cooling  methods  in  Chapter  6;  hyper¬ 
sonic  aerodynamic  forces  and  stability  considera¬ 
tions  in  Chapter  7;  and  techniques  and  problems 
of  aerodynamic  experimentation  in  Chapter  8. 

Basic  references  on  aerodynamics  and  tabula¬ 
tions  in  general  use  are  listed  at  the  end  of  this 
chapter. 
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CHAPTER  2 

ATMOSPHERIC  AND  GEOPHYSICAL  PROPERTIES 


2-1  INTRODUCTION 

The  aerodynamic  characteristics  and  dynamic 
behavior  of  a  vehicle  in  the  vicinity  of  a  planet  such 
as  the  earth,  depend  upon  the  nature  of  the  atmos¬ 
phere  and  upon  geophysical  parameters  such  as  the 
planet’s  gravitational  field. 

Standards  for  the  principal  atmospheric  and 
geophysical  properties  cf  the  earth  have  been  set  up 
by  the  Air  Research  and  Development  Command 
(ARDC),  the  U.  S.  Air  Force1'2  and  more  recently 
by  the  U.S.  Committee  on  Extension  to  the  Standard 
Atmosphere  (COESA)20  after  extensive  study  of 
existing  information.  Reviews  are  made  annually  as 
new  information  is  acquired.  These  Standards  pro¬ 
vide  a  unified  basis  for  reference  and  are  recom¬ 
mended  for  initial  design. 

In  actuality,  deviations  from  standard  condi¬ 
tions  occur  due  to  such  factors  as  geographic  loca¬ 
tion,  seasonal  changes,  time  of  day  and  weather. 
Their  effect  may  be  significant  for  purposes  cf 
final  design  and  detailed  performance.  However, 
their  importance  depends  upon  the  magnitude  of 
their  effect  on  vehicle  parameters,  the  required 
reliability  cf  the  vehicle  and  the  probability  that 
deviations  from  standard  will  occur. 

At  altitudes  below  46  km  (150,000  ft)  atmos¬ 
pheric  data  are  obtained  quite  accurately  with  the 
use  of  balloons.  At  higher  altitudes,  data  are  obtained 
by  means  of  rockets,  satellites  and  remote  observa¬ 
tions.  Theoretical  analyses  are  used  to  correlate  the 
data  and  provide  extrapolations  to  extremely  high 
altitudes  (700  km). 

The  1959  ARDC  Standard  takes  into  account 
density  data  obtained  by  rocket-borne  ionization 
gages  in  the  110  to  220  km  range.  It  also  considers 
the  more  uncertain  densities  inferred  from  the  change 
in  period  of  satellites  having  perigee  (minimum) 
altitudes  of  170  to  650  km.  Below  20  km,  the  1959 
Standard  agrees  with  prior  Standards  such  as  that 
of  NACA3.  Below  53  km  this  Standard  agrees  with 
the  1956  ARDC  model.  Between  90  and  150  km 
the  values  of  density  are  about  one-half  the  1956 
values.  At  200  km,  the  density  is  about  5  times 
the  prior  Standard,  and  at  600  km  the  1959  value 
is  about  20  times  the  1956  Standard. 


The  U.S.  Standard  Atmosphere  is  based  on  more 
recent  data  acquired  in  such  additional  diverse 
experiments  as  rocket  grenade  launchings,  falling 
sphere  density  measurements,  etc.  Below  20  km  the 
1962  U.S.  Standard  Atmosphere  agrees  with  the  1959 
ARDC  Atmosphere  and  the  International  Civil 
Aviation  Organization  (ICAO)  Atmosphere.  Up  to 
80  km  the  1959  Standard  is  still  in  agreement  with 
the  1962  Standard.  Above  80  km  the  differences 
between  the  two  are  depicted  graphically  in  Figs. 
2-1  (Q)  and  2-1  (R). 

Probable  errors  in  pressure  and  density  are 
estimated  to  be  less  than  10  percent  at  altitudes 
below  60  km  (200,000ft),  and  to  increase  to  perhaps 
200-300  percent  at  200  km4. 

2-2  PROPERTIES  OF  THE  EARTH’S 
ATMOSPHERE 

For  the  purposes  of  this  handbook,  the  prop¬ 
erties  of  the  earth’s  atmosphere  as  developed  for  the 
1962  U.S.  Standard  Atmosphere  and,  secondarily, 
the  1959  ARDC  Standard  Atmosphere  will  be 
adopted.  Table  2-1  presents  a  list  of  the  notations 
used  in  both  Standards2  ,20.  Sea-level  values  are  given 
in  Table  2-2  in  both  metric  and  English  units. 
Values  of  the  physical  constants  used  in  deriving 
the  Standard  are  shown  in  Table  2-3.  Pertinent 
conversion  factors  are  given  in  Table  2-4.  Table  2-5 
lists  the  numerous  altitude-dependent  atmospheric 
properties  which  are  examined  in  the  1962  Standard. 
The  variation  of  these  properties  with  altitude  is 
illustrated  in  Figure  2-1  (scaled  in  metric  units). 

2-2.1.  Acceleration  of  Gravity 

The  variation  of  acceleration  of  gravity  g  with 
actual  geometric  altitude  Z  is  classically  expressed 
by  the  equation : 


where  g0  is  a  standard  value  at  zero  altitude  (Table 
2-2),  and  r  is  an  effectiveradius  cf  the  earth  adjusted 
to  provide  for  its  rotation  and  nonspherical  shape 
(Table  2-3).  This  relation  may  be  applied  to  other 
planets  if  suitable  values  of  g0  and  r  are  employed. 
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2-2.2.  Modified  Geopotential  Altitude  H*  Z 

*  r  ~  T  +  Z 

The  modified  geopotential  altitude  H  is  simply 

a  transformed  altitude  scale  with  the  units  of  length.  The  equation  cf  static  equilibrium  of  a  fluid, 

defined  to  facilitate  the  solution  of  the  basic  equation  sometimes  termed  the  "barometric  equation,”  is 

of  static  equilibrium  of  fluids  in  the  atmosphere  by 

absorbing  the  effect  of  variations  in  g.  It  is  defined  ^P  =  —  PO  dZ 

by  the  differential  relation  Substitution  of  Eq.  2-2  into  Eq.  2-4  produces 

gadH*  _  g(Z)  dZ  (2-2)  an  equation  resembling  that  which  holds  for  g  =  g0  = 

.  .  .  .  „  „  ,  .  ,  constant,  namely 

which,  when  combmed  with  hq.  2-1  and  integrated, 

yields  dP  =  -  pg0  till*  (2-5) 

TABLE  2-1. 

LIST  OF  SYMBOLS  AND  ABBREVIATIONS  USED  IN 

STANDARD  ATMOSPHERE  TABLES* 

a 

Acceleration 

lbf 

Pound  (force) 

a 

Radius  cf  the  earth  at  the  equator 

M 

Mean  molecular  weight  cf  air 

b 

Subscript  indicating  base  or  reference  level 

m 

Meter 

BTU 

British  Thermal  Unit 

in' 

Standard  geopotential  meter 

“C 

Degrees,  in  thermodynamic  Celsius  scale 

m 

Mass 

c. 

Speed  cf  sound 

mb 

Millibar 

cal 

Calorie 

mks 

Meter-kilogram-second  system  cf  units 

cm 

Centimeter 

N 

Avogadro's  number 

E 

Energy 

n 

Number  density 

“F 

Degrees,  in  thermodynamic  Fahrenheit  scale 

», 

Loschmid’s  number 

/ 

Ellipsoid  flattening 

nt 

Newton 

F 

Force 

0 

Subscript  indicating  sea-level  value 

fps 

Foot-pound-second  system  cf  units 

p 

Pressure 

ft 

Foot 

pdl 

Poundal 

ft’ 

Standard  geopotential  foot 

°R 

Degrees,  in  thermodynamic  Rankine  scale 

G 

Dimensional  constant  in  gcometric-geopotential  re¬ 
lationship 

R* 

r 

Universal  gas  constant 

Effective  radius  of  earth 

9 

Acceleration  of  gravity,  effective  value 

s 

Sutherland’s  constant 

9c 

Acceleration  of  gravity  at  the  equator 

sec 

Second 

gm 

Gram 

T 

Temperature  in  absolute  thermodynamic  scales 

gm-mol 

H 

Gram  mole 

Altitude  in  geopotential  measure 

Ti 

Ice  point  temperature  in  absolute  thermodynamic 
scales 

Hg 

H, 

Mercury 

Scale  height 

tm 

Molecular-scale  temperature  in  absolute  thermo¬ 
dynamic  scales 

H', 

Geopotential  scale  height 

t 

time 

H* 

Modified  geopotential  altitude 

t 

Temperature  in  nonabsolute  thermodynamic  scales 

i 

Subscript  indicating  ice  point  value 

thsd  ft 

Thousands  of  feet 

in. 

i.  n.  mi 

Inch 

International  nautical  mile 

ti 

Ice  point  temperatnre  in  nonabsolute  thermo¬ 
dynamic  scales 

“K 

Degrees,  in  thermodynamic  Kelvin  scale 

V 

Particle  speed  (arithmetic  average) 

k 

kg-cal 

Thermal  conductivity 

Kilogram-calorie 

V 

Mole  volume  cf  air  under  existing  conditions  cf  T 
and  P 

kg 

Kilogram  (mass) 

z 

Altitude  in  geometric  measure 

kgf 

Kilogram  (force) 

ft 

Constant  used  in  Sutherland’s  viscosity  equation 

kg-mol 

Kilogram-mole 

Y 

Ratio  cf  specific  heats  (Cp/C,.) 

km 

Kilometer 

V 

Kinematic  viscosity 

km’ 

Standard  geopotential  kilometer 

Coefficient  cf  viscosity 

kwhr 

Kilowatt  hour 

V 

Collision  frequency 

L 

Mean  free  path 

P 

Mass  density 

LjV 

Molecular-scale-temperature  gradient 

a 

Effective  collision  diameter  cf  a  mean  air  molecule 

l 

Length 

+ 

Latitude  of  the  earth 

lb 

Pound  (mass) 

w 

Specific  weight 

For  numerical  values  of  physical  constants  see  Table  2-3. 
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TABLE  2-2.  SEA-LEVEL  VALUES  OF  ATMOSPHERIC  PROPERTIES20 


Symbol 

Metric  units  ( mks ) 

English  units  ( fi-lb-sec ) 

Po 

1.013250X  105  newtons  m-2 

2116.22  Ibf  ft-2 

Po 

1 .2250  kg  in-* 

0.076474  lb  ft-3 

to 

15’  C 

59.0'  F 

ffo 

9.80665  m  sec 

32.1741  ft  sec  2 

s 

110.4°  K 

198.72'  R 

Ti 

273.15'  K 

491.67'  R 

(3 

1.458X10-8  kg  sec-1tn_1(°K)_1/2 

7.3025X10-’  lb  ft->sec  i(°R)-”2 

Y 

1.40  (dimensionless) 

1.40  (dimensionless) 

<T 

3.65X10-10  m 

1.1975X10"4  ft 

N 

6.02257 X  1028  (kg-mol)-‘ 

2.73179  X  1028  (Ib-mol)-' 

R* 

8.31432  joules  (°K)_1  mol-1 

1545.31  ft  lb  (lb-mol)-1  (°R)-> 

Differences  between  H*  and  Z  are  generally 
rather  small  except  at  very  high  altitudes. 

2-2.3.  Geopotential  Altitude 

The  geopotential  altitude  H  is  numerically  equal 
to  H*,  but  is  expressed  in  units  of  energy  per  unit 
mass,  or  velocity  squared,  rather  than  in  units  of 
simple  length.  It  is  given  by  the  relation 

H  =  (g0/G)II*  (2-6) 

where  g0  is  defined  in  Table  2-3  and  G  has  the 
value  9.80665  m2  sec-2  m'~x  in  the  metric  system 
and  32.1741  ft2  sec-2  ft'"1  in  the  English  system 
of  units.  The  units  of  H  are  geopotential-mcters  m\ 
where  m'  is  9.80665  (m/sec)2.  Physically,  the  geo¬ 


potential  meter  expresses  the  fact  that  if  1  kg  of 
mass  is  raised  through  one  m',  its  potential  energy 
is  increased  by  9.80665  joules  (see  Table  2-4  for 
conversions  to  English  units). 

2-2.4.  Molecular-Scale  Temperature 

The  molecular-scale  temperature  T  M  is  defined 
so  as  to  simplify  the  equation  of  state  of  a  perfect 
gas,  which  is  assumed  to  apply  throughout  the 
atmosphere,  by  weighting  the  effects  of  changes  in 
average  molecular  weight  at  the  higher  altitudes. 
T M  may  be  related  uniquely  to  the  density,  pressure 
or  speed  cf  sound  without  specification  cf  the  mo¬ 
lecular  weight.  The  density  and  sound  speed  are 
quantities  often  measured  experimentally  in  high 
altitude  tests. 


TABLE  2-3.  NUMERICAL  VALUES  USED  IN  1962  U.  S.  STANDARD  ATMOSPHERE20 


Symbol 

Metric  units  (mks) 

English  units  (ft-lb-sec) 

C..0 

340.294  m  sec-' 

1116.45  ft  sec-1 

*9o 

9.80665  nt  sec-2 

32.1741  ft  sec-2 

flp.o 

8434.5  m 

27,672  ft 

k0 

6.0530X  10_#  kg-cal  m-1  sec_1(°K)-1 

4,0674X10-'  BTU  ft"1  sec-‘(°R)-1 

6.6328  X10-*  m 

2.1761  X10-’  ft 

Mo 

28.9644  (dimensionless) 

28.9644  (dimensionless) 

71q 

2.5471X  1025  m-3 

7.2  1  27  X  1023  ft-3 

*po 

1.013250X  10s  newtons  m-2 

2116.22  Ibf  ft-3 

*TO 

288.15°  K 

518.67°  R 

VO 

458.94  in  sec-1 

1505.7  ft  sec-1 

Vo 

1.4607  XI 0-8  m2  sec-1 

1.5723  X 10-1  ft2  sec-1 

MO 

1.7894X10-5  kg  m_1  sec-1 

1. 2024X10-'  lb  ft"1  sec-1 

*>0 

6.9193 X10»  sec-1 

6.9193  X  10s  sec-' 

*PO 

1 .2250  kg  m-! 

0.076474  lb  ft"3 

*<tq 

3.65  X  10- 10  m 

11.975X10-“  ft 

Wo 

12.013  kg  m-2  sec-2 

2.4605  lb  ft-2  sec-2 

•These  values  are  adopted  for  purposes  of  computation.  The  rcmainingvalues  are  derived  from  the  adopted  values. 
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TABLE  2-4.  CONVERSION  FACTORS' 


(A)  Metric  to  English  Conversions  cf  Units 
cf  Length,  Mass  and  Geopotential 


(D)  Thermal  to  Mechanical  Units, 
Metric-English 


A  Defined  relations  (the  constants  are  adopted  as  being 


exact) : 

1  foot 
li.n,  mi 
1  pound 

1  standard  geopotential 
foot 

B.  Derived  relations: 

1  meter 
1  meter 
1  kilogram 
li.n.  mi 
1  foot 

1  standard  geopotential 
meter 


=  0.3048  meter 
=  1.852  meters 
=  0.45359237  kilogram 
=0.3048  standard  geopo¬ 
tential  meter 

=  3.2808399...  feet 
=  5.3996680...  X10-M.n.  mi 
=  2.2048226...  pounds 
=6,076.1155...  feet 
=  1.6457883...  X10"<  i.n.  mi 
=  3.2808399.  .  .  standard  geo¬ 
potential 
feet 


(B)  Metric  to  English  and  Absolute  to 
Non-Absotute  Conversions  cf 
Temperature  Units 


A.  Defined: 

<(°C)  =T(°K)  — Ti(°K),  where  T,(°K)  = 

273.15(°K) 

=1.8T(°K) 

t(°F)  — i,-(°F)  =  T(°R)-T,(,,R),  where  <i(°F)=32(°F) 

B.  Derived  relations: 

h(°  C)  =0(°C) 

T,(°R)  =491.670(°R) 

<(4C)  =(T(°R)-  Tl(° R)[/ 1 .8  =  ff(°F)— fi(° F)]/ 1.8 

T(°R)  =1.8[t(°C)  +  273.15(°C)]=<(°F)  -1,(°F)  + 

491.670(°R) 

t(°F)  — 32(°F)  =  1.8l(°C)  =  1.8[T(°K)  -  273.15(°K)] 


A.  Defined  relations: 


1  kg-cal  =  gijj  kw-hr  (exact) 


1  kg-cal  =  0  45^89237BTU  =  3.9683207  BTU 

1  joule  =  1  watt-sec 
B.  Derived  relations: 

1  kw-hr=3.6X  10*  watt  sec  =  3.6 X  10*  joules 

3  6V  |08 

1  kg-cal  =  -^-  joules  =  4, 186.0465  joules 
=  4,186.0465  kg  m2  sec-2 

Q  V  1  f|® 

1  ks~cal  =860 x 9.80665™  gf  =  426 .85795  m  kgf 

_ _ ft  thf 

1  kg-cal  =  860X9. 80665X0. 45359237X0.3048 
=3087.4696  ft  Ibf 

0.45359237  ...  ,  , 

1  BTU  = - — ~ - kg-cal  =  0. 26199676  kg-cal 

1 . 0 

1  RTT ’  —  3.6 X  108  .. - 

860  X  0.3048  X  9.80665X  1.8 

=778.02922  ft  lbf 

- ^ _ lb  ft>  sec-2 

1  BTU  =860X  (0.3048)2  X  1.8 

=  25032.349  lb  ft2  see’2 

“Thccaloric  used  herds  Ihc  International  Steam  Table  caloric  and  the  joule 
Is  ihc  mean  international  joule. 


( C)  Absolute  Systems  cf  Units  to 
Absolute-Force,  Gravitational  System  cf 
Units,  Metric-English 

A.  Defined: 

lforceunit  =  l  mass  unit  X  g0 

B.  Derived  relations: 

1  kgf  =9.80666  kg  m  sec"2 

1  kg  =g-g^jkgf  sec2  m"‘  =  0.10197162  kgf  sec2  m"1 

1  lbf  =0.45359237  kgf 
1  lbf  =32.174049  lb  ft  sec-2 
1  lb  =0.031080950  lbf  sec2  ft"1 
=  0.031080950  slug 
1  slug  =  32. 174049  lb 


The  molecular-scale  temperature  is  defined  by 
the  relation 


P/p  =  (R*/M)T 

-  (R  /M0)Tm  p-  m  -  ^  (2-8) 

Thus  Tm  is  seen  to  be  directly  proportional  to  p/p 
which  will  be  seen  to  be  proportional  to  the  square 
cf  the  sound  speed  C,. 

It'  it  is  assumed  that  T M  is  a  function  cf  H* 
which  can  be  estimated  or  measured,  Eqs.  2-5  and 
2-8  can  be  combined  to  yield  the  expression 


TV  =  T_ 
Mo  M 


(2-7) 


where  M0  is  a  standard  molecular  weight  defined 
in  Table  2-3. 

The  equation  of  state  is  expressed  as  follows: 


where 


(2.10) 


In  the  1959  ARDC  Standard,  the  ratio  Tu/T0 
is  the  basic  function  which  is  assumed  as  a  function 
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TABLE  2-4.  CONVERSION  FACTORS20  (continued) 

(E)  Defining  Properties  of  the 
Standard  Atmosphere 


Altitude, 

H,km 

Molecular 
scale  tem¬ 
perature, 
Tm,  "K 

gradient, 
L'  y, 
°K/km 

Molecular 

weight. 

M 

Kinetic 

temper¬ 

ature, 

T,  °K 

0.000 

288.  15 

-6.5 

28.9644 

288.  15 

11.000 

216.65 

0.  0 

28.9644 

216.  65 

20.000 

216.  65 

+1.  0 

28.9644 

216.  65 

32.  000 

228.65 

+  2.  8 

28.9644 

228.65 

47.000 

270.65 

0.  0 

28.9644 

270.65 

52.  000 

270.  65 

-2.0 

28.9644 

270.  65 

61.000 

252.65 

-4.  0 

28.9644 

252.65 

79.000 

180.65 

0.  0 

28.9644 

180.  65 

88.  743 

180.65 

28.9644 

180.65 

Altitude. 

Z,  km 

vtoleeular- 
cale  tem¬ 
perature, 
Ty,  °K 

iradient, 

Ly, 

°K/km 

Molecular 

weight. 

Kinetic 

temper¬ 

ature, 

T,°K 

90 

180.65 

+  3 

28.  9644 

180.65 

100 

210.65 

+  5 

28.  88 

210.02 

110 

260.  65 

+  10 

28.  56 

257.  00 

120 

360.  65 

+20 

28.  07 

349.  49 

150 

960.65 

+  15 

26.92 

892.  79 

160 

1.110.65 

+  10 

26.  66 

022.  2 

170 

1,210.  65 

+  7 

26.40 

l,  103.4 

190 

1,350.  65 

+  5 

25.  85 

l,  205.  4 

230 

1,550.  65 

+  4 

24.  70 

1,  322.  3 

300 

1,830.  65 

+3.  3 

22.  66 

1.432.  1 

400 

2,  160.  65 

+2.  6 

19.94 

l,  487.  4 

500 

2,  420.  65 

+  1.  7 

17.94 

1,  499.  2 

600 

2,  590.  65 

+  1.  1 

16.84 

l,  506.  1 

700 

2,  700.  65 

16.  17 

1,  507.  6 

oiH*  m  order  to  match  test  data.  ML  other  param¬ 
eter  variations  follow  as  outlined  above. 

2-2.5.  Scale  Height 

Scale  heights  are  indicative  of  the  magnitude 
of  vertical  pressure  gradients  in  the  atmosphere  and 
may  be  defined  in  geometric  units  (H,,  meters)  and 
in  modified  geopotential  units  ( H *,  meters)  as 
follows : 


where 

i_=-fzV«P)  and  (2-12) 

2-2.6.  Speed  cf  Sound 

The  speed  of  sound  C.  is  a  scalar  parameter 
defined  here  formally  as  follows: 


where  y  is  the  iscntropic  exponent  or,  in  perfect 
gases,  the  ratio  of  specific  heats.  Actually  y  is  a 
function  of  the  state  and  molecular  structure  cf  the 
gas.  Here  it  is  assumed  that  y  =  1 .4  (exactly  for 
altitudes  below  90  km).  Above  this  altitude,  in¬ 
formation  on  the  sound  speed  is  given  by  neither 
the  1962  U.S.  nor  the  1959  ARDC  Standards. 

The  value  cf  C,  is  not  necessarily  the  actual 
speed  with  which  sound  is  transmitted  over  long 
distances,  although  in  usual  cases  cf  high  density 
gas  and  moderate  wave  frequency,  the  actual  long 
range  propagation  speed  coincides  with  the  value 
of  C,.  In  fact,  the  definition  used  here  (Equation  2- 13] 
stems  from  this  classical  case.  However,  a  real  dis¬ 
turbance  traveling  an  appreciable  distance  from  its 
generator  is  attenuated  or  damped  to  a  degree  which 
increases  with  increasing  wave  frequency  and  w  ith 
decreasing  density'  cf  the  medium.  In  moving  fluids, 
particularly  those  in  which  chemical  reactions  are 
occurring,  still  other  factors  influence  the  actual 
propagation  cf  disturbances.  For  this  reason  the 
parameter  C  loses  its  significance  at  very  high 
altitudes  and  is  replaced  for  reference  purposes  by 
the  mean  speed  cf  the  gas  particles,  a  quantity 
closely  related  to  the  formal  definition  of  C,.  It  is 
noted  that  even  in  the  formal  definition  of  C,  given 
by  Ref.  13,  the  value  cf  y  is  no  longer  considered 
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TABLE  2-5. 

ALTITUDE-DEPENDENT  ATMOSPHERIC  PROPERTIES 


Kinetic  temperature 
Molecular-scale  temperature 
Mean  molecular  weight 
Pressure 
Mass-density 
Specific  weight 
Acceleration  of  gravity 
Mole  volume 


Scale  height 
Number  density 
Mean  free  path 
Collision  frequency 
Speed  of  sound 
Thermal  conductivity 
Kinematic  viscosity 
Viscosity  coefficient  (absolute) 


constant  at  extreme  altitudes  due  to  the  changing 
physical  and  chemical  structure  of  the  atmosphere. 


2-2.7.  Air  Particle  Speed 


The  mean  air  particle  speed  V  is  the  arithmetic 
average  of  the  distribution  of  speeds  of  all  air 
particles  within  a  unit  volume,  assuming  that  all 
air  particles  in  the  mixture  have  an  average  mass 
associated  with  the  mean  molecular  weight.  The 
quantity  has  statistical  meaning  only  if  a  sufficiently 
large  number  cf  particles  is  present  so  that  their 
velocities  follow  a  “Maxwellian”  statistical  dis¬ 
tribution5  and  provided  that  variations  in  P  and 
p  are  negligible  within  the  volume.  The  following 
definition,  derived  from  the  classical  kinetic  theory 
of  gases5,  is  introduced: 


V  = 


This  yields 


where 


(2-14) 


(2-15) 


Although  the  physical  meaning  of  V  also  de¬ 
creases  at  extremely  high  altitudes,  it  is  utilized  on 
the  basis  of  its  formal  definition. 


2-2.8.  Specific  Weight 

The  specific  weight  o>  is  defined  as  follows: 

u  =  pg  (2-16) 

This  simply  Newton's  second  law  referred  to 
a  unit  volume,  where  the  mass  density  is  assumed 
to  be  defined  and  the  force  (weight)  derived. 

2-2.9.  Molecular  Weight 

The  molecular  weight  M  of  a  compound  is 
defined  from  a  chemical  standpoint  as  a  dimension¬ 


less  number  whose  magnitude  is  16  times  the  ratio 
of  the  average  mass  of  a  molecule  of  the  compound 
to  the  average  mass  of  an  oxygen  atom.  It  is  assumed 
that  the  oxygen  and  the  compound  have  their 
natural  distribution  of  isotopes  and  that  the  average 
is  to  be  construed  as  the  arithmetic  mean. 

A  unit  such  as  the  gram-mole  is  defined  as  the 
quantity  of  mass  in  grams  which  is  numerically  equal 
to  the  molecular  weight  of  the  substance.  That  is, 
1  gram-mole  =  M  gm.  The  molar  concentration  is 
defined  as  the  number  of  moles  of  a  component 
(in  metric  units,  gram-mole),  per  unit  volume  of 
gas  mixture.  A  mole-fraction  is  defined  as  the  number 
cf  moles  of  a  component  per  mole  of  mixture. 

Standard  air  is  assumed  to  be  dry.  Its  standard 
sea-level  molar  composition  is  shown  in  Table  2-6. 
This  composition  is  assumed  to  be  constant  below 
90  km.  Therefore  the  following  standard  molecular 
weight  for  the  air  mixture  applies  in  this  regime: 
Mi,  =  28.966.  For  higher  altitudes  the  mean  molec¬ 
ular  weight  variation  shown  in  Fig.  2-1  (F)  is  as¬ 
sumed  as  a  standard,  based  upon  experimental 
evidence. 

2-2.10.  Mole  Volume 

The  mole  volume  v  is  defined  by  the  relation 
v  =  M/p  (volume/mole) 

R*T 

=  - -  for  perfect  gases  (2-17) 


TABLE  2-6. 


SEA  LEVEL  ATMOSPHERIC  COMPOSITION 
FOR  A  DRY  ATMOSPHERE* 


Constituent  Gas 

Mol.  Fraction 
(Percent) 

Molecular  Weight 
(0  =  16.000) 

Nitrogen  (N2) 

78.09 

28.016 

Oxygen  (02) 

20.95 

32.0000 

Argon  (A) 

0.93 

39.944 

Carbon  dioxide  (C02) 

0.03 

44.010 

Neon  (Ne) 

1.8  X  10-3 

20.183 

Helium  (He) 

5.24  X  10-‘ 

4.003 

Krypton  (Kr) 

1.0  x  10-4 

83.7 

Hydrogen  (H2) 

5.0  X  lO”6 

2.0160 

Xenon  (Xe) 

8.0  X  10"6 

131.3 

Ozone  (03) 

1.0  X  10-« 

48.0000 

Radon  (Rn) 

6.0  X  10“la 

222.0 

*  These  values  are  taken  as  standard  and  do  not  neces¬ 
sarily  indicate  the  exact  condition  cf  the  atmosphere.  Ozone 
and  Radon  particularly  are  known  to  vary  at  sea  level  and 
above,  but  these  variations  would  not  appreciably  affect  the 
value  cf  M 0. 
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At  given  values  of  P  and  T,the  mole  volume  is  the 
same  for  all  perfect  gases. 

2-2.11.  Number  Density 

The  number  density  n  of  a  gas  denotes  the 
number  of  gas  particles  per  unit  volume,  neglecting 
electrons  or  other  sub-atomic  particles.  Since 
Avogadro's  number  N,  defined  as  the  number  of 
particles  in  one  mole  cf  gas,  is  a  universal  constant, 
it  follows  that 

N 

n  =  -y  particles/vol  (2-18) 

where  N  =  6.024  X  1026  particles/kg-mole  of  gas 

2-2.12.  Collision  Frequency 

By  classical  methods  of  the  kinetic  theory  of 
gases6  an  estimate  may  be  made  of  the  frequency  v 
with  which  a  given  gas  particle  will  collide  succes¬ 
sively  with  other  particles.  In  the  simplest  case, 
all  participating  particles  are  characterized  by  the 
same  dimension,  a  diameter  a,  which  expresses  the 
effective  extent  of  a  particle's  force  field.  Moreover, 
the  gas  is  considered  to  behave  as  a  perfect  gas  w  ith 
cr  small  compared  to  distances  between  particles. 
The  result  is  the  follow  ing  expression  for  the  collision 
frequency  : 

v=V2m*V  (2-19) 

=  collisions/sec,  or  particles  colliding 
with  a  given  particle/sec 

Although  care  must  be  used  in  ascribing  too 
much  significance  to  the  dimension  a,  because  of 

the  actual  field-like  nature  of  a  particle's  range  of 
influence,  the  expression  for  v  is  carried  over  formally 
to  the  standard  air  mixture.  For  this  purpose  an 
average  effective  collision  diameter  is  defined  as 
follows : 

cr  =  3.65  x  l(r10  meters  =  1.199  X  10~9  ft  (2-20) 

2-2.13.  Mean  Free  Path 

A  mean  free  path  length  L  which  would  be 
traversed  by  a  particle  traveling  with  mean  speed  V 
and  participating  in  v  collisions  per  second,  is  given 
by  the  expression 

L  =  —  =  (V2W)-1  (2-21) 


cules  and  atoms  comprising  the  atmosphere  at  a 
given  altitude.  Usually  the  random  translational 
degrees  of  freedom  furnish  the  kinetic  energy  by 
which  the  temperature  is  defined.  For  a  homogeneous 
gas  with  a  "Maxwellian"  particle  velocity  distribu¬ 
tion  the  previously  given  expression  for  the  mean 
random  speed14  yields  a  relation  between  the  random 
kinetic  energy  and  the  temperature,  namely, 


|  MV2  =  -R*T  (2-22) 

or  in  terms  of  the  mean  square  speed  V2, 

\MV2  =  |  R*T  (2.23) 

The  kinetic  temperature  T  is  not  necessarily 
the  temperature  which  would  be  acquired  by  a  body 
at  rest  within  the  gas  or  moving  through  it,  since 
the  ultimate  temperature  of  a  body  depends  upon 
a  balance  between  various  modes  of  energy  transfer. 

2-2.15.  Viscosity 


The  absolute  viscosity  m  and  the  kinematic 
viscosity  n  =  p/p  are  obtained  by  use  of  the  Suther¬ 
land  equation,  which  is  based  on  the  kinetic  theory 
cf  gases  and  is  adjusted  empirically,  as  follows: 


nm3/2 

M  =  T  +  S 


(2-24) 


where 


(3  —  1.458  X  10-6  kg/sec 'in  -  \/°K  and 
s  =  110.4°K. 


The  use  cf  this  relation  should  be  limited  to  altitudes 

below  40  km. 

2-2.16.  Thermal  Conductivity 


A  modified  Sutherland  type  cf  relation  is  used 
for  the  thermal  conductivity  k,  namely, 


k  =  T  +  ,S\  X  ft)-12/r  kg'cal/m'sec  °K  (2’25> 
where 


&  =  6.235  X  10"7 
Sl  =  245. 4°K 

and  the  relation  is  also  limited  to  altitudes  below 
40  km. 


2-2.14.  Temperature  (Real  Kinetic) 

The  temperature  T  is  a  quantity  related  to  the 
kinetic  energy  per  unit  mass  of  the  random  motion 
(as  distinct  from  ordered  mass  motion)  of  gas  mole¬ 


2-2.1 7.  Simplified  Exponential  Approximation  of 
Density  Variation 

Both  the  pressure  and  density  vary  approx¬ 
imately  exponentially  with  altitude  in  the  atmos- 
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phere,  each  decreasing  by  a  factor  cf  about  10 
roughly  every  50,000  ft.  In  trajectory  analyses  it 
has  been  found  useful  to  approximate  the  1956 
ARDC  Standard  density  as  follows6,7  for  altitudes 
below  150  km  (490,  000  ft.) : 

=  exp  (-M 

Pr 

where  pT  is  a  reference  density  evaluated  to  achieve 
the  best  fit  cf  the  actual  density  in  the  altitude 
range  cf  prune  interest,  and  dr  is  a  local  decay  factor. 
A  fit  cf  the  1956  Standard  is  show  n  in  Figs.  2-2  and 
2-3  where  the  following  values  have  been  used  : 

pr  =  0.0027  slugs/ft3 
1  =  23,500  ft 

Pr 

0  r(r  +  Z)  =  900 

A  series  cf  local  values  of  0,  estimated'  in  the  70 
to  169  km  range  (230,000  to  540,000  ft)  are  also 


shown  marked  on  Fig.  2-2.  To  obtain  these  values, 
values  of  pr  were  selected  at  the  altitudes  indicated, 
and  corresponding  values  of  l/0r  varying  between 
18.500  and  29,800  ft  were  derived  for  the  intervals 
delineated. 

2-2.18.  Extreme  Deviations  from 
Standard  Conditions 

Deviations  from  standard  temperatures  are 
found  to  approach  practical  extremes  under  the 
following  general  circumstances:'-  arctic  winter, 
arctic  summer  and  tropical  day.  Values  of  these 
temperature  deviations,  compared  with  the  ARDC 
Standard  for  altitudes  up  to  100,000  ft  are  shown 
in  Figs.  2-4  and  2-5. 

The  corresponding  pressures  and  densities  are 
found  by  solution  cf  the  hydrostatic  equations  for 
perfect  gases,  as  previously  discussed  in  this  hand¬ 
book,  assuming  an  invariant  atmospheric  composi¬ 
tion.  Values  of  these  extreme  conditions  are  tabulated 
in  Ref.  1  and  are  also  discussed  in  Refs.  9  and  10. 


Figure  2-2.  Comparison  cf  Exponential  Approximation  with  ARDC  Model  cf  Earth  Atmosphere 
(1956).  (Solid  circles  indicate  1962  U.  S.  Standard  values.) 
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2-3  PROPERTIES  OF  TVE  ATMOSPHERE 
IMPORTANT  FOR  ELECTROMAGNETIC 
WAVE  TRANSMISSION 

2-3.1.  Index  of  Refradion  ( D ) 


The  index  cf  refraction  cf  air  for  radio  waves 
is  defined  as  the  ratio  of  the  propagation  velocity 
in  real  air  to  the  propagation  velocity  which  would 
be  attained  hi  a  vacuum,  that  is,  the  speed  of  light. 
This  index  is  essentially  independent  of  frequency, 
except  possibly  for  wavelengths  of  less  than  one  cm, 
and  is  also  independent  of  polarization.  It  is  primarily 
a  function  cf  the  air  temperature,  pressure  and 
moisture  content”  as  mdicated  by  the  following 
empirical  formula : 


(D  -  1)  X  108 


79 P 
T 


11 A 
T 


3.8  X  106A 


where 


Figure  2-3.  Dimensionless  Parameter  s/jf  for  P  =  barometric  pressure  (millibars) 

ARDC  Model  of  Earth  Atmosphere”  r„  /OT,. 

J  1  T  =  temperature  (°K) 


2-2.19.  Atmospheric  Shells 


A  =  water  vapor  pressure  (millibars) 


A  number  cf  classifications  cf  altitude  ranges 
in  the  atmosphere  have  been  proposed.  Several  are 
indicated  in  Figures  2-6(A)  and  (B)." 


Figure  2-4.  Temperature-Height  Profiles  of  Extreme 
and  ARDC  Model  Atmospheres  to  100,000ft 1 


For  many  purposes  it  is  sufficiently  accurate  to 
assume  that  D  decreases  linearly  with  altitude. 


TFHPERATURC  *  KILVIN 
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Figure  2-5.  Hot  and  Cold  Temperature  Extremes 
and  ARDC  Model  Atmosphere  Temperature- 
Height  Profile  to  100, 000 f  ( 
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(A)  Schematic  Representation  of  Proposed  Atmospheric  Nomenclature 
Figure  2-6.  Systems  of  Nomenclature' 


2-3.2.  Electrical  Conductivity  (a) 

The  current  in  a  gaseous  medium  can  be  shown 
to  be  given  by 


where 


I  =  current  (amps) 

N  =  number  of  electrons/cm3 
q  =  electric  unit  charge  (esu  in  cgs  electrostatic 
units) 

v  —  velocity  of  a  charged  particle  (em/see) 

G  =  average  frictional  constant 

i  =  V--1 

w  =  angular  frequency  of  electric  field  (radian/ 

sec) 

m  =  mass  of  charged  particle  (gm) 

E„,  =  amplitude  of  applied  electric  field  (stat- 
volts/cm) 

/  =  time  (sec) 


Only  the  current  due  to  electrons  need  be 
considered  because  of  the  large  mass  of  other 
particles. 


The  conductivity  is  then  a  complex  number 
given  as  follows: 


where 

r  s  m/G,  a  characteristic  time,  for  practical  pur¬ 
poses  the  mean  free  time  between  electron 
collisions. 

Thus  it  is  seen  that  the  population  cf  electrons 
N  in  the  air.  primarily  in  the  ionosphere,  is  a  quantity 
which  is  significant  for  wave  transmission. 

2-3.3.  Electron  Density  in  the  Atmosphere  (Ar) 

The  levels  of  maximum  electron  density  at 
middle  latitudes13  are  given  in  Table  2-7  and  Fig.  2-7. 
Figure  2-8(B)  illustrates  the  estimated  electron  dis¬ 
tribution  in  the  Van  Allen  belts  located  at  altitudes 
of  several  earth  radii.14 
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It  is  of  interest  to  note  that  the  electron  density 
and  resulting  conductive  properties  of  air  at  very 
high  temperatures  have  become  cf  interest  because 
cf  the  high  temperatuers  generated  in  regions  very 
near  the  surfaces  of  high-speed  aircraft.  Since  the 
transmission  of  electromagnetic  waves  through  this 
hot  gas  blanket  may  be  necessary  ,  the  properties 
which  govern  the  transmission  characteristics,  pri¬ 


marily  N  and  r,  must  be  obtained  for  the  thermo¬ 
dynamic  state  of  the  hot  gas  blanket.  Once  these 
properties  are  evaluated,  the  conductivities  and 
other  characteristics  may  be  estimated  as  for  any 
other  transmitting  medium.  This  also  holds  true  for 
transmission  through  flaming  exhaust  gases.  A 
further  discussion  of  these  problems  will  be  given 
in  a  later  paragraph. 


NAME 

DESCRIPTION 

TEMPERATURE 

TROPOSPHERE 

The  region  nearest  the  surface,  having  a  more  orless  uniform  decrease  of  tempera¬ 
ture  with  altitude.  The  nominal  rate  of  temperature  decrease  Is  6.  5°K/km,  but  in¬ 
versions  are  common.  'The  troposphere,  the  domain  of  weather,  is  in  convective 
equilibrium  with  the  sun- warmed  surface  of  the  earth.  Thelropopause,  which  occurs 
at  altitudes  between  6  and  18  kilometers  (higher  and  colder  over  the  equator),  is  the 
domain  cf  high  winds  and  highest  cirrus  clouds. 

STRATOSPHERE 

The  region  next  above  the  troposphere  and  having  a  nominally  constant  temperature. 

The  stratosphere  is  thicker  over  the  poles,  thinner  or  even  nonexistent  over  the 
equator.  Maximum  cf  atmospheric  ozone  found  near  stratopause.  Rare  nacreous 
clouds  also  found  near  stratopause.  Stratopause  is  at  about  25  kilometers  in  middle 
latitudes.  Stratospheric  temperatures  areinlhe  order  of  arctic  winter  temperatures. 

MESOSPHERE 

The  region  of  the  first  temperature  maximum.  The  mesosphere  lies  above  the  stra¬ 
tosphere  and  below  the  major  temperature  minimum,  which  la  found  near  80  kilo¬ 
meters  altitude  and  constitutes  the  mesopause.  A  relatively  warm  region  between 
two  cold  regions;  the  region  of  disappearance  of  moel  meteors.  The  mesopause  is 
found  at  altitudes  of  from  70  to  85  kilometers.  Mesosphere  la  in  radiative  equilibrium 
between  ultraviolet  ozone  healing  fcy  the  upper  fringe  of  ozone  region  and  the  infrared 
ozone  and  carbon  dioxide  cooling  by  radiation  to  space. 

THERMOSPHERE 

The  region  of  rising  temperature  above  the  major  temperature  minimum  around  80 
kilometers  altitude.  No  upper  altitude  limit.  The  domain  cf  the  aurorae.  Tempera¬ 
ture  rise  at  base  cf  thermosphere  attributed  to  tew  Infrequent  collisions  among  mole¬ 
cules  to  maintain  thermodynamic  equilibrium.  The  potentially  enormous  infrared 
radiative  cooling  by  carbon  dioxide  hs  not  actually  realized  owing  to  inadequate  colli¬ 
sions. 

COMPOSITION 

HOMOSPHERE 

The  region  cf  substantially  uniform  compoeilion.  in  the  sense  cf  constant  mean 
molecular  weight,  from  the  surface  upwards.  The  homopause  is  found  at  altitudes 
between  80  and  100 kilometers.  The  compoeilion  changes  here  primarily  because  of 
dissociation  cf  oxygen.  Mean  molecular  weight  decreases  accordingly.  The  ozonos- 
phere,  having  its  peak  concentration  near  stratopause  altitude,  does  not  change  the 
mean  molecular  weiaht  of  the  atmosphere  significantly. 

HETEROSPHSRE 

The  region  cf  significantly  vary  ing  composition  above  the  homosphere  and  extending 
indefinitely  outwards.  The  "molecular weight"  of  air  diminishes  from  29  at  about  90 
kilometers  to  16  at  about  500  kilometers.  Well  above  the  level  of  oxygen  dissocia¬ 
tion,  nitrogen  begins  to  dissociate  and  diffusive  separation  (lighter  atoms  and  mole¬ 
cules  rising  to  the  lop)  sets  in. 

IONIZATION 

IONOSPHERE 

The  region  of  sufficiently  large  electron  density  to  affect  radio  communication.  How¬ 
ever.  only  about  one  molecule  in  1000  in  the  F2  region  to  one  in  100,000,000  in  the 

D  region  is  ionized.  The  bottom  of  the  ionosphere,  the  D  region,  is  found  at  about  80 
kilometers  during  the  day.  At  night  the  D  region  disappears  and  the  bottom  of  the 
ionosphere  rises  to  100  kilometers.  The  top  of  the  ionoephere  is  not  well  defined  but 
has  often  been  taken  as  about  400  kilometers.  The  recent  extension  upward  to  1000 
km  based  on  satellite  and  rocket  data  is  shown. 

CHEMICAL  REACTIONS 

CHEMOSPHERE 

The  region  where  chemical  activity  (prtmarLly  photochemical)  is  predominant.  The 
chemosphere  ia  found  within  the  altitude  limits  of  about  20  to  110  kilometers. 

MOLECULAR  ESCAPE 

EXOSPHERE 

The  region  wherein  molecular  escape  from  the  earth's  atmosphere  is  significant. 

The  base  of  the  exosphere,  the  critical  level  is  thought  to  be  at  analttfude  above 

300  kilometers,  possibly  as  high  as  1000  kilometers.  Satellite  data  indicating  higher 
densities  at  these  altitudes  favor  higher  exosphere  levels.  Lighter  atoms  and  mole¬ 
cules  can  escape  at  lower  altitudes  than  heavier  ones.  The  earth's  magnetic  field  ef¬ 
fectively  prevents  the  escape  of  charged  particles,  however. 
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(B)  Salient  Features  cf  Various  Atmospheric  Regions 
Figure  2-6.  Systems  of  Nomenclature'  (continued) 


ELECTRON  NUMBER  DENSITY  (CM1) 

Figure  2-7.  Average  Variation cf  Electron  Densify 
with  Altitude  in  the  I)  Ionic  Layer  During  a 
Summer  Day‘ 


TABLE  2-7. 

LEVELS  OF  MAXIMUM  ELECTRON  DENSITY 
AT  MIDDLE  LATITUDES10 


Approximate  N 
(electrons  per  cm:) 


Region 

Level 

(km) 

Range 

(km) 

(T>ay) 

(Night) 

D 

60 

60-100 

1.5  X  W 

— 

E 

100 

70-120 

1.5  X  10' 

1  X  104 

E, 

120 

75-130 

3.0  X  10' 

— 

Fi 

200 

160-220 

2.5  X  10' 

— 

Ft 

300 

250-500 

1.5X  10' 

2.5  X  10® 

G 

500 

400-700 

2.0  X  10' 

— 

2  4  WINDS 

2-4-1.  Introduction 

Although  the  presence  of  winds  and  gusts  in 
the  atmosphere  must  be  considered  in  formulating 
the  structural  loads  on  a  missile  in  flight,  these 
factors  are  even  more  important  for  determining  the 
characteristics  of  flight  control  systems  and  disper¬ 
sion  in  ballistic  flight.  For  ballistic  vehicles,  control 


IOS  10*  I01 

PROTON  FLUX  (proton* /cm*-  »*c ) 


(A)  Particle  Flux  versus  Altitude  in  the  Equatorial  Plane 
Figure  2-8.  Particle  Flux  in  Van  Allen  Radiation  Belt 21 
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LATITUDE  RELATIVE  TO  INTEGRAL  INVARIANT  EQUATOR 


(B)  Quiet  Day  Flux  Contours  of  Electrons  with  Energies  Greater  than  20  kev  in  Electrons /cm2 -sec 
(Proton  Flux  Contours  are  shown  for  reference.  Accuracy  is  only  to  several  orders  of  magnitude. 
Data  shown  is  as  of  1961.  New  information  is  continually  being  evolved  and  should  be  consulted 
for  quantitative  utility.) 


Figure  2-8.  Particle  Flux  in  Van  Allen  Radiation  Belt 21  (continued) 


systems  are  particularly  important  for  establishing 
accurately  the  terminal  properties  of  the  boost 
stages.  Launching  inaccuracies  cannot  be  corrected 
in  the  ballistic  flight.  Therefore,  design  atmospheres 
for  wind  magnitude  profiles  and  the  rate  of  change 
of  wind  magnitude  with  altitude,  that  is,  “wind 
shear”  distributions,  must  be  established  in  connec¬ 
tion  with  a  given  missile  design. 

Figures  2-9(A)  and  2-9(B)  illustrate  synthetic 
wind  speed  profiles  for  winter  conditions  in  the 
northeastern  and  northwestern  parts  of  the  United 
States. 

The  main  factorswhich  are  influenced  by  winds, 
and  which  are  important  for  determining  control 
system  features  and  ultimate  dispersion  of  ballistic 
stages,  are:  maneuver  loads  due  to  winds,  possible 
dispersion  due  to  winds  and  the  effect  of  wind  shear 
on  a  system  in  which  control  of  the  longitudinal 
Mach  number  involves  employment  of  on-board 
aerodynamic  sensors. 


Wind  effects  are  relatively  unimportant  during 
the  high  speed  flight  of  heavy,  nonmaneuvering 
ballistic  stages  such  as  modern  re-entry  nose  cones 
cf  long  range  missiles.  However  they  may  be  signif¬ 
icant  for  new  maneuverable  re-entry  concepts  and 
for  ballistic  interceptors  which  employ  lift  (glide) 
after  thrust  is  terminated.  These  subjects  are  dis¬ 
cussed  in  detail  in  References  13,  15,  16  and  17. 

2-4.2.  loads  Due  to  Winds  on  a 
Maneuvering  Missile 

A  ballistic  missile  is  controlled  during  its  boost 
and  aiming  stages.  A  completely  guided  missile  is 
controlled  during  the  time  it  closes  in  on  the  target 
as  well.  The  control  load  problem  is  the  same  pro¬ 
vided  the  missile  is  commanded  to  follow  a  pro¬ 
grammed  attitude  (angle  of  attack)  over  a  period 
cf  tune.  The  presence  of  winds  during  this  period 
induces  an  additional  angle  of  attack  over  the  no- 
wind  angle.  This  results  in  an  additional  aerodynamic 
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ALTITUDE  IN  THOUSANDS  OF  FSET 


O  50  100  150  200  250  300 

WIND  SPEED  IN  FT/SEC 

Figure  2-9(A).  Synthetic  Wind  Speed  Profiles  Ex¬ 
ceeded  1%,  5%,  10%,  20%  and  50%  of  the 
Winter  for  the  Windiest  Area  (Northeastern  Part) 
of  the  United  States' 

maneuver  load  dependent  on  both  wind  velocity 
and  shear,  which  the  control  system  must  handle. 

ft'  a  missile  is  aerodynainically  stable,  the  wind 
load  increment  will  set  up  a  moment  which  causes 
the  missile  to  decrease  its  angle  cf  attack  and  thereby 


Figure  2-10.  Rate  of  Heat  Gain  by  Solar  and 
Nocturnal  Radiation 18 


WIND  SPEED  IN  FT/SEC 

Figure  2 -9(B).  Synthetic  Wind  Speed  Profiles  Ex¬ 
ceeded  1%,  5%,  10%,  20%  and'  50%  of  the 
Winter  for  the  Calmest  Area  (Northwestern  Part) 
of  the  United  States' 

reduce  the  load  increment.  However  if  the  missile 
is  changing  altitude,  the  nature  cf  the  wind  shear 
and  the  missile's  time  rate  of  change  cf  altitude 
may  be  such  that  the  missile  response  to  wind  loads 
does  not  lead  to  a  load  reduction.  For  the  Navaho 
SM-64A  missile,  a  guided  robot  plane  with  ramjet 
engines,  solid-fuel  launching  boosters  and  5000-mi 
range,  a  large  portion  of  the  lateral  and  longitudinal 
maneuver  loads  during  boost  and  terminal  dive-in 
were  directly  attributable  to  winds.  Critical  trajec¬ 
tory  conditions  existed  in  the  25,000  to  45,000  ft 
altitude  band,  a  region  of  severe  wind  velocities  and 
wind  shears.  In  this  band,  during  boost,  winds  con¬ 
tributed  nearly  all  the  lateral  load  requirement 
(±90,000  lb)  and  the  largest  dispersion  in  longitu¬ 
dinal  maneuver  loads  (±30.000 lb). 

2-4.3.  Wind  Dispersions 

Winds  contribute  strongly  to  aerodynamic  and 
flight  control  system  dispersions  ;  the  critical  condi¬ 
tion  for  dive  dispersions  is  that  of  maximum  wind 
velocity  during  the  entire  dive.  Mach  number  fluc- 
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tuations  due  to  changing  winds  may  engender  large 
control  load  requirements  in  systems  which  are 
commanded  by  means  of  Mach  number  sensing 
devices.  This  is  true  particularly  if  the  detected 
Mach  number  is  too  low. 

2-4.4.  Control  of  Horizontal  Flight 

If  part  of  a  controlled  trajectory  is  horizontal, 
if  the  control  system  is  actuated  by  a  Mach  number 
sensing  device  and  if  speed  corrections  are  made  by 
climbing  or  diving,  the  control  sy  stem  will  be 
strongly  sensitive  to  wind  shear  which  results  in 
changes  in  relative  velocity.  Consequently,  if  the 


missile  acceleration  capabilities  are  insufficient  to 
decrease  the  indicated  error  signal,  the  missile  sys¬ 
tem  may  become  destabilized. 

2-5  RADIATION  FROM  SOLAR  SYSTEM 
AND  SURROUNDING  ATMOSPHERE 

The  heat  flux  q  in  Btu/ft2  sec  due  to  direct 
solar  irradiation,  due  to  nocturnal  irradiation  from 
below  the  missile  (mainly  caused  by  the  earth)  and 
that  due  to  nocturnal  irradiation  from  above  (mainly 
caused  by  the  gaseous  atmosphere)  are  shown  in 
Fig.  2-10. 
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CHAPTER  3 

TYPICAL  BALLISTIC  TRAJECTORIES 


3-1  INTRODUCTION 


atmospheric  density  by  the  simple  relation 


Although  a  detailed  discussion  of  trajectories 
is  given  in  AMCP  706-284  (C),  Engineering  Design 
Handbook,  Ballistic  Missile  Series,  Trajectories  (U), 
a  number  of  interesting  trajectory  characteristics 
which  can  be  represented  in  relatively  simple  fashion 
are  restated  here  for  convenience.  A  ready  reference 
for  this  simplified  presentation  is  Reference  1  which 
is  based  in  past  on  the  discussions  cf  References 
2  and  3. 


3-2  BALLISTIC  FLIGHT  AND  RE-ENTRY 


The  powered  flight  of  launching  usually  termi¬ 
nates  above  the  effective  atmosphere.  At  the  power- 
off  or  burnout  point,  the  vehicle’s  flight  is  charac¬ 
terized  by  a  velocity  F<  and  an  inclination  angle  0, 
between  the  velocity  vector  and  local  horizontal. 
Subsequent  flight  outside  the  atmosphere  occurs  in 
an  elliptic  trajectory  with  the  earth’s  center  as  one 
focus.  Thus  the  vehicle  again  approaches  the  atmos¬ 
phere  with  velocity  V,  =  Vt  and  an  angle  9,  =  ©,. 

For  minimum  expenditure  of  energy  to  achieve 
a  given  range  (4>  =  range  -t-  earth’s  radius)  the 
injection  and  re-entry  velocities  are  given  by 


V<  =  V,  =  V . 


2  sin  (<t>/2) 

.1  +  sin  (-f>/2)  J 


o  <  $  <  n,  (3-i) 

$  measured  in  radians 


where  V,  is  the  orbiting  speed  near  the  earth’s 
surface,  equal  to  'V gr0  or  25,930  ft/see, 

g  =  32.2  ft/sec2,  and 
n0  =  earth’s  radius  in  feet. 


p/po  =  e  Py  (3-3) 

where 

Po  =  0.0034  slugs/ft3 
1/d  =  22,000  ft 
y  =  altitude  in  ft. 

With  this  relation  and  neglecting  lift,  variations 
in  the  drag  coefficient  CD,  and  gravity  forces  as 
compared  to  aerodynamic  forces,  one  obtains  the 
following  expression  for  velocity  history  during  flight 
within  the  atmosphere  : 

F 

Y  =  exp 

where 

B  is  the  ballistic  coefficient,  W/CDA 
W  is  the  vehicle  weight, 

A  is  the  area  upon  which  CD  is  based,  and 
p  is  the  local  density. 

The  value  of  CD  for  many  blunt  bodies  at  high 
speed  approximates  that  of  a  sphere  or  a  value  on 
the  order  of  0.9.  For  satellite  configurations  and 
weights,  values  of  B  are  on  the  order  cf  10  lb/ft2. 
For  a  representative  body  such  as  an  iron  sphere 
3.6  in.  diameter,  B  =  100 lb/ft2;  for  a  3-ft  diameter 
iron  sphere,  B  =  1000  lb/ft2. 

The  deceleration  of  the  vehicle  measured  in 
units  cf  g  can  be  shown  to  be  equal  to  the  ratio 
of  drag  to  weight.  This  in  turn  can  be  expressed 
by  pV2/2B. 

Examination  of  deceleration,  characteristics 
(Figure  3-1)  shows  the  following: 


-P0 


2/3B  sin  J 


(3-4) 


The  corresponding  launching  and  re-entry  angles  are 


For  a  6000  mile  range  V{  ~  23,000  ft/sec  and 
0;  ~  23°,  whereas  for  a  2000  mile  range  F,  ~ 
16,000  ft/sec  and  0,  ~  38’’. 

In  the  analysis  of  vehicle  re-entry  it  is  con¬ 
venient  to  approximate  the  vertical  variation  of 


a.  Peak  deceleration  occurs  when  V  =  0.61  V{ 
regardless  cf  B. 

b.  Deceleration  is  negligibly  small  at  altitudes 
above  250,000  ft  (less  than  1  g  from  B  >  10) 
unless  one  considers  satellite  flight  over  long 
periods  of  time.  This  indicates  that  rarefied 
gas  effects  are  relatively  unimportant  insofar 
as  aerodynamic  forces  cf  ballistic  vehicles  are 
concerned. 
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A.  MAX.  DECELERATION  AT  V*  0.61  V| 

6.  MAX.  AVERAGE  HEATING  RATE  AT  V-0.72  Vj 

C.  MAX.  STAGNATION  POINT  HEATING  RATE  AT  V-0.85  Vi 


Figure  3-1.  Altitude  and  Velocity  for  Maximum 
Deceleration  and  Heating  cf  Ballistic  Missiles 


The  glide  velocity  history  as  a  function  of  the 
partial  range  T  (angular  displacement  from  launch 
site  swept  out  by  radius  vector  joining  vehicle  and 
earth’s  center)  is  given  by 


—  exp 


—  2(g>  ~  gT 

L  Lift  . 


or,  as  a  function  of  altitude. 

fZY  _  i  +  pVI  . 

\V.)  K 

£.  =  e-ftr 
Po 


(3-5) 


(3-6) 


where 


'  =  J*L. 

CtA  ’ 


CL  is  the  lift  coefficient. 


c.  Maximum  deceleration  cf  58  g  occurs  for  a 
minimum  energy  range  of  4700  mi.  regardless 
cf  B. 

It  can  be  shown  that  within  reasonable  engineering 
accuracy  the  heat  transfer  rate  to  the  stagnation 
point  of  a  blunt  body  is  proportional  to  \/pF3, 
whereas  the  overage  heating  rate  of  the  entire 
surface  of  blunt  or  sharp  vehicles  is  roughly  propor¬ 
tional  to  pV3.  For  a  ballistic  missile,  the  maxima 
of  these  rates  occur  at  F/F,-  =  0.85  (stagnation 
point)  and  F/F,-  =  0.72  (average  heating). 

3-3  LIFTING  RE-ENTRY  (GLIDE) 

Basic  characteristics  of  re-entry  for  unpowered 
vehicles  whose  lift  (force  component  normal  to  the 
flight  path)  cannot  be  neglected  may  be  deduced 
in  a  simplified  manner  by  the  use  of  a  reasonable 
set  of  approximations,  namely  : 

a.  Sum  of  lift  and  centrifugal  force  is  equal  to  the 
vehicle  weight. 

b.  Angle  between  flight  path  and  local  horizontal 
is  small. 

c.  Centrifugal  force  is  approximately  WV2/gr0. 

d.  Lift-to-drag  ratio  L/D  is  constant. 

The  above  assumptions  are  poorest  near  impact. 

Skip  vehicles,  which  enter  and  leave  the  atmos¬ 
phere  again  one  or  more  tunes,  are  not  considered 
here.  Despite  the  fact  that  these  vehicles  are  efficient 
in  the  conversion  cf  velocity  to  range,  they  have 
a  number  of  disadvantages  related  to  repeated  and 
large  decelerations  and  heating  cycles. 


L-  -Cc.  , 

D  ~  C'D  'an<l 

0  <  Ar  <  <f>. 

The  initial  velocity  required  for  a  given  range 
is  shown  plotted  for- values  cf  L/D  in  Figure  3-2, 
based  on  the  relation 


The  desirability  cf  having  large  L/D  is  clearly 
indicated.  For  low  values  of  L/D,  longer  range  is 
achieved  with  a  nonlifting  vehicle  for  a  given  F,-. 


Figure  3-2.  Initial  Velocity  cf  Glide  Vehicles  Re¬ 
quired  to  Achieve  a  Given  Range,  Showing  the 
Effect  cf  the  Lift-to-Drag  Ratio 
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CHAPTER  4 

PHYSICAL  AND  THERMODYNAMIC  PROPERTIES  OF  AIR 


4-1  EQUILIBRIUM  THERMODYNAMIC  AND 
TRANSPORT  PROPERTIES 

In  Chapter  2  the  properties  of  air  are  given 
from  the  standpoint  of  their  variation  with  altitude. 
In  this  chapter  the  equilibrium  thermodynamic  and 
transport  properties  of  air  are  presented  in  more 
general  terms  for  use  in  flow  calculations  (See 
References  25,  28-38  in  Chapter  1).  It  is  recalled 
that  the  thermodynamic  properties  such  as  enthalpy, 
specific  heat  and  composition,  etc.,  are  characteristics 
of  the  medium  which  are  not  dependent  upon  flow 
or  time.  On  the  other  hand,  the  transport  properties 


such  as  viscosity',  thermal  conductivity  and  diffusion 
coefficients  are  each  essentially  proportionality  fac¬ 
tors  between  gradients  of  quantities  such  as  velocity, 
temperature  and  component  concentrations  and 
their  respective  fluxes  such  as  shear  stress,  heat 
flux  and  component-diffusion  flux. 

The  thermodynamic  properties  of  equilibrium 
air  presented  here  in  tabular  and  graphical  form. 
Tables  4-1  through  4-4  and  Figures  4-1  through  4-10, 
(from  Reference  31  in  Chapter  l)may  also  be  rep¬ 
resented  in  the  form  cf  Mollier  diagrams.  Such 
Mollier  diagram  are  available  elsewhere  (References 
25,  28,  29  in  Chapter  1). 


TABLE  4-1.  TRANSPORT  PROPERTIES  OF  AIR  FOR  VARIOUS  TEMPERATURES 

AND  PRESSURES 


(A)  Coefficient  d  Viscosity 


1°  K 

Ratio  vh0t 

Reference  coefficient, 

v0i  in  10-« 

Pressure,  atmospheres 

100 

10 

1.0 

0.01 

0.001 

0.0001 

lb  sec 

gm 

ft j 

cm  sec 

500 

1.000 

L000 

1.000 

1.000 

1.000 

LOOO 

1.000 

0.558 

267 

KMX) 

1.000 

1.000 

KMX) 

1.  (MX) 

1.000 

LOOO 

LOOO 

.868 

416 

1500 

1.000 

1000 

1.000 

LOOO 

1.000 

LOOO 

LOOO 

1.100 

527 

2000 

1000 

1000 

1.000 

l.(XX) 

LOOO 

LOOO 

LOOO 

1.293 

619 

2500 

1.000 

1.000 

1.000 

1  (XX) 

LOOO 

LOOO 

LOOO 

1.461 

700 

3000 

1000 

looo 

1.000 

LOOO 

LOOO 

1.000 

LOOO 

1.612 

772 

3500 

1.000 

1.001 

1.003 

1.006 

1.010 

1.010 

1.011 

1.751 

838 

4000 

1.003 

1.008 

1.016 

1.020 

1.022 

1.024 

1.032 

1.879 

899 

4500 

1.010 

1.022 

1.029 

1.033 

1.038 

1.055 

1.096 

1.999 

957 

5000 

1.022 

1.036 

1.043 

1.051 

1.074 

1.128 

1.181 

2.11 

ion 

5500 

1.036 

1.052 

1.060 

1.086 

1.146 

1.209 

1.227 

2.22 

1062 

6000 

1.050 

1.067 

1.090 

1.148 

1.228 

1.257 

1.256 

2.32 

1112 

6500 

1.072 

1.090 

1.139 

1.229 

1.276 

1. 286 

1.271 

2.42 

1159 

7000 

1.089 

1.124 

1.208 

1.294 

1.317 

1.303 

1.264 

2.52 

1204 

7500 

1.112 

1.175 

1.283 

1.332 

1.337 

1.307 

1.210 

2.61 

1247 

3000 

1.143 

1.238 

1.342 

1.371 

1.347 

1.280 

1.072 

2.69 

1289 

8500 

1.185 

1.307 

1.386 

1.386 

1.343 

1.207 

.826 

2.78 

1330 

9000 

1.238 

1.368 

1.425 

1.396 

1.314 

1.068 

.517 

2.86 

1370 

9500 

1.298 

1.418 

1.438 

1.393 

1.251 

.853 

.261 

2.94 

1408 

10000 

1.361 

1  468 

1.445 

1.375 

1.143 

.595 

.118 

3.02 

1446 

10500 

1.418 

1.496 

1.448 

1.335 

.  983 

.361 

.055 

3.10 

1482 

11000 

1.467 

1.501 

1.442 

1.267 

.782 

.200 

.029 

3.17 

1518 

11500 

1.509 

1.511 

1.424 

1.168 

.571 

.  108 

.018 

3.24 

1552 

12000 

1.549 

1.520 

1.394 

1.040 

.387 

.063 

.012 

3.31 

1586 

12500 

1.577 

1.516 

1.342 

.881 

.249 

.036 

.009 

3.38 

1620 

13000 

1. 581 

1.508 

1.274 

.  711 

.  158 

.024 

.008 

3.45 

1652 

13500 

1.594 

1.492 

1.187 

.  547 

.  100 

.018 

.007 

3.52 

1684 

14000 

1.599 

1.468 

1.082 

.408 

.  067 

,015 

.007 

3.58 

1716 

14500 

1.601 

1.415 

.940 

.  268 

.042 

.013 

.008 

3.65 

1747 

15000 

1.604 

1.387 

.828 

.  212 

.016 

.012 

.008 

3.71 

1777 
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TABLE  4-1.  TRANSPORT  PROPERTIES  CF  AIR  FOR  VARIOUS  TEMPERATURES 
AND  PRESSURES  (continued) 

(B)  Coefficient  tf  Thermal  Conductivity 


T,°K 

Ratio  kjkt 

Reference  coefficient, 
k0,  in  10-4 

Pressure  atmospheres 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

Btu 

watt 
cm  °K 

ft  sec  °E 

500 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

5.84 

364 

1000 

1. 100 

1. 100 

1.100 

1  100 

1.100 

1.  100 

1. 100 

9.  10 

567 

1500 

1.150 

1.150 

1.150 

1  150 

1.150 

1.150 

1.150 

11.  53 

719 

2000 

1.177 

1.177 

1.177 

1.251 

1.460 

2.09 

3.99 

13.55 

844 

2500 

1.256 

1.317 

1.619 

2.50 

4.  63 

7.67 

5.50 

15.31 

954 

3000 

1.421 

1.928 

3.20 

5.48 

5.02 

2.  19 

1.465 

16.90 

1053 

3500 

1.941 

3  15 

4.72 

3.96 

1.719 

Z  11 

3.71 

18.35 

1143 

4000 

2.69 

3.94 

2.  99 

1.600 

2.91 

6.04 

15.03 

19.69 

1227 

4500 

3.22 

3.06 

1.714 

3.32 

7.34 

17. 65 

30.5 

21.0 

1305 

5000 

3.07 

1.997 

3.29 

7.  18 

16.63 

25.8 

11.84 

22.  1 

1379 

5500 

2.46 

2.91 

5.99 

13.71 

22.2 

11.40 

3.54 

23.  3 

1449 

6000 

1.930 

4.  53 

10.19 

18.74 

13.09 

3.96 

6.  14 

24.3 

1516 

6500 

3.35 

6.  98 

14.50 

15.39 

5.  49 

5.92 

12.99 

25.  4 

1580 

7000 

4.69 

9.97 

15.69 

8.32 

3.28 

10.95 

26.9 

26.4 

1642 

7500 

6. 31 

12.48 

12.24 

5.92 

8.62 

19.97 

51.9 

27.3 

1701 

8000 

8.  21 

13.19 

7.  80 

3.42 

13.99 

34.7 

84.8 

28.2 

1759 

8500 

9.86 

11.55 

5.  10 

9.  72 

22.  1 

56.3 

107.6 

29.1 

1814 

9000 

10.90 

8.79 

3.26 

14.03 

34.1 

78.8 

93.8 

30.0 

1868 

9500 

10.88 

6.38 

6.  66 

20,3 

49.  1 

92.6 

54.9 

30.8 

1921 

10000 

9.87 

4.07 

10.95 

28.2 

64.8 

84.8 

25.6 

31.6 

1972 

10500 

8.33 

6.50 

14.60 

37.6 

76.6 

59.9 

11.45 

32.4 

2020 

11000 

6.84 

8.47 

19.17 

48.  3 

78.6 

35.2 

5.81 

33.2 

2070 

11,500 

5.  59 

10.48 

24.  5 

58.2 

68.6 

19.04 

3.40 

34.0 

2120 

12000 

4.  79 

12.96 

30.9 

65.4 

51.6 

10.75 

2.29 

34.7 

2160 

12500 

3.34 

15.55 

37.0 

67.3 

34.8 

6.09 

1.807 

35.5 

2210 

13000 

5.78 

18.77 

43.3 

63.4 

22.5 

4.08 

1.622 

36.2 

2250 

13500 

9.95 

22.3 

48.7 

54.  5 

14.28 

3.18 

1.461 

36.9 

2300 

14000 

11.44 

25.9 

52.3 

43.  5 

9.58 

2.76 

1.586 

37.6 

2340 

14500 

13.73 

31.4 

54.2 

30.4 

6.25 

2.43 

1.670 

38.2 

2380 

15000 

14.95 

33.8 

52.4 

24.3 

3.04 

Z33 

1.754 

38.9 

2420 

TABLE  4-1.  TRANSPORT  PROPERTIES  CF  AIR  FOR  VARIOUS 
TEMPERATURES  AND  PRESSURES  (continued) 

(C)Prandtl  Number 


T,  °K 

Pressure!  atmospheres 

100 

10 

L0 

0.1 

0.01 

0.001 

0.0001 

500 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

1000 

.  756 

.756 

.  756 

.756 

.756 

.  756 

.756 

1500 

.767 

.767 

.767 

.767 

.767 

.767 

.  767 

2000 

.773 

.773 

.773 

.  706 

.724 

.  668 

.  614 

2500 

.  762 

.751 

.  696 

.645 

.611 

.654 

.771 

3000 

.  740 

.680 

.627 

.  636 

.740 

.745 

.714 

3500 

.678 

.  631 

.660 

.744 

.  737 

.658 

.600 

4000 

640 

.  662 

.  762 

.  759 

.019 

.580 

,587 

4500 

.654 

.743 

.752 

.  010 

.578 

.  611 

.764 

5000 

.  702 

.707 

.  611 

.581 

.624 

.799 

.  993 

5500 

.748 

.  620 

.  583 

.  617 

.785 

.  989 

.871 

6000 

.  763 

.  592 

.  602 

.  736 

.969 

.  891 

.  455 

6500 

.  610 

.592 

.  67  3 

.  906 

.  955 

.464 

.392 

7000 

.593 

.  620 

.  796 

.  986 

.830 

.404 

.361 

7500 

.595 

.  688 

.927 

.  969 

.424 

.371 

.  342 

8000 

.  620 

.  788 

.983 

.  648 

,  387 

.351 

.322 

8500 

.666 

.891 

.943 

,  411 

.363 

.  335 

.  279 

9000 

.730 

.961 

.807 

.382 

.348 

.316 

.  200 

9500 

.806 

.966 

.497 

.  364 

.336 

.279 

.  114 

10000 

.  886 

.872 

.429 

.  348 

.319 

.  216 

.0576 

10500 

.937 

.532 

.404 

.339 

.  295 

.  145 

.0314 

11000 

.  955 

.463 

:382 

.327 

.254 

.0877 

.0213 

11500 

.947 

.434 

.369 

.312 

.201 

.  0524 

.0167 

12000 

.90S 

.412 

.  355 

.  292 

.  146 

.0346 

.0143 

12500 

,  728 

.396 

.343 

.263 

.101 

.  0238 

.0129 

13000 

.  525 

,383 

.  333 

.227 

.  0688 

.  0190 

.0121 

13500 

,  438 

.369 

.319 

.  185 

.0470 

.0162 

.0110 

14000 

.421 

.360 

.302 

.  144 

.0345 

.0149 

.0108 

14500 

.  401 

,349 

.277 

.  0986 

.  0245 

.0130 

.0109 

15000 

.394 

.341 

.253 

.0819 

.  0129 

.0120 

.0110 

TABLE  4-2.  THERMAL  PROPERTIES  CF  AIR  FOR  VARIOUS  TEMPERATURES  AND  PRESSURES 


(A)  Specific  Heat 


(B)  Thermal  Conductivity  Ratio ,  fe„/fe0 


T, 


(§+') 


“K 


Pressure,  atmospheres 


Pressure,  atmospheres 


100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

3.  59 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1000 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

1.100 

1.100 

1.  100 

1.100 

1.100 

1.100 

1.100 

1500 

4.20 

4.20 

4.  20 

4.20 

4.20 

4.20 

4.  20 

1.150 

1.150 

1.150 

1.150 

1.  150 

1.150 

1.150 

2000 

4.33 

4.33 

4.33 

4.33 

4.33 

4.33 

4.34 

1,177 

1.177 

1.177 

1.177 

1.  177 

1.178 

1.181 

2500 

4.41 

4.41 

4.41 

4.41 

4.42 

4.45 

4.47 

1.195 

1.195 

1.196 

1.198 

1. 206 

1.225 

1.252 

3000 

4.46 

4.47 

4.47 

4.48 

4.51 

4.52 

4.  53 

1.  205 

1.208 

1.215 

1. 235 

1.268 

1.285 

1.289 

3500 

4.  50 

4.  52 

4.52 

4.  54 

4.  56 

4.55 

4.  56 

1.218 

1.230 

1.259 

1.289 

1.310 

1.314 

1.320 

4000 

4.  53 

4.55 

4.  57 

4.  ffc 

4.  58 

4.60 

4.65 

1.238 

1.270 

1.312 

1.329 

1.340 

1.356 

1.406 

4500 

4.57 

4.  59 

4.60 

4.61 

4.  64 

4.  74 

4.97 

1.269 

1.319 

1.346 

1.366 

1.396 

1.482 

1.674 

5000 

4.60 

4.61 

4.64 

4.  68 

4.81 

5.12 

5.46 

1.311 

1.357 

1.389 

1.428 

1.541 

Mil 

1.983 

5.500 

4.62 

4.66 

4.  71 

4.85 

5.  20 

5.  62 

5.76 

1.351 

1.406 

1.448 

1.569 

1.829 

2.08 

2.  14 

6000 

4.  64 

4.  71 

4.85 

5.  20 

5.73 

5.  94 

6.00 

1.383 

1.457 

1.566 

1.824 

2.  14 

2.  25 

2.50 

6.500 

4.  72 

4.82 

5.  12 

5.71 

6.  1 1 

6.  21 

6.  25 

1.459 

1.543 

1.767 

2.  14 

2.33 

2.  69 

2.  67 

7000 

4.78 

5.00 

5.53 

6.  16 

6.35 

6.44 

6.  56 

1.519 

1.685 

2.05 

2.38 

2.47 

2.86 

2.  79 

7500 

4.  89 

5.28 

6.02 

6.40 

6.58 

6.68 

6.97 

1.606 

1.889 

2.33 

2.  52 

3.14 

3.04 

2.90 

8000 

5.04 

5.68 

6.42 

6.  66 

6.  78 

6.  97 

7.  57 

1.733 

2.  14 

2.55 

2.  66 

3.33 

3.11 

2.75 

8500 

5.27 

6.06 

6.  66 

6.91 

6.97 

7.  34 

8.41 

1.898 

2.  40 

2.70 

3.64 

3.47 

3.29 

2.29 

9000 

5.  55 

6.42 

6.87 

6.98 

7.20 

7.  85 

9.35 

2.  10 

2.61 

2.83 

3.85 

3.56 

2.95 

1.771 

9500 

5.85 

6.  66 

6.98 

7.11 

7.  49 

8.54 

10.11 

2.32 

2.  78 

2.  96 

4.  03 

3.62 

2.  66 

1.307 

10000 

6.  15 

6.89 

7.07 

7.26 

7.86 

9.29 

10.52 

2.  53 

2.93 

3.03 

4.21 

3.54 

2.  23 

1.105 

10500 

6.41 

7.02 

7.  12 

7.  45 

8.34 

9.95 

10.73 

2.  71 

3.  15 

3.04 

4.34 

3.38 

1 .883 

1.061 

1 1000 

6.  59 

7.05 

7.20 

7.  67 

8.91 

10.41 

10.84 

2.  86 

3.  18 

3.06 

4.41 

3.  10 

1.643 

1.092 

11.500 

6.71 

7.06 

7.33 

7.  97 

9.  49 

10.67 

10.91 

2.  97 

3.24 

3.02 

4.  45 

2.  84 

1.582 

1.152 

12000 

6.76 

7.07 

7.39 

8.33 

9.  99 

10.82 

10.94 

3.06 

3.29 

2.99 

4.  39 

2.  59 

1.573 

1.223 

12500 

6.88 

7.  10 

7.  54 

8.75 

10.37 

10.90 

10.98 

3.  17 

3.26 

2.88 

4.28 

2.44 

1.643 

1.301 

13000 

6.86 

7.  11 

7.71 

9.20 

10.62 

10.96 

10.99 

3.37 

3.27 

2.  76 

4.  13 

2.  42 

1.702 

1.381 

13500 

6.91 

7.  17 

7.95 

9.03 

10.77 

10.99 

ii.  ai 

3.51 

3.27 

2.  58 

3.  97 

2.  43 

1.814 

1.461 

14000 

6.  90 

7.23 

8.  24 

10.01 

10.87 

11.00 

1 1.02 

3.  55 

3.20 

2.41 

3.90 

2.  49 

1.892 

1.586 

14500 

6.90 

7.37 

8.66 

10.37 

10.94 

11.01 

11.02 

3.  56 

3.12 

2  11 

3.  79 

2.56 

2.01 

1.670 

15000 

6.  88 

7.44 

8.  88 

10.52 

11.02 

11.02 

11.03 

3.59 

3.09 

1.930 

3.84 

2.  64 

2.  13 

1.  754 

4-3 


TABLE  4-3.  THERMODYNAMIC  PROPERTIES  CF  AIR  FOR  VARIOUS  TEMPERATURES  AND  PRESSURES 


(A)  Compressibility ,  Z 


(B)  Dimensionless  Enthalpy,  ZH/RT 


T.°  K 

Pressure,  atmospheres 

Pressure,  atmospheres 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

1  000 

1000 

1000 

1.000 

1000 

1000 

1.000 

3.52 

3.52 

3.52 

3.  52 

3.52 

3.52 

3.52 

1000 

1000 

1000 

1000 

1000 

1000 

1000 

1  000 

3.65 

3.65 

3.65 

3.65 

3.65 

3.65 

3.65 

1500 

1.000 

1000 

1.000 

1000 

1000 

1.000 

1.000 

3.80 

3.80 

3.80 

3.80 

3.80 

3.80 

3.80 

2000 

1000 

1000 

1.000 

1.001 

1.002 

1.005 

1.016 

3.92 

3.92 

3.92 

3.93 

3.97 

4.07 

4.41 

2500 

1000 

1.001 

1.004 

1.011 

1.033 

1.088 

1.163 

4.01 

4  03 

4  09 

4.27 

4.81 

6.16 

8.02 

3000 

1.003 

1.009 

1.026 

1.072 

1.149 

1.192 

1.200 

4.13 

4.25 

4.61 

5.55 

7.13 

8  02 

8.19 

3500 

1.012 

1.035 

1.092 

1.167 

1.197 

1.203 

1.211 

4.34 

4.75 

5.75 

7.08 

7.62 

7.77 

8.03 

4000 

1.033 

1.089 

1.165 

1.198 

1.208 

1.228 

1.287 

4.70 

5.56 

6.74 

7.28 

7.53 

8.09 

9.82 

4500 

1.071 

1.149 

1.196 

1.213 

1.245 

1.337 

1.577 

5.20 

6.29 

6.98 

7.33 

8.14 

10.55 

16.80 

5000 

1.118 

1.186 

1.214 

1.252 

T.  359 

1.622 

1.910 

5.73 

6.62 

7.10 

7.96 

10.48 

16.68 

23.46 

5500 

1.159 

1.208 

1.348 

1.348 

1.599 

1.898 

1.990 

6.13 

6.80 

7.58 

9.73 

15.14 

21.58 

23.58 

6000 

1.189 

1.235 

1.316 

1.529 

1.849 

1.9S3 

2.008 

6.38 

7.11 

8.  70 

12.93 

19.30 

21.97 

22.  54 

6500 

1.214 

1.279 

1.437 

1.752 

1.961 

2.006 

2.032 

6.62 

7.72 

10.64 

16.46 

20.35 

21.24 

21.93 

7000 

1.243 

1.351 

1.607 

1.904 

1.997 

2.027 

2.088 

6.95 

8.76 

13.20 

18.34 

20.01 

20.69 

22.29 

7500 

1.284 

1.457 

1.778 

1.971 

2.017 

2.067 

2.210 

7.44 

10.24 

15.48 

18.66 

19.54 

20.72 

24.26 

8000 

1.341 

1.590 

1.896 

2.001 

2. 044 

2.144 

2.446 

8.16 

11.99 

16.73 

18.43 

19.34 

21.65 

28.65 

8500 

1.418 

1.727 

1.959 

2.023 

2.090 

2. 284 

2.826 

9.10 

13.63 

17.09 

18.17 

19.60 

23.85 

35.75 

9000 

1.512 

1.838 

1.993 

2.050 

2.166 

2.510 

3.282 

10.20 

14.79 

17.04 

18.09 

20.49 

27.66 

43.74 

9500 

1.616 

1.914 

2.018 

2.090 

2.286 

2.832 

3.645 

11.36 

15.40 

16.91 

18.29 

22.17 

33.00 

49.15 

10000 

1.718 

1.962 

2.042 

2.149 

2.462 

3.202 

3.843 

12.42 

15.61 

16.84 

18.85 

24.78 

38.79 

50.96 

10500 

1.807 

1.993 

2.071 

2.234 

2.700 

3.526 

3.932 

13.23 

15.64 

16.90 

19.84 

28.28 

43.28 

50.64 

11000 

1.876 

2.018 

2.111 

2.351 

2.983 

3.745 

3.969 

13.77 

15.60 

17.13 

21.31 

32.31 

45.57 

49.48 

11500 

1.927 

2.042 

2.163 

2  505 

3.272 

3.867 

3.985 

14.08 

15.58 

17. 57 

23.28 

36.13 

46.09 

48.07 

12000 

1.965 

2.067 

2.232 

2.694 

3.520 

3.931 

3.993 

14.22 

15.62 

18.24 

25.69 

39.01 

45.64 

46.64 

12500 

1.993 

2.098 

2.318 

2.910 

3.700 

3.963 

3.996 

14.28 

15.74 

19.16 

28.36 

40.66 

44.74 

45.26 

13000 

2.017 

2.135 

2.426 

3.135 

3.818 

3.979 

3.998 

14.30 

15.96 

20.32 

30.99 

41.26 

43.69 

43.97 

13500 

2.039 

2.180 

2.553 

3.347 

3.889 

3.988 

3. 999 

14.31 

16.28 

21.72 

33.27 

41.17 

42.61 

42.76 

14000 

2.062 

2.233 

2.700 

3.527 

3.932 

3.993 

3.999 

14.34 

16.71 

23.29 

34.97 

40.69 

41.55 

41.64 

14500 

2.086 

2.297 

2.861 

3.667 

3.957 

3.996 

4000 

14.40 

17.26 

24.98 

36.02 

40.01 

40.53 

40.58 

15000 

2.113 

2.372 

3.028 

3.769 

3.973 

3.997 

4000 

14.49 

17.  92 

26.66 

36.53 

39.24 

39.57 

39.60 

(C)  Dimensionless  Entropy,  ZS/R 


(D)  Dimensionless  Specific  Heat  at 
Constant  Pressure,  ZCV/R 


500 

21.1 

23.4 

25.7 

28.  0 

30.3 

32.6 

34.9 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

1000 

23.7 

26.0 

28.3 

30.6 

32.9 

35.2 

37.5 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

1500 

25.3 

27.6 

29.9 

32.2 

34.5 

36.8 

39.2 

4.20 

4.20 

4.20 

4.20 

4.21 

4.23 

4.29 

2000 

26.6 

28.9 

31.2 

33.5 

35.8 

38.2 

40.  9 

4.34 

4  35' 

4.41 

4  57 

5.10 

6.73 

11.57 

2500 

27.  5 

29.9 

a2.2 

34.7 

37.6 

41.4 

45.9 

4.51 

4.74 

5.46 

7.63 

13.58 

24.00 

20.03 

3000 

28.4 

30.8 

33.5 

36.9 

41.0 

44.6 

47.5 

5.03 

6.23 

9.62 

16.64 

17.66 

7  88 

5.41 

3500 

29.3 

32.0 

35.5 

39.4 

42.7 

45.6 

48.6 

6.27 

9.50 

14.90 

12.31 

6.51 

6.64 

10.69 

4000 

30.2 

33.5 

37.3 

40.6 

43.6 

46.  9 

51.6 

8.24 

12.47 

11.05 

6.91 

8.49 

16.46 

40.78 

4500 

31.3 

34.9 

38  3 

41.5 

45.1 

50.5 

60.0 

10.12 

11.06 

7.63 

9.61 

19.68 

49.06 

101.49 

5000 

32.4 

36.0 

39.2 

42.9 

48  4 

58.0 

68.9 

10.53 

8.57 

9.58 

18.98 

46.22 

87.18 

47.49 

5500 

33.4 

36.8 

40.4 

45.5 

54.3 

61.8 

71.3 

9.56 

9.06 

15.85 

37.24 

72.70 

44.58 

12.82 

6000 

34.2 

37.7 

42.2 

49.7 

59.9 

€1.1 

72.2 

9.05 

12.38 

26.95 

57.60 

49.37 

15.16 

11.43 

6500 

34.9 

38.9 

44.9 

54.4 

62.6 

68.1 

73.4 

10.05 

18.31 

41.13 

54.37 

20.60 

11.76 

19.13 

7000 

35.8 

40.5 

48.3 

57.6 

63.8 

69.1 

75.4 

12.59 

26.50 

49.61 

31.18 

12.89 

16.38 

36.82 

7500 

36.8 

42.7 

51.6 

59.2 

64.6 

70.5 

79.0 

16  48 

35.19 

42.69 

17.32 

13.99 

27.02 

70.02 

8000 

38.0 

45.1 

53.9 

60.1 

65.7 

72.8 

85.0 

21.42 

40.33 

28.27 

13.68 

19.40 

45.62 

121.60 

8500 

39.4 

47.6 

55.3 

61.0 

€7.1 

76.4 

94.1 

26.70 

38.14 

18.58 

14.91 

28.85 

74.32 

173.36 

9000 

41.1 

49.  5 

56.2 

61.9 

69.2 

81.6 

104.3 

31.05 

30.38 

14.83 

18.90 

43.11 

111.46 

173.32 

9500 

42.  8 

51.0 

57.0 

63.1 

72.0 

88.6 

112.2 

32  90 

22.50 

14.65 

25.26 

62.69 

143.94 

114.86 

10000 

44.5 

52.0 

57.8 

64.6 

75.8 

96.2 

116.6 

31.43 

17.46 

16.53 

34.11 

86.43 

147.10 

60.  29 

10500 

45.9 

52.8 

58.7 

66.6 

80.6 

102  7 

118.8 

27.46 

15.19 

19.86 

45.53 

109.34 

114.76 

31.99 

11000 

47.1 

53.4 

59.7 

69.0 

86.0 

107.1 

119.9 

22.82 

14.83 

24.46 

59.25 

121.78 

73.60 

19.92 

11500 

48.0 

54. 1 

60.9 

71.9 

91.4 

109.6 

120.7 

18.96 

15.71 

30.27 

74.10 

115.48 

44.18 

14.93 

12000 

48.8 

54.8 

62.4 

75.4 

95.8 

111.1 

121.3 

16.39 

17.45 

37.23 

87.57 

93.38 

27.91 

12.82 

12500 

49.4 

55.6 

64.0 

79. 1 

99.1 

112.1 

121.8 

15.01 

19.88 

45.17 

96.01 

67.36 

19.67 

11.89 

13000 

50.0 

56.4 

66.0 

82.9 

101.3 

112.8 

122.3 

14.57 

22.89 

53.70 

96.15 

46.35 

15.56 

11.47 

13500 

50.5 

57.3 

68.2 

86.4 

102.8 

113.3 

122.1 

14.80 

26.42 

62.12 

87.59 

32.30 

13.49 

11.26 

14000 

51.1 

58.4 

70.5 

89.4 

103.8 

113.8 

123.1 

15.52 

30.41 

69.41 

73.42 

23.72 

12.41 

11.16 

14500 

51.  7 

59.5 

73.1 

91.7 

104.5 

114.2 

123.5 

16.62 

34.77 

74.36 

58.  06 

18.66 

11.83 

11.10 

15000 

52.2 

60.8 

75.6 

93.4 

105.1 

114.6 

123.9 

18.02 

39.37 

75.95 

44.65 

15.70 

11.51 

11.08 
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TABLE  4-3.  THERMODYNAMIC  PROPERTIES  OF  AIR  FOR  VARIOUS  TEMPERATURES  AND  PRESSURES 

(continued) 


(E)  Ratio  <£  Specific  Heats,  y 


( F)  Dimensionless  Speed  of  Sound 
Parameter,  a2  p/ p 


T, 

°K 

Pressure,  atmospheres 

Pressure,  atmospheres 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

1.387 

1.387 

L  387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1.387 

1000 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1.337 

1500 

1.312 

1.312 

1.312 

1.312 

1.312 

1.310 

1.306 

1.312 

1.312 

1.312 

1.312 

1.312 

1.310 

1.306 

2000 

1.300 

1.299 

1.296 

1.286 

1.260 

1.209 

1.153 

1.300 

1.299 

1.296 

1.285 

1.259 

1.206 

1.144 

2500 

1.288 

1.277 

1.249 

1.202 

1.161 

1.152 

1.157 

1.288 

1.276 

1.247 

1.196 

1.144 

1.119 

1.132 

3000 

1.266 

1.235 

1.195 

1.178 

1.181 

1.239 

1.304 

1.265 

1.229 

1.181 

1.147 

1.150 

1.232 

1.302 

3500 

1.241 

1.211 

1.202 

1.212 

1.270 

1.252 

1.176 

1.234 

1.192 

1.166 

1.187 

1.265 

1.250 

1.171 

4000 

1.230 

1.223 

1.230 

1.260 

1.213 

1,150 

1.133 

1.212 

1.187 

1.204 

1.254 

1.208 

1.137 

1.097 

4500 

1.240 

1.243 

1.251 

1.204 

1.154 

1.155 

L  190 

1.208 

1.210 

1.241 

1.196 

1.133 

1.101 

1.092 

5000 

1.256 

1.252 

1.212 

1.166 

1.172 

1.203 

1.168 

1.217 

1.233 

1.202 

1.143 

1.111 

1.103 

1.124 

5500 

1.262 

1.231 

1.183 

1.182 

1.214 

1.  183 

1.257 

1.229 

1.217 

1.161 

1.124 

1.113 

1.133 

1.249 

6000 

1.253 

1.206 

1.190 

1.221 

1.202 

1.237 

1.266 

1.230 

1.186 

1.141 

1.124 

1.135 

1.225 

1.263 

6500 

1.235 

1.201 

1.220 

1.228 

1.217 

1.265 

1. 188 

1.214 

1.165 

1.137 

1.136 

1.193 

1.260 

1.183 

7000 

1.223 

1.217 

1.246 

1.216 

1.2.58 

1.210 

1.155 

1.195 

1.155 

1.142 

1.167 

1.249 

1.205 

1.140 

7500 

1.223 

1.243 

L  244 

1.237 

1.237 

1.173 

1.164 

1.182 

1. 154 

1.156 

1.216 

1.231 

1.162 

1.128 

8000 

1.235 

1.264 

1.235 

1.252 

1.201 

1.168 

1.201 

1.175 

1.159 

1.181 

1.242 

1.193 

1.143 

1.130 

8500 

1.255 

1.267 

1.243 

1.235 

1.183 

1.188 

1.242 

1.174 

1.169 

1.214 

1.228 

1.168 

1.140 

1.136 

9000 

1.275 

1.260 

1.252 

1.213 

1.185 

1.224 

1.244 

1.176 

1.185 

1.237 

1.203 

1.157 

1. 144 

1.145 

9500 

1.288 

1.255 

1.248 

1.201 

1.203 

1.256 

1.216 

1.181 

1.206 

1.237 

1.185 

1.155 

1.151 

1.157 

10000 

1.291 

1.259 

1.236 

1.201 

1.232 

1.263 

1.211 

1.190 

1.228 

1.225 

1.176 

1.159 

1.159 

1.185 

10500 

1.287 

1.262 

1.226 

1.213 

1.263 

1.244 

1.2.56 

1.201 

1.242 

1.212 

1.173 

1.165 

1.169 

1.244 

11000 

1.282 

1.261 

1.222 

1.233 

1.281 

1,230 

1.339 

1.216 

1.246 

1.202 

1.175 

1.173 

1.187 

1.334 

11500 

1.280 

1.256 

1.226 

1.258 

1.280 

1.243 

1.427 

1.232 

1.242 

1.198 

1.180 

1.181 

1.221 

1.424 

12000 

1.282 

1.252 

1.236 

1.283 

1.267 

1.288 

1.491 

1.247 

1.236 

1.197 

1.187 

1.192 

1.276 

1.489 

12500 

1.284 

L  250 

1.251 

1.301 

1.257 

1.352 

L  528 

1.257 

I.  231 

1.198 

1.194 

1.207 

1.346 

1.528 

13000 

1.285 

1.251 

1.271 

1.307 

1.263 

1.419 

1.  ,548 

1.263 

1.227 

1.202 

1.202 

1.232 

1.415 

1.548 

13500 

1.284 

1.257 

1.291 

1.303 

1.288 

1.472 

L558 

1.265 

1.226 

1.208 

1.210 

1.269 

1.470 

1.  558 

14000 

1.284 

1.266 

1.311 

1.295 

1  329 

L  509 

1.563 

1.265 

1.226 

1.214 

L  221 

1.317 

L  508 

1.  .563 

14500 

1.284 

1.278 

1.326 

1.290 

1.377 

1.532 

L  565 

L  264 

1.228 

1.221 

1.235 

1.370 

1.532 

1.565 

15000 

1.286 

1.293 

1.336 

1.263 

1.425 

1.547 

1.  567 

1.263 

1.232 

1.228 

1.254 

1.420 

1.546 

1.567 
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TABLE  4-4.  CONVERSION  TABLE  FOR  OBTAINING  THERMODYNAMIC  PROPERTIES 
EKCM  THE  DIMENSIONLESS  QUANTITIES 


Multiply 

By 

To  get 

In 

ZE  ZH 
BTm  RT 

0-  0686  T  (°K) 

0.287  T  (°K) 

0.0686  T  (°R) 

1716  T  (°R) 

Energy  or  enthalpy 

cal/gm 

joule/gm 

Btu/lb  mass 
ft-lb/slug 

Z8  ZC. 

R  '  R  ' 

7.C.  ZCZ 

nr 01 -ir 

0.0686 

0.287 

0.0686 

1716 

Entropy  or  specific  heat 

cal/gm  °K 
joule/gm  °K 
Btu/lb  mass  °R 
ft-lb/slug  °R 

i/i« 

1.462X10-'  T'/»  (l+~)_1 
(Tin°K) 

2. 28X10-*  7’1/a  (j+^r)’1 

Coefficient  of  viscosity 

gm/cm-sec 

II)  sec/ft' 

(Tln/R) 

4.76X10“*  7’1/a  (l-t-^)~‘ 

cal 

cm-sec  °K 

(T  In  °K) 

1.994X10-*  TO*  (l+'p)1 

watt 
tan  °K 

(T  In  °K) 

t  km 
-r-OT  T" 

k$ 

2.39X10-'  7'1/a 

Coefficient  of  thermal  conductivity 

Btu 

ft-sec  °R 

(rin°R) 

1.856X10-*  Ti/a  (i+^r)"1 

ft-lb 

ft-sec  °R 

(Tin  °R) 

4-6 


P» it  tude 


Figure  4-1.  Stagnation  Temperature  and  Pressure  Figure  4-2.  Domains  cf  Pressure  and  Temperature 
in  Air  as  a  Function  of  Altitude  and  Velocity  for  the  Major  Chemical  Reactions  in  Air 
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Figure  4-5.  Energy  of  Air  as  a  Function  of  Temperature 
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Figure  4-6.  Zero  Frequency  Speed  d  Sound  Parameter  for  Air  as  a  Function  of  Temperature 


Figure  4-7.  Compressibility  d  Air  as  a  Function  cf  Temperature 
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Figure  4-8.  Ratio  tf  the  Coefficient  cf  Thermal  Conductivity  cf  Air  to  the 
Reference  Coefficient ,  fc0,  as  a  Function  of  Temperature 
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Figure  4-10.  Prandtl  Number  for  Air  as  a  Function  of  Temperature 


CHAPTER  5 

PRESSURE  DISTRIBUTIONS 


5-1  INTRODUCTION 

The  pressure  distribution  over  the  surface  of 
a  vehicle  must  be  known  to  permit  an  evaluation 
of  the  skin  friction,  heat  transfer,  and  other  diffusive 
phenomena  which  occur  in  the  boundary  layer  over 
the  vehicle,  and  to  permit  an  evaluation  of  the  forces 
other  than  those  due  to  skin  friction.  A  detailed 
determination  of  the  pressure  distribution  requires 
complex  calculations  of  the  flow  field  about  the  body. 
Methods  by  which  such  complete  analyses  may  be 
carried  out  are  the  subject  of  numerous  standard 
works  on  aerodynamics  and  are  outside  the  scope 
of  this  work.  However,  approximate  methods  have 
been  developed  which  give  estimates  with  accuracy 
sufficient  for  present  purposes  and  indeed  for  pre¬ 
liminary  design.  The  main  part  of  this  chapter  will 
be  devoted  to  these  approximate  procedures  to 
provide  a  “feel”  for  the  possible  values  which  may 
be  encountered. 

5-2  “NEWTONIAN”  PRESSURE 
DISTRIBUTION 

Primary  interest  for  ballistic  missile  applications 
centers  on  the  hypersonic  speed  range,  which  for 
present  purposes  can  be  defined  as  being  in  excess 
cf  5000  ft/sec.  In  this  regime  surfaces  which  make 
significant  contributions  to  the  forces  acting  on  the 
body  are  usually  inclined  at  angles  8  with  the  flight 
velocity  vector  V  which  are  not  “very  small”.  This 
is  particularly  true  of  blunt  noses  which  are  roughly 
spherical  in  shape,  of  conical  surfaces  whose  angles 
are  not  small,  and  of  the  surfaces  of  inclined  bodies. 
Moreover,  on  bodies  of  the  type  usually  employed, 
surfaces  which  are  “shaded”  from  the  oncoming 
velocity  vector  (i.e.,  those  whose  surface -normals 
make  oblique  angles  with  V)  usually  make  minor 
contributions  to  the  forces  and  heat  inputs  to  the 
body. 

As  a  result,  the  following  empirical  estimate  can 
be  applied  as  a  first  approximation  to  surfaces  which 
are  “seen”  by  the  vector  V : 


Here  p,  denotes  the  undisturbed  pressure  at  a 
given  flight  altitude  and  p ,  is  the  maximum  pressure 
to  be  found  on  the  body  surface  arid  is  usually  equal 
to  the  stagnation  pressure  on  the  body  (i.e.,  the 
pressure  at  a  point  where  the  velocity  is  zero  on 
a  body  of  revolution).  A  reasonable  estimate  for 
p,  at  hypersonic  speeds  is  given  as  follows: 

p,  ~pF2»p/  (5-2) 

where  p  is  the  density  at  the  flight  altitude  and 
V  is  the  flight  velocity. 

The  pressure  relation  cited  above  is  termed  the 
“Newtonian”  distribution  because  of  associated  rea¬ 
soning  cf  Newton  who  deduced  a  form  of  this  kind 
by  assuming  that  the  component  cf  momentum 
normal  to  the  surface  in  a  molecular  gas  is  completely 
lost  on  impact.  Many  theoretical  discussions  cf  the 
rigorous  range  and  conditions  of  validity  of  this 
relation  can  be  found  in  the  literature,  e.g.,  Ref.  7 
cf  Chapter  1.  However,  it  is  presented  here  in  an 
empirical  sense  without  further  justification. 

It  is  seen  that  when  8  — >  0  the  Newtonian  rela¬ 
tion  predicts  that  p  — >  p„  a  result  which  is  not 
quite  correct  because  of  the  presence  cf  a  shock  wave 
around  the  body.  As  a  result,  another  method  is 
often  used  to  continue  the  estimate  cf  pressures  into 
regions  of  small  surface  inclination.  This  method  cf 
estimation  employs  the  “Prandtl-Mcyer  Expansion” 
relation  which  is  discussed  in  Paragraph  5-3. 

The  pressure  alone  is  not  sufficient  to  permit 
the  calculation  of  the  velocity  and  other  properties 
at  the  body  surface  which  are  required  for  estimating 
the  heat  transfer  and  other  boundary  layer  char¬ 
acteristics.  This  procedure  requires  a  knowledge  cf 
the  equations  cf  motion  governing  isentropic  flow 
and  may  be  found  in  most,  textbooks  of  aero¬ 
dynamics.  An  indication  of  the  general  nature  cf 
the  procedure  can  be  given  for  the  highly  simplified 
conditions  wherein  there  are  no  chemical  reactions 
occurring  in  the  gas,  and  when  the  body  is  blunted 
so  that  a  stagnation  point  exists.  In  this  case 

7T  =  /( 7,  Ml)  (5-3) 

r  • 

where 
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Ml  is  the  local  Mach  number  VL/aL, 

V L  is  the  local  velocity  on  the  surface, 
aL  ~  50  Tl  is  the  local  sound  speed  in  ft/sec, 

TL  is  the  local  temperature  in  °R,  and 
y  is  the  ratio  of  specific  heats  for  the  air  mixture 
(1.4  at  low  temperatures). 

Thus  a  knowledge  of  p/p,  and  y  determines 
Ml.  Next,  energy  conservation  considerations  give 
a  relation 

=  g(yMD  (5-4) 

where 

T,  is  termed  the  stagnation  temperature  and  has 
the  approximate  value  V2/2Cv, 

CP  is  the  specific  heat  at  constant  pressure,  and 
V  is  the  flight  velocity. 

In  this  manner,  the  local  properties  of  the  flow 
at  the  surface  of  the  body  can  be  determined.  For 
flows  with  chemical  reactions  the  estimation  of  local 
properties  may  be  carried  out  with  the  use  of  a 
Mollier  chart  when  the  flow  is  in  chemical  equilib¬ 
rium.  Otherwise,  the  equations  governing  chemical 
rate  processes  must  be  coupled  with  the  flow  equa¬ 
tions  to  permit  the  remaining  properties  to  be 


evaluated.  Pertinent  tabulations  of  flow  properties 
have  been  cited  in  Chapter  1  as  References  22-25. 

A  final  word  of  caution  is  in  order  concerning 
the  value  of  p,  which  is  required  in  conjunction  w  ith 
p  to  give  local  flow  properties.  For  a  blunt  body  in 
which  the  viscous  boundary  layer  is  thin  the  fore¬ 
going  procedure  can  be  applied  with  good  accuracy. 
However,  on  a  sharp-nosed  body,  where  the  stagna¬ 
tion  point  is  enveloped  within  the  boundary  layer 
near  the  sharp  apex,  the  value  cf  p,  may  be  some¬ 
what  higher  than  pV 2  Indeed  it  rigorously  would 
be  the  ‘"stagnation”  pressure  behind  the  inclined 
bow  shock  along  a  streamline  which  pierces  the 
boundary  layer  at  the  local  point  of  interest  on  the 
body.  For  example,  on  a  sharp  cone,  the  value  cf 
p,  to  be  used  in  Equation  5-3  is  that  corresponding 
to  the  oblique  shock  about  the  body. 

5-3  “PRANDTL-MEYER”  EXPANSION  AND 
FLOW  FELDS  OF  CONES 

Exact  solutions  of  the  equations  of  motion  can 
be  carried  out  for  special  ty  pes  of  flow.  These  solu¬ 
tions  are  useful  in  estimating  pressures  over  bodies. 
One  cf  these,  the  Prandtl-Meyer  expansion,  provides 
a  relation  between  the  local  surface  inclination  5 
and  the  local  Mach  number  M L  for  a  given  initial 


Ml 

Figure  5-1.  Prandtl-Meyer  Expansion  Relation 
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set  of  values  5,  and  ML{  and  for  a  given  y.  There¬ 
fore  if  one  ends  an  estimate  of  p  at  a  given  point 
on  the  surface  where  6  and  ML  are  known,  one  can 
continue  the  estimate  for  usual  values  cf  6  by  means 
of  the  Prandtl-Meyer  relation. 

A  plot  of  this  relation  for  y  =  1.4  is  given  in 
Figure  5-1.  Tabulations  for  other  values  cf  y  are 
available  in  the  literature.  Some  calculations  are 
also  available  which  take  into  account  equilibrium 
high  temperature  air  chemistry  . 


Tabulations  of  the  properties  of  flows  about 
sharp  cones  which  represent  another  useful  set  cf 
exact  aerodynamic  calculations  can  be  found  in 
NACA  TR  1135  (Reference  4)  and  are  described 
in  standard  texts  of  supersonic  aerodynamics.  The 
pressures  on  cones  are  sometimes  used  empirically 
to  give  the  pressure  on  a  body  of  revolution  whose 
inclination  y  is  assumed  to  be  part  of  a  fictitious 
cone  cf  half  angle  6. 

An  example  of  how  some  of  these  methods  are 


Figure  5-2.  The  Cone-Cylinder- Flare  Configuration,  Overall  Mode? 


SURFACE  DISTANCE,  S 

Figure  5-3.  Analytical  and  Experimental  Pressure  Distributions3 
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applied  to  a  typical  missile  shape  to  permit  a  de¬ 
termination  of  theoretical  heat  rates  is  given  in 
Reference  3.  For  the  cone-cylinder-flare  configura¬ 
tion  shown  in  Figure  5-2,  pressures  were  measured 
at  Mf  =  8.  The  experimentally  determined  pressures 
are  shown  in  Figure  5-3  together  with  approximate 
theoretical  estimate.  Methods  used  are  indicated  on 
the  figures. 
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CHAPTER  6 


THERMAL  EFFECTS 


6-1  INTRODUCTION 

In  high  speed  flight,  high  temperatures  are 
generated  around  a  missile  and  consequently  within 
its  structural  components.  Conditions  on  the  exposed 
surface  of  the  missile  form  the  connecting  link  be¬ 
tween  the  air  flow  and  the  interior  thermal  state 
of  the  structure.  Therefore,  interfacial  conditions 
such  as  surface  temperature  or  rate  of  heat  flow 
through  the  surface  must  be  expressed  in  terms  cf 
the  surface  contour  and  flight  parameters  (e.g., 
flight  velocity,  etc.).  On  the  other  hand,  the  surface 
conditions  must  be  related  to  properties  within  the 
structure  so  that  the  internal  temperature  distribu¬ 
tion  can  be  established.  The  joining  of  these  two 
classes  of  information  through  their  common  bound¬ 
ary  condition  provides  the  relation  between  the 
flight  situation  and  the  thermal  characteristics  of 
the  structure. 

The  first  step  in  estimating  structural  effects, 
such  as  thermal  stresses,  consists  of  determining  the 
internal  temperature  level  or  temperature  distribu¬ 
tion.  Theoretically,  this  problem  is  governed  by  the 
classical  heat  conduction  equation  which  expresses 
the  energy  balance  within  the  material.  The  follow¬ 
ing  parameters  influence  the  temperature  distribu¬ 
tion  :  (i)  structural  geometry;  (ii)  material  properties, 
i.e.,  both  terrnal  conductivity  k(T)  and  thermal 
diffusivity  a(T)  =  k/wc  (where  w  is  the  specific 
weight  and  c  the  specific  heat)  in  time-dependent 
cases,  and  only  k(T)  in  steady-state  cases;  (iii)  rate 
of  energy  transfer  through  the  boundaries  or  energy 
generated  within  the  material;  or  alternatively,  the 
surface  temperature,  which  implies  the  heat  transfer 
necessary  to  maintain  this  temperature. 

Using  the  known  temperature  distribution  or 
temperature  level  in  the  structure,  one  can  examine 
excessive  temperatures  or  other  destructive  condi¬ 
tions,  probable  cooling  requirements,  possibility'  of 
surface  melting  or  ablation  (mass  loss)  which  may 
be  useful  for  cooling,  as  well  as  evaluating  thermal 
stresses  and  required  thermal  materials.  Changes  in 
body-geometry  or  flight-programs  may  be  employed 
to  alleviate  or  circumvent  an  unfavorable  situation. 


6-2  CONDITIONS  FOR  APPLICATION  OF 
STEADY  FLOW  ANALYSIS 

It  is  worthwhile  to  recall  here  that  flow  and 
transfer  processes  are  usually  viewed  from  a  coor¬ 
dinate  system  fixed  to  the  body.  In  this  frame  of 
reference  a  body  moving  through  stationary  air  is 
aerodynamically  equivalent  to  the  same  body  held 
fixed  with  respect  to  the  ground  and  subjected  to 
an  airstream  ;  the  ground-reference  approximates  an 
inertial  coordinate  system.  The  airstream  is  initially 
uniform,  undisturbed,  parallel  to  the  flight  velocity 
vector,  and  equal  in  velocity  to  the  flight  velocity. 
The  state  parameters  of  this  relative  flow,  that  is, 
its  pressure,  temperature  and  density,  are  equal  to 
those  of  the  stationary  air  at  a  given  altitude.  In  this 
coordinate  system  the  flow  may  be  regarded  as 
steady  or  at  least  quasi-steady  for  all  but  the  most 
refined  heat  transfer  purposes.  By  quasi-steady  we 
mean  that  time-dependent  flows  may  be  treated  as 
a  continuous  series  of  instantaneously  steady  states, 
and  that  explicit  time-derivatives  of  the  flow  field 
are  negligible,  with  respect  to  spatial  derivatives. 

Some  effects  of  unsteadiness,  such  as  those  due 
to  the  acceleration  of  a  missile,  can  be  taken  into 
account  by  considering  the  consequent  inertia  forces 
to  act  as  additional  body  forces,  augmenting  the 
body  forces  such  as  those  cf  gravity.  Only  recently 
have  studies  of  lime-dependent  boundary  layer  and 
heat  transfer  phenomena  received  appreciable  atten¬ 
tion.  For  example,  special  types  of  oscillating  flows 
have  been  studied1-4  and  have  indicated  effects  cf 
unsteadiness  on  shear  and  heat  transfer  which  are 
important  only  in  a  secondary  way.  More  general 
treatments  of  the  effects  cf  time-dependent  flight 
on  boundary  layer  phenomena  including  heat  trans¬ 
fer  have  recently  appeared5-7  and  will  be  discussed 
further  in  a  later  chapter. 

A  useful  approximate  criterion  for  judging  when 
the  direct  effects  of  unsteadiness  may  be  neglected 
in  boundary  layers  is  that  the  following  inequality 
be  satisfied: 

#W«1’  (6-D 
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where  ut  is  a  characteristic  velocity,  (e.g.,  the  veloc¬ 
ity  external  to  the  boundary  layer  at  the  station  x), 
x  a  streamwise  distance  and  t  denotes  time.  These 
conditions  are  usually  satisfied  in  current  missile 
applications.  For  example,  for  a  body  100  ft  long, 
traveling  at  1000  ft/sec  with  an  acceleration  of  30  g 

-2  <■<  1 ,  for  7i—l  (6-2) 

u,  at  10 

Shock-tube  experimentation  has  indicated  that 
flows  over  forebodies  become  established  in  the  order 
of  tens  of  microseconds  in  high-speed  flow.  Wakes 
in  the  base  region  take  somewhat  longer  to  become 
steady,  perhaps  on  the  order  of  milliseconds.  These 
times  are  negligible  from  the  flight  standpoint. 

It  is  well  to  note  here  that  thermodynamic  rela¬ 
tions  and  other  physical  laws,  such  as  the  equation 
of  state  of  a  gas  and  the  law  of  heat  conduction, 
which  strictly  apply  only  to  gases  macroscopically 
“at  rest,”  are  usually  applied  without  alteration  to 
gases  in  motion.  This  procedure  yields  excellent 
results  in  practice. 

6-3  PROCESSES  OF  ENERGY  TRANSFER 
NEAR  A  MISSILE 

A  high  speed  missile  imparts  energy  to  a  limited 
mass  of  air  in  the  vicinity  of  the  vehicle.  Here  "air” 
is  used  to  denote  the  ambient  atmosphere  regardless 
of  its  actual  composition).  Primarily  this  energy 
transfer  takes  place  by  compression  of  the  displaced 
air  and  by  heat-dissipation  due  to  the  shear  stresses 
in  the  boundary  layer.  Part  of  the  compression- 
energy  is  lost  to  the  surrounding  atmosphere  through 
the  dissipative  action  of  shock  waves  which  extend 
outward  from  the  body.  These  shock  losses  increase 
as  the  bluffness  or  pressure  drag  of  the  body  in¬ 
creases,  and  as  the  shock  wave  compression  becomes 
correspondingly  stronger.  In  ballistic  flight  this  in¬ 
crease  in  drag  leads  to  additional  deceleration. 

A  further  loss  cf  energy  to  the  surroundings 
takes  place  as  the  energized  boundary  layer  is  shed 
rearward,  while  a  much  smaller  amount  is  lost  by 
radiation  from  the  hot  gas  surrounding  the  body. 
Energy  also  may  be  gained  or  lost  by  the  missile 
surface  by  radiation  which  is  not  directly  concerned 
with  the  surrounding  air,  for  instance,  by  radiation 
to  the  surface  directly  from  the  sun  or  earth,  and 
by  radiation  from  the  surface  itself  to  the  local 
surroundings.  This  radiation  transfer  is  influenced 
by  the  atomic-scale  characteristics  of  the  air  sur¬ 
rounding  the  missile,  among  other  things,  and  will 
be  discussed  in  a  later  section. 

The  portion  of  the  energy  transmitted  by  the 


missile  to  its  surroundings  and  not  permanently  lost 
to  the  atmosphere  is  returned  to  the  missile  surface 
by  several  possible  heat  transfer  mechanisms  de¬ 
scribed  below. 

6-3.1.  Conduction  and  Convection8  22 


Conduction  is  a  primary  energy  transfer  process 
characterized  by  the  proportionality  between  the 
heat  flux  q  and  the  temperature  gradient  at  a  point. 
Within  the  boundary  layer  the  governing  relation  is 


Btu 


dT 

dy 
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where  q  is  transmitted  along  y,  i.e.,  the  outward 
normal  to  the  surface.  The  thermal  conductivity 
IC  is  a  property  of  this  mode  of  energy  transport 
on  a  molecular  scale.  With  this  sign  convention, 
negative  q  indicates  the  heat  transfered  to  the  body 
from  the  gas. 

Forced  convection  denotes  mechanically-enforced 
(e.g.,  by  pressure  gradients)  fluid  motion  which 
brings  fluid  masses  into  contact  permitting  conduc¬ 
tion  to  take  place  more  readily.  Free  convection 
denotes  fluid  motion  induced  by  body  forces  such 
as  gravity;  it  has  a  negligible  direct  effect  on  heat 
transfer  to  missiles  although  body  forces  may  in¬ 
fluence  the  transition  from  laminar  to  turbulent  flow 
and  thereby  stimulate  increased  heat  transfer  rates. 
A  special  class  of  body  forces  resulting  from  the 
presence  of  electromagnetic  fields  may  influence  the 
motion  cf  conducting  fluids  such  as  ionized  air,  or 
air  artificially  "seeded”  with  materials  which  en¬ 
hance  its  conductivity.  As  a  result,  the  convection 
properties  of  the  fluid  may  be  altered;  these  latter 
effects  fall  into  the  category  of  “magneto-hydro¬ 
dynamics”. 

In  turbulent  flew  random  motion  exists  on  a 
macroscopic  scale  as  distinguished  from  the  molec¬ 
ular  activity  which  leads  to  the  basic  conductive 
transport  previously  described.  Here  turbulent  eddy¬ 
ing  motion  results  in  violent  and  initimate  mixing 
of  small  discrete  fluid  segments  cf  different  momenta 
and  temperature  and  consequently  cause  a  significant 
increase  in  both  the  heat  transport  and  the  effective 
(turbulent)  shear  stress.  By  analogy  with  laminar 
flowfi.e.,  flows  devoid  cf  random  eddies,  attempts 
are  made  to  relate  the  turbulent  heat  flux  to  the 
local  gradient  of  the  mean  temperature  by  an  em¬ 
pirical  proportionality  factor.  This  factor  is  often 
termed  the  turbulent  eddy  conductivity  [kT\.  Unlike 
its  laminar  counterpart  k,  the  turbulent  factor  kT 
is  not  simply  a  physical  property  cf  the  fluid,  but 
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is  quite  sensitive  to  the  nature  of  the  turbulent  flow. 
For  example,  the  value  cf  kT  may  be  quite  different 
in  a  boundary  layer  near  a  surface  than  in  a  free  jet. 
The  analogous  turbulent  shear  parameter,  the  tur¬ 
bulent  eddy  viscosity  [er],  is  sometimes  considered 
to  be  dependent  on  the  local  derivatives  of  the 
mean  velocity  profile  such  that  the  turbulent  shear 
is  considered  to  be  proportional  to  the  square  of 
the  mean-velocity  derivative  rather  than  to  the 
first  power  of  this  derivative,  as  in  laminar  flow. 

6-3.2.  Mass  Diffusion23  33 


In  some  cases  the  air  flowing  over  the  surface 
may  not  be  of  a  homogeneous  nature.  This  occurs 
when  air  is  dissociated  into  a  mixture  of  atoms  and 
molecules  at  high  temperature,  or  when  foreign 
coolants  or  material  ablated  from  the  solid  surface 
mix  into  the  airstream.  In  nonhomogeneous  mixtures 
currents  are  set  up  which  tend  to  promote  homo¬ 
geneity.  The  phenomena  associated  with  this  mixing 
process  are  known  as  "diffusion”  processes.  In  addi¬ 
tion  to  the  diffusion  of  the  various  components 
through  each  other  (i.e.,  "mass  diffusion”)  certain 
secondary  processes  occur.  These  are  known  as 
“temperature  diffusion”  and  "pressure  diffusion.” 
They  are  observed  experimentally  in  nonhomo¬ 
geneous  mixtures  which  are  at  rest  on  the  average; 
when  the  different  species  diffuse  through  each  other; 
and  when  the  molecular  motion  sets  up  temperature 
and  pressure  differences33.  For  example,  the  mean 
diffusion  velocity  of  species  1  in  a  binary  (two 
component)  mixture  can  be  expressed  as  follows: 


Vi  =  -Dtt 


grad  (In  C.)  +  U  '  (r  —  C2  grad  (In  P ) 


+  Cr 1  (jjf)  grad  (In  T)  (6-4) 


where  W t  is  the  molecular  weight  cf  species  i, 

W  =  [Z<VTF,]-\ 

and 

C,  =  p 

An  identical  relation  is  obtained  for  V2  if  C, 
and  Ci  are  interchanged.  The  diffusion  coefficient 
Du  is  virtually  independent  of  concentration. 
Thermal  diffusion  (characterized  by  thermal  diffu¬ 
sion  coefficient  Dr)  tends  to  concentrate  lighter 
and  smaller  species  in  the  warmer  region.  The 
thermal  term  is  usually  less  than  10%  of  the  leading 
term,  while  the  pressure  term  is  negligible  (on  the 
order  of  the  inverse  of  the  momentum-thickness 
Reynolds  number)  compared  to  the  concentration 


gradient  term.23  Thus  only  the  concentration  grad¬ 
ient  term  remains  for  present  application,  and  the 
mass  transport  is  usually  expressed  as  follows  for 
species  1 : 

PC1Vl  =  -PDn  grad  (Cj)  (6-5) 

This  relation  is  known  as  Fick’s  law. 

This  transport  influences  the  gas  composition 
as  well  as  the  flow  properties  at  a  point,  and  thereby 
affects  transport  properties  such  as  thermal  con¬ 
ductivity  and  viscosity  coefficient. 

It  is  noted  that  a  "diffusion  process”  can  be 
defined  as  any  process  in  which  the  net  transport 
of  a  substance  is  directly  proportional  to  the  gradient 
of  a  given  property.  The  proportionality  constant 
is  the  “diffusion  coefficient.”  The  term  “mutual 
diffusion”  describes  the  intermixing  of  2  different 
gases,  while  “self  diffusion”  describes  the  inter¬ 
mixing  of  2  types  (e.g.,  isotopes)  cf  the  same  gas. 

The  energy  absorption  and  release  by  gas  com¬ 
ponents  which  dissociate  and  recombine,  or  which 
react  chemically  with  other  components,  are  linked 
with  the  diffusion  process.  For  instance,  absorbed 
energy  may  be  released  near  a  surface  due  to  the 
diffusion  of  the  reacting  component  across  the  bound¬ 
ary  layer. 

In  most  technical  applications  it  will  suffice  to 
consider  only  the  diffusion  in  binary  (2-component) 
mixtures.  For  example,  dissociated  air  can  be  con¬ 
sidered  for  some  purposes  to  behave  as  a  mixture 
of  2  species:  molecules  and  atoms.  The  differences 
in  the  atomic  weights  of  oxygen  (atomic  weight  16) 
and  nitrogen  (atomic  weight  14)  in  these  cases 
can  be  neglected  within  a  good  approximation.  In 
this  connection  it  will  be  seen  that  binary  diffusion 
can  be  accounted  for  in  a  rather  simple  way  when 
the  “Lewis  number”  (Le  =  pDl2CJk)  which  char¬ 
acterizes  the  effectiveness  of  diffusion,  as  compared 
to  conduction,  in  the  energy  transfer  is  equal  to 
unity. 

We  note  here  that  the  shear  stresses  set  up  by 
the  intermingling  of  species  in  a  mixture  are  negligible 
with  respect,  to  the  usual  viscous  stresses  under  the 
high  Reynolds  number  conditions  of  practical  flight. 
For  a  more  advanced  treatment  cf  chemical  effects 
in  flow  systems  see  References  32  and  33. 

6-3.3.  Heat  Transport  in  a  Binary  Mixture 

Consider  dissociated  air  as  a  binary  mixture 
of  atoms  and  molecules  (species  1  and  2,  respec¬ 
tively).  The  enthalpy  of  each  species  is  then 

h,  =  f  CVi  dT  "I"  h!(  i  =  1,2  (6-6) 

0  • 
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where  h°  is  the  heat  of  formation  of  the  ith  species. 
We  may  take 

for  atoms  h°  =  5-10  ev 

for  molecules  h°2  =  0 

Thus  we  may  write  the  heat  flux  q  in  a  direction 
n  as  follows  (in  practice  n  will  be  the  outward  normal 
to  the  body  surface) : 


+  ,6-7) 


where  F,  is  given  by  Equation  6-5  neglecting  pres¬ 
sure  and  temperature  diffusion,  and  body  force 
differences  on  the  specieg^TJius  we  obtain 
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where 


ture-dependent,  being  proportional,  e.g.,  to  Ts.  Thus 
the  radiation  flux  from  the  gas  can  well  vary  like 
T'2  and  become  appreciable.  Conditions  of  this  kind 
may  be  encountered  during  satellite  and  inter¬ 
planetary  re-entry.34-36 

Although  the  detailed  mechanism  of  radiation 
energy'  transfer  is  quite  complex,  its  engineering 
treatment  may  be  somewhat  simplified.  Radiation 
energy  is  transferred  by  electromagnetic  waves  or 
the  equivalent  quanta.  Most  monatomic  and  di¬ 
atomic  gases  such  as  helium,  hydrogen,  oxygen  and 
nitrogen  are  transparent  to  (i.e.,  emit  but  do  not 
appreciably  absorb)  thermal  radiation,  except  at  very 
high  temperatures;  while  polyatomic  gases  such  as 
water  vapor  and  carbon  dioxide  may  absorb  or  emit 
in  certain  wave  lengths.  A  cogent  introduction  to 
the  role  of  radiative  effects  in  aerodynamics  may 
be  found  in  Ref.  37. 


the  static  enthalpy  h  =  S  Cfli, 

1.2  _ 

the  mass  averaged  specific  heat  Cv  =  ^2,  C,CPl  and 

1 . 2 

the  Lewis-Semenov  number  or  Lewis  number 
Le  =  pDU2C„/k. 

For  Le  =  1,  a  reasonable  approximation  for 
gas  mixtures  in  which  the  molecular  weights  of  the 
components  differ  by  a  factor  cf  only  2  or  3  as 
encountered  in  most  missile  problems,  Equation  6-8 
becomes  simply 


This  relation  is  identical  to  the  expression  for 
the  heat  flux  in  the  simple  case  where  no  dissociation 
or  other  chemical  reactions  occur,  and  where  the 
gas  is  homogeneous.  It  represents  a  major  simplifica¬ 
tion  in  the  analytic  treatment  of  flows  with  chemical 
reactions  since  for  this  case  the  heat  transfer  rate 
is  independent  of  the  mechanism  of  heat  transfer 
and  explicitly  independent  of  chemical  reaction  rates 
to  an  excellent  approximation.  Of  course,  the  fluid 
properties  IC  and  CT  are  dependent  on  the  state  and 
concentration  cf  species  at  a  point.  However,  Equa¬ 
tion  6-9  reduces  the  role  of  the  temperature  (ex¬ 
plicitly)  in  influencing  the  heat  flux  in  flows  with 
chemical  reactions  of  which  dissociation  is  a  special 
case.  It  is  emphasized  that  ionization  (freeing  of 
electrons)  has  not  been  considered  as  yet. 

6-3.4.  Radiation34  37 

Radiation  heat  transfer  rates  are  usually  pro¬ 
portional  to  T*.  However,  at  very  high  gas  temper¬ 
atures  the  gas  emissivity  becomes  sharply  tempera¬ 


6-4  PARAMETERS  FOR  CALCULATION  OF 
CONVECTIVE  HEAT  RATES 

Heat  flux  to  a  surface  can  be  viewed  as  occurring 
under  the  driving  force  of  a  temperature  or  enthalpy 
difference  termed  the  "potential."  The  propor¬ 
tionality  factor  between  the  flux  and  potential  is 
termed  the  heat  transfer  coefficient.  This  concept 
cf  the  heat  flux  is  borne  out  by  theoretical  analysis 
and  experiments  which  are  based  on  the  Prandtl 
boundary  layer  concept,  and  show  that  the  following 
more  definite  quantities  are  of  prime  importance  in 
determining  the  heat  rates. 

6-4.1.  Stagnation  Enthalpy  (hs) 

This  is  the  value  of  enthalpy  which  would  be 
acquired  by  a  moving  air  mass  f  an  adiabatic 
deceleration  to  zero  velocity  were  to  occur.  It  is 
derived  from  the  conservation  of  energy  principle 
which  states  that  under  adiabatic  conditions  the 
sum  cf  the  local  enthalpy  and  kinetic  energy  per 
unit  mass  remains  constant,  i.e., 

h  +  2 JJ  =  h’  (6_10) 

where  J  is  the  mechanical  equivalent  of  heat  (778 
ft-lb/Btu)  and  g  is  the  gravitational  acceleration. 

For  constant  values  of  specific  heat,  Equation 
6-10  becomes 

Y  =  1  +  (t  -  l)M2/2  (6-11) 

For  conditions  in  the  atmosphere  where 
T  500°R  and  h  <  125  Btu/lb,  the  dominant 
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factor  in  determining  h,  at  speeds  above  say  6000 
ft/sec  is  the  right  side  of  Equation  6-10;  that  is. 
the  effect  of  altitude  variations  on  stagnation  en¬ 
thalpy  is  negligible.  This  is  not  the  case  at  lower 
velocities. 

6-4.2.  Adiabatic  Surface  Enthalpy  (haio) 

This  is  the  value  which  hypothetically  would 
occur  at  a  totally  insulated  (zero  net  heat  transfer) 
surface  in  the  absence  of  any  heating  or  cooling 
devices  or  substances  except  those  naturally  in  the 
air.  It  is  the  value  which  that  portion  of  the  boundary 
layer  fluid  m  immediate  contact  with  the  surface 
would  acquire  in  the  absence  of  heat  sources  and 
sinks.  It  differs  from  h,  because  cf  viscous,  conduc¬ 
tive  and  heat  storage  effects  in  the  fluid,  and,  in 
cast'?  of  low-density  or  rarefied  gas  flow,  because 
the  fluid  may  not  adhere  perfectly  to  the  surface. 
It  :?  given  with  excellent  accuracy  by 

—  hr  |  (Pr) 1/2  for  laminar  flow  (6-12) 

A,  —  hr  ~~  (Pr)'/3  for  turbulent  flow 

win  re  Pr  is  the  Prandtl  number.  The  form  cf 
Equation  6-12  for  constant  C„  in  terms  of  T,  is 
envious.  The  quantity  h,  denotes  the  local  enthalpy 
directly  outside  the  boundary  layer. 

6-4.3.  Equilibrium  Surface  Enthalpy  (hC(J) 

The  actual  temperature  of  a  point  on  the  surface 
of  a  body  (and  its  corresponding  enthalpy  in  the 
fluid  directly  adjacent  to  the  surface)  is  established 
by  the  heat  between  the  heat  input  to  the  surface 
by  all  possible  means,  the  heat  re-emitted,  e.g.,  by 
radiation,  and  the  heat  transferred  to  the  interior 
of  the  body.  Under  steady  state  conditions  the  final 
surface  temperature  is  termed  the  “equilibrium 
temperature”  T„  having  a  corresponding  en¬ 
thalpy  hcq. 

For  example,  if  a  coolant  is  injected  into  the 
flow  at  one  point  on  an  otherwise  insulated  body, 
the  surface  acquires  a  temperature  level  and  varia¬ 
tion  different  from  haw.  The  same  is  true  for  a 
completely  insulated  surface  which  loses  heat  by 
radiation. 

6-5  CONVECTIVE  HEAT  RATES  FOR 
CONSTANT  PRESSURE 
AND  TEMPERATURE 

6-5.1.  Two-Dimensional  Laminar  Flow 

The  convective  heat  rates  for  two-dimensional 
laminar  flow  at  a  station  s  along  a  surface  having 
constant  pressure  and  temperature  are  described  by 


the  following  equation: 


(«v)T/2(Pr)I/5 

where 

Nn  =  -t - ^ -  is  the  Nusselt  number, 

~r  ( K  -  K,„) 

vz  pe 

Rk  =  P‘U'S  is  the  Reynolds  number  based  on  S,  and 

(6-13) 

C1 

Pr  =  P".  ”*  is  the  Prandtl  number. 

it  is  often  specified’  that,  for  good  agreement 
with  experiments,  the  “subscript-e”  quantities  in 
Equation  6-13  should  be  evaluated  at  a  “reference 
enthalpy”  given  by 

h*  =  h„  +  0.5(ft„  -  h„)  +  0.22 (h„  -  h„)  (6-14) 

However,  in  Reference  9,  Gold  shows  that,  for 
temperatures  below  the  dissociation  level,  Equation 
6-13  is  virtually  independent  cf  the  evaluation  tem¬ 
perature  and  may  be  expressed  in  a  simplified  form 
as  follows: 

f  ~l5  T  =  2 .00  X  HT6(^)  7  (6-15) 

where  all  quantities  are  in  British  Units.  A  nomo¬ 
graph  for  Equation  6-16  is  given  in  Figure  6-1. 

Furthermore,  shock-tube  experiments  under  condi¬ 
tions  of  dissociated  inviscid  flow  (Reference  31)  have 
shown  that  Equation  6-13  represents  the  data  within 
the  same  order  cf  accuracy  as  the  “reference  en¬ 
thalpy”  method.  Equation  6-13  applies  when  the 
outer  flow  is  not  dissociated  although  the  flow 
within  the  boundary'  layer  may  be  dissociated. 

For  laminar  flow  over  cones  the  constants  in 
Equations  6-13  and  6-15  are  multiplied  by  y/z 
to  give:’ 

(R^'% Prj~2  =  0-55  (616) 

T  ~  T  =  3-56  x  10-eN^Y  2  (6-17) 

Figure  6-1  also  shows  a  nomograph  for  Equa¬ 
tion  0-17. 10 

6-5.2.  Two-Dimensional  Turbulent  Flow 

For  this  type  cf  flow  the  equation  for  convective 
heat  transfer  is 

(R^fer)1'3  =  00296  (6"18) 
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Figure  6-1.  Flat  Plate  Heat  Transfer  Coefficient  (hL),  Laminar  Boundary  Layer" 


In  this  expression  the  "reference  enthalpy"  evalua¬ 
tion  should  be  used  for  best  accuracy.  A  simplified 
form  has  been  derived'  for  Equation  6-18,  but  it 
is  not  independent  of  temperature.  The  simplified 
form  is 

T  =  10.25  X  1CT6  (0-19) 

where  T*  is  the  temperature  obtained  by  replacing 
the  h  values  by  T  in  Equation  6-14.  A  nomograph 
for  Equation  6-19  is  shown  in  Figure  6-2. 10 

For  a  cone  in  turbulent  flow,  the  right  sides 
of  Equations  6-18  and  6-19  must  be  multiplied 
by  1.15." 

It  has  also  been  shown12  that  Equation  6-19 
can  be  applied  locally  as  a  good  approximation  in 
turbulent  flow  over  blunt  bodies  in  hypersonic  flow, 
a  case  of  utmost  importance  in  missile  design. 
Evaluations  of  this  formula  for  a  range  of  local 


surface  pressures,  altitudes,  and  Mach  number  are 
given  in  Figure  6-3.  These  charts  take  into  account 
the  effects  cf  dissociation  equilibrium. 

The  charts  are  expressed  in  terms  of  the  local 
value  of  the  ratio  P/P.,  where  P  is  the  local  surface 
pressure  and  P .  is  the  local  surface  stagnation 
pressure. 

A  critical  point  of  heating  on  a  spherical  nose 
of  radius  R„  over  which  the  flow  has  turned  turbulent 
occurs  in  the  neighborhood  cf  the  sonic  point, 
generally  found  at  an  angular  position  of  about  40° 
from  the  stagnation  point  (i.e.,  at  about  0.7  R„  in 
arc  length  along  the  body  from  the  stagnation  point). 
At  this  point  the  following  rough  estimate  may  be 
deduced  from  the  "local  flat  plate"  formula  for 
hypersonic  conditions. 

qR\  £*  AP°/V,  (6-19a) 

,  it*. 
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where 


A  =  9.6  when  R„  is  in  ft, 

Pt  in  slugs/ft3, 
uf  in  kilo-ft/sec  and 
q  is  in  Btu/ft2-sec. 

In  most  hypersonic  cases  hm/h „  «  1. 

6-6  CONVECTIVE  HEAT  RATE  AT  THE 
STAGNATION  POINT  OF  A 
BLUNT  BODY 

In  practice,  the  stagnation-point  laminar  heat 
rates  are  among  the  most  important  parameters  for 
design.  Since  the  growth  of  the  boundary  layer  is 
initiated  in  the  stagnation  region,  the  flow  will  al¬ 
ways  be  laminar  in  this  area,  at  least  for  a  short 
distance  from  the  stagnation  point  or  line. 

The  stagnation  region  is  an  area  of  maximum 
heat  flux  when  the  flow  over  the  body  is  completely 


laminar.  If,  however,  the  boundary  layer  becomes 
turbulent  at  a  downstream  position  the  heat  flux 
in  some  downstream  area,  e.g.,  near  a  sonic  point, 
may  be  larger  than  at  the  stagnation  point.  This 
peak  turbulent  flux  may  be  on  the  order  of  twice 
the  stagnation  point  laminar  value.  Nevertheless, 
the  calculation  of  the  heat  flux  to  a  blunt  body 
always  starts  with  the  laminar  stagnation  point 
computation.  This  is  so  for  two-dimensional  (wing- 
like)  bodies ;  axially-symmetric  bodies  (bodies  of 
revolution  at  zero  angle  of  attack);  bodies  with 
almost  spherical  noses  even  when  at  an  angle  of 
attack  (provided  the  stagnation  point  remains  on 
the  sphere) ;  and  nonaxially  symmetric  bodies  (such 
as  those  with  ellipsoid-like  noses).  In  nonaxially 
symmetric  cases  the  heat  rate  is  not  single-valued 
at  the  stagnation  point,  but  takes  on  different  values 
along  each  of  the  streamlines  emerging  from  the 
stagnation  point. 


t  =1.15  (h-^) 
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Figure  6-2.  Flat  Plate  Heat  Transfer  Coefficient  (hT),  Turbulent  Boundary  Layer" 
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( D )  Maximum  Heating  Rates  as  a  Function  cf  Flight 
Mach  Number  for  Various  Altitudes  and  Wall 


(B)  Altitude:  100,000ft  Tw  =  500° R  Temperatures 

Figure  6-3.  Turbulent  Heat  Transfer  Rates  as  a  Function  oEStatic-to-Stagnation-Pressure  Ratio 


The  stagnation  region  heat  rates  depend  strongly 
on  the  magnitude  of  the  velocity  gradient  ( dujds ) 
along  streamlines  emerging  from  the  stagnation 
point.  In  fact,  they  are  directly  proportional  to 
(dujds)1/2.  Thus  for  nonaxially  symmetric  bodies 
the  heat  rates  will  change  along  the  various  stream¬ 
lines  as  dujds  changes. 

A  very  accurate  relation  for  the  laminar  stagna¬ 
tion  point  heat  rate  had  been  obtained  theoret¬ 
ically23,31  and  verified  experimentally.  This  relation 
takes  into  account  the  “real  gas”  effects  of  molecular 
dissociation  in  the  two  limiting  cases  of  “very  fast” 
and  “very  slow”  rates  of  dissociation  and  brackets 
the  values  for  intermediate  reaction  rates.  It  is 
given  by 


(h..  -  hw) 


where  hD 


\pSMa  / 


X  1  +  ( Lem 


(tr 


(6-20) 


51  CJ  —  hi)  is  the  chemical  formation 


also  suggests  that  this  value  of  A  may  probably 
apply  to  the  atmosphere  of  Venus  and  Mars  as  well 
as  the  Earth. 

As  previously  mentioned.  Equations  6-20  and 
6-21  also  apply  for  the  stagnation  point  or  line  of 
bodies  having  nearly  spherical  or  nearly-cylindrical 
noses  at  angles  of  attack  provided  the  stagnation 
region  remains  on  the  nose,  i.e.,  for  moderate  angles 
of  attack. 

A  nomograph  giving  the  stagnation-point  heat 
transfer  rates  for  ballistic  re-entry  of  a  body  with 
an  approximately  spherical  nose  is  given  in  Fig.  6-4. 12 

More  recently,  Hoshizaki31,  presented  a  simple 
relation  for  determining  the  heat  transfer  rate  at 
an  axisymmetric  stagnation  point  in  the  velocity 
range  of  6  to  50  kilo-feet/sec.  This  expression  in¬ 
cludes  the  effect  of  variable  Prandtl  and  Lewis 
number  as  well  as  ionization,  for  equilibrium  chem¬ 
istry,  and  has  the  following  form 


q,  a/AL  =  C„,> 


X  du, "1 1  /2  vJuf  \3A Y,  _  K\ 
,Uf  ds  _  Pr  \1()4/  \  hj 


enthalpy  produced  by  dissociation  of  molecules  and  where 


=  0  for  two-dimensional  flow 
=  1  for  axially-symmetric  flow'. 

Here  the  constant  m  =  0.52  for  very  rapid  dissocia¬ 
tion  (i.e.,  for  the  case  of  chemical  equilibrium)  and 
m  =  0.63  for  very  slow  dissociation  (i.e.,  for  the 
case  of  constant  or  “frozen”  chemical  composition) 
with  a  catalytic  surface.  According  to  the  modified 
Newtonian  outer  flow  approximation,  the  velocity 
gradient  in  the  stagnation  region  is 

iu  =  £-1  2(p,  -  Vf ) 
ds  L  p’ 


For  values  of  Le  1.4,  as  is  the  case  for 

dissociated  air,  the  term  in  square  brackets  in 
Equation  6-20  takes  on  a  maximum  value  of  about 
1.2.  Thus,  Equation  6-20,  which  was  derived  in  a 
very  complex  fashion,  agrees  quite  closely  with  the 
earlier  results  of  more  approximate  calculations. 

A  simplified  engineering  form  for  Equation  6-20 
would  be31 


For  R  in  feet,  the  flight  conditions  p,  and  u, 
in  slugs/ft3,  and  in  kiloft/sec,  respectively,  and  q,  in 
Btu/ft2-sec,  the  constant  A  is  equal  to  15.5. 

The  value  of  A  is  4710  when  R ,  is  in  cm,  p,  in 
gm/cm3,  u;  in  km/sec,  and  q,  in  watts/cm2.  Lees 


Cair  =  2.59  X  104  Btu-ft'/2/lb1/2-sec2 

A  comparison  of  experimental  and  theoretical 
results  for  stagnation  point  heat  transfer  at  high 
velocities  is  given  in  Figure  6-5. 

A  nonspherical  stagnation  region  with  two  planes 
of  symmetry  may  be  approximated  using  an  ellip¬ 
soidal  local  body  geometry.  Thus,  Equations  6-20 
and  6-21  may  be  applied  in  this  case  with  slight 
modifications. 13  For  this  approximation,  the  value 
cf  dujds  in  Equation  6-20,  or  the  corresponding 
value  of  R.  in  the  Newtonian  approximation,  must 
be  interpreted  as  the  velocity  gradient  and  radius 
cf  curvature  corresponding  to  the  plane  of  the 
smaller  principle  radius  of  curvature,  i.e.,  (dujds)l 
and  R„.u  respectively.  On  the  right,  side  of  Equation 
6-20  the  factor  (2)’l/2  must  be  replaced  by  (1  '  B)l/2 
where  B  is  the  ratio  of  the  two  perpendicular  velocity 
gradients,  i.e.,  B  =  (dujds) J (dujds) , .  Similarly 
in  the  Newtonian  approximation  [Equation  6-21] 
B  =  R„,JR„,2.  Note  that  for/?„i2  =  Rn, i,  B  =  1, 
thus  giving  the  axially  symmetric  case;  while  for 
R„,i  —*■  ,  B  — >  0  giving  the  two-dimensional  case. 

6-7  LAMINAR  HEAT  TRANSFER  ALONG  A 
STREAMLINE  OF  A  BLUNT  BODY 

With  good  approximation,  particularly  for  high 
supersonic  and  hypersonic  flow  conditions  over  blunt 
bodies,  we  may  consider  the  flow  along  an  inviscid 
streamline  emerging  from  the  stagnation  region  as 
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if  that  streamline  were  everywhere  a  local  body  of  be  shown  to  hold  for  the  aforementioned  cases 
revolution.  The  pressure  variation  along  this  stream-  w  ithin  good  engineering  accuracy.23,24 
line  must  be  known  apriori,  of  course.  Then  the 

following  ratio  of  local  heat  transfer  to  the  stagna-  6-8  MASS  TRANSFER  AND  ABLATION 
tion  point  heat  rate  can  be  obtained.  COOLING 


r-w  r  f  (^Won-<feT 

\as/,  LJo  Vp.M,/  J 


where 


(6-22) 


r0  =  radius  of  cross  section  cf  bodies  of  revolution 
and 

n  =  1  for  bodies  which  are  like  bodies  of  revolu¬ 
tion. 

Equation  6-22  was  derived  by  Lees14  strictly 
for  the  case  of  axially  symmetric  or  two-dimensional 
flow  assuming  “local  similarity.”  However,  it  can 


The  “first  generation”  cf  ballistic  missiles  have 
generally  employed  “heat-sink”  methods  cf  main¬ 
taining  sufficiently  cool  structural  temperatures. 
That  is,  a  mass  cf  metal,  otherwise  structurally 
extraneous,  was  applied  to  the  body  usually  in  the 
form  cf  a  sheath.  Theis  material  absorbed  the  heat 
generated  during  flight  and  thereby  maintained 
tolerable  surface  temperatures.  Copper  was  found 
to  be  a  particularly  desirable  heat  sink  material 
because  cf  its  conductive  and  absorptive  properties. 

An  upper  limit  to  the  utility  of  metallic  heat 
smks  is  given  by  their  maximum  absorptive  capacity 
prior  to  melting  or  structural  deterioration.  For 
copper  this  absorptive  capacity  is  on  the  order  of 
400  Btu/lb.  For  missile  performance  requiring  total 
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heat  inputs  over  the  trajectory  greater  than  several 
hundred  Btu/lb  other  cooling  methods  than  pure 
absorption  must  be  utilized. 

In  the  “second  generation”  missiles  which  re¬ 
quire  higher  impact  velocities  than  can  be  achieved 
with  simple  nondestructive  heat-sink  methods,  ma¬ 
terials  have  been  utilized  which  melt  or  evaporate 
in  a  relatively  orderly  fashion.  These  materials, 
primarily  plastic  in  nature  (e.g.,  bonded  resins  such 
as  Fiberglas,  or  materials  like  Teflon)  have  been 
used  as  coatings  or  external  caps  on  the  nose  cone. 
Their  latent  heat  of  melting,  vaporization  or  sub¬ 
limation  augment  their  effective  heat  absorptive 
capacity,  increasing  it  to  the  order  of  several  thou¬ 
sand  Btu/lb.  These  materials  are  expendable  and 
their  removal  due  to  thermal  effects  during  flight  is 
termed  “ablation.” 

Of  course,  ablating  materials  must  have  not 
only  desirable  absorptive  capabilities  but  also  must 
have  desirable  structural  properties  to  resist  surface 
shear  by  the  air  and  other  aerodynamic  loads. 

An  additional  desirable  characteristic  of  abla¬ 
tion  cooling  is  derived  from  the  fact  that  the  material 


which  is  sublimed  or  evaporated  from  the  surface 
acts  as  a  gaseous  coolant  blocking  the  usual  entry 
of  heat  to  the  surface,  such  as  if  a  cooling  fluid  were 
being  blown  out  cf  the  surface  itself.  Also  accom¬ 
panying  this  process  may  be  chemical  reactions 
between  the  gaseous  ablated  material  and  the  high 
temperature  air  around  the  body.  It  has  been  found 
that  these  chemical  reactions  may  actually  be 
beneficial  for  the  cooling  process. 

Naturally,  these  complicated  processes  are  some¬ 
what  difficult  to  treat  analytically.  However,  some 
progress  has  been  made  toward  an  understanding 
and  simplified  theoretical  treatment  of  boundary 
layer  flows  subjected  to  mass  addition  and  chemical 
reaction.  (References  22,  23,  24) 

Both  theory  and  experiments  indicate  that 
Teflon  placed  in  the  stagnation  region  of  a  spherical 
nose  in  hypersonic  flow  has  an  effective  absorptivity 
or  “effective  heat  of  ablation”  which  varies  linearly 
from  about  3,000  Btu/lb  when  the  flight  velocity 
is  15,000ft/sec  to  about  5,000 Btu/lb  at  a  velocity 
of  20,000  ft/sec.  This  dependence  upon  velocity  and 
position  of  the  material  on  the  body  reflects  the 
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complex  fluid-mechanical  and  physico-chemical  proc¬ 
esses  which  continue  to  give  the  engineering  quantity 
termed  the  “effective' heat  of  ablation”  denoted 
by  Htfr. 

For  engineering  purposes  one  estimates  the 
mass  lost  from  an  ablating  material  by  computing 
the  theoretical  heat  input  to  the  body  as  if  there 
were  no  mass  addition  due  to  material  loss  (i.e.,  the 
value  q  as  previously  discussed  in  this,  chapter). 
Then  the  quantity  q/H.,,  yields  the  lb  of  material 
lost  per  ft2  of  exposed  surface  per  second  during 
the  course  cf  exposure. 

For  these  purposes  it  is  of  interest  to  note  that 
the  order  of  magnitude  of  the  average  heat  input 
to  the  entire  surface  cf  a  spherical  nose  is  character¬ 
ized  by  the  stagnation  point  heating  value  when 
the  flow  is  laminar  and  by  the  sonic  point  value 
when  the  flow  is  turbulent.  Integration  of  the  heat 
transfer  rate  over  a  hemispheric  surface  at  a  given 
instant  of  time  shows  a  mean  heat  rate  equal  to 
about  5  the  peak  heat  rate.  On  pointed  or  slender 
bodies  with  large  amounts  of  lateral  surface,  the 
flat-plate  or  cone  heat  rate  levels  are  more  charac¬ 
teristic  of  the  overall  heat  inputs  to  the  body  than 
are  the  local  peak  values  which  occur  at  the  tip 
and  which  may  require  special  local  cooling  methods 
for  the  preservation  thereof. 

An  idea  of  the  theoretical  methods  by  which 
the  actual  heat  transfer  to  an  ablating  material  may 
be  studied  in  detail,  and  whereby  the  effective  heat 
of  ablation  may  be  estimated,  is  given  in  the  para¬ 
graphs  which  follow. 

A  simplified  expression  for  oxygen-and-nitrogcn- 
limited  chemical  reactions  with  a  pure  ablating  sur¬ 
face  material  E  is  treated  in  detail  by  Lees24'25. 
Here  the  salient  points  of  interest  are  outlined.  If 
q,  is  the  net  rate  of  convective  heat  transfer  to  the 
solid  interior  available  for  heating  fresh  surface 
material,  where  subscript  s  denotes  solid  and  g 
denotes  gas,  then  for  reactions  occurring  between 
atomic  and  molecular  oxygen  and  nitrogen  and  a 
surface  material  E  (here  an  element)  there  results 
the  expression 


q.  = 


i  dT  .  „  y  ,  dCj 

^  dn  +  p°12  ?  hi  dn 


-  (pv)whjg) 


+  (pv)*,hE,w(s). 


(6-23) 


For  pDl2CJk  =  1  this  yields 

q,  =  P'U'CH(Ah  —  B'Le)  (6-24) 


Ah  =  £  -  hiiW) 

1 

+  4+  I lh{.„[Ct.,  -  (B'  +  1)C,.„] 

+  [S'  -  (B'  +  1  )CE,w]hE.„{g)  (6-25) 


__  (pr’)M 
p'UeCH 


(6-26) 


and  Le  is  the  latent  heat  of  sublimation  per  gram 
cf  E.  For  oxygen-nitrogen-limited  reactions,  wherein 
the  mass  fractions  of  individual  chemical  elements 
C  are  introduced  rather  than  each  species  C„  this 
expression  may  be  simplified  to 

q,  =  P'U'Cn{h,,t  Cj.thj,w 

j 

+  YICJ.,AQE.J-  B’Le\  (6-27) 


where  AQ  denotes  the  heat  released  in  a  chemical 
reaction  per  unit  mass,  and 

V  Wt 

where  mj  ,*  is  the  number  of  atoms  (or  molecules) 
cf  the  element  J  in  the  species  k. 

For  the  complete  combustion  of  carbon  by 
oxygen  and  nitrogen  the  added  chemical  enthalpy 
potential  is  about  24,000  Btu/lb,  about  twice  the 
potential  without  combustion.  However,  this  is  re¬ 
lieved  by  the  heat  absorption  due  to  sublimation, 
about  28,500  Btu/lb,  about  ten  times  that  of  un¬ 
burnt  carbon.  Thus  the  net  effect  is  not  catastrophic 
but  may  be  somewhat  beneficial  in  reducing  heat 
inputs. 

The  “heat  blockage”  or  reduction  in  Cff  brought 
about  by  the  mass  transfer  can  be  regarded  as  an 
increase  in  effective  heat  sink  capacity,  i.e., 

Erff  _  (I  H.n 

Le  +  Lr  C  u 

where  CHi 0  is  the  Stanton  number  for  zero  mass 
transfer.  The  ratio  CH,0/CH  depends  only  on 

g,  _  Ah,,f 

Le  +  Lt 


as  given  by  Fig.  6-6.  Here  LT  is  the  thermal  capacity 
cf  the  solid  material  from  some  interior  temperature 
Ti  up  to  the  surface  temperature  T„,  i.e., 

T  =  g.  ~  <?S 

(pv)w 


where  CH  is  the  Stanton  number, 


where  qR  is  the  net  radiative  heat  transfer  away 
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Figure  6-6.  Effective  Heat-Sink  Capacity  with  Mass 
Addition  '5 


from  the  wall.  Furthermore,  the  term  A he!t  is  an 
"effective"  enthalpy  potential  defined  by 


6-9  RADIATION  HEAT  TRANSFER 

The  einissivity  £  for  air  considered  as  an 
optically  thin  layer  (e  <  l)has  been  tabulated  in 
Ref.  39  and  approximately  correlated  in  Ref.  24 


as  follows: 


where  L  is  a  characteristic  length  cf  the  body 
(in  cm)  and  T ,,  is  the  stagnation  temperature 
external  to  the  boundan  layer  (in  °K).  If  chemical 
equilibrium  is  assumed,  one  estimates,  for  the  flight 
regimes  cf  interest: 


—  ~  17,  T„  ~  0.32w./ 
Pr 


where  uf  is  the  flight  velocity  (in  ft/sec). 

With  an  approximate  shock  layer  geometry  and 
data  for  equilibrium  states  from  Fief.  39,  estimates 
of  the  heat  flux  to  the  stagnation  region  of  a  blunt 
body  have  been  made  in  Ref.  3,5.  The  graphical 
representation  of  such  estimates  is  presented  in 
Fig.  6-7.  These  results  may  be  correlated  adequately 
by  the  expression  developed  in  Ref.  40.  In  this  case, 
i.e.,  e  <  1,  the  heat  flux  due  to  radiation  is  given  by 


(\L.60/  \  8 . 5 

—£L—)  (t?f) 

P sea  level'  \1U/ 


where 


qR  is  the  radiative  heat  flux  (in  Btu/ft2  —  sec), 
R„  is  the  nose  radius  (in  ft)  and 
pf  is  the  density  at  flight  altitude. 


Figure  6-7.  Radiative  Heating  as  a  Function  of  a  Flight  Altitude  and  Velocity™ 
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However,  if  the  air  is  not  in  chemical  equilibrium 
in  the  shock  layer,  the  local  temperatures  may  be 
much  higher  than  those  obtained  assuming  chemical 
equilibrium.  As  an  approximation,  a  first  estimate 
of  the  radiant  heat  flux  can  be  made  if  the  non- 
equilibriumt  emperature  is  known,  by  evaluating 
the  expression  for  e/L  given  above  at  the  non¬ 
equilibrium  state  and  using  the  relation  qR  ~  eaT\ 
where  a  is  the  Stefan-Boltzmann  constant. 

Additional  work  relative  to  radiative  heat 
transfer  effects  is  given  in  References  29,  36  and  37. 


6-10  CONCLUDING  REMARKS 

There  are  numerous  special  problems  concerning 
the  heat  transfer  to  missiles  which  have  not  been 
treated  explicitly  in  this  Chapter.41  For  example, 
the  augmented  heat  rates  which  may  be  engendered 
by  protuberances  from  a  surface  have  been  con¬ 
sidered  in  Reference  42,  and  the  heat  transfer  and 
pressures  over  the  blunt  base  of  a  body  has  been 
studied  in  Reference  43.  The  reader  should  consult 
the  references  for  the  solution  of  other  heat  transfer 
problems. 
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CHAPTER  7 

HYPERSONIC  AERODYNAMIC  FORCES 


7-1  INTRODUCTION 

Considerable  information  has  been  compiled 
concerning  forcesand  moments  on  bodies  in  sub¬ 
sonic  and  supersonic  flight.  Much  of  this  information 
concerns  wings,  wing-like  shapes,  and  sharp-nosed 
slender  bodies  in  the  low  or  moderate  supersonic 
range,  e.g.,  below  Mach  4.  For  the  purpose  of  this 
discussion,  the  speed  range  below  roughly  Mach  4 
will  be  termed  “low-speed”  and  the  range  above 
Mach  4  will  be  referred  to  as  “hypersonic.” 

The  utility  of  low  speed  aerodynamic  properties 
in  connection  with  the  consideration  of  modern 
ballistic  missiles  is  somewhat  limited  inasmuch  as 
low  speed  flight  regimes  are  confined  to  the  boost 
phase  and  the  final  50,000  ft  (roughly)  of  re-entry. 
In  the  boost  phase  the  vehicle  weight  and  engine 
thrust  far  exceed  the  aerodynamic  forces;  hence 
the  aerodynamic  forces  are  of  secondary  importance 
for  preliminary  analyses.  This  is  not  to  say  that 
they  do  not  have  an  important  role  in  the  final 
phase  of  design,  especially  in  connection  with  stabil¬ 
ity  and  control.  However,  the  considerations  in¬ 
volved  are  of  a  relatively  detailed  nature,  and  may 
be  thought  of  as  lying  within  the  scope  of  more 
‘  ‘conventional”  aerodynamics. 

In  the  final  phase  of  re-entry  which  may  involve 
low  speeds  if  the  re-entry  vehicle  is  sufficiently  light, 
the  details  of  forces  do  not  have  a  great  influence 
on  preliminary  design  since  their  dispersive  effects 
cannot  be  great  and  since  critical  regions  cf  peak 
heating  and  deceleration  occur  well  above  the  low 
speed  part  cf  the  trajectory. 

Low-speed  problems  therefore  will  not  be  con¬ 
sidered  explicitly  in  this  chapter.  Properties  of  this 
regime  are  developed  in  sufficient  detail  in  standard 
texts  and  in  applied  design  works  such  as  References 
1  and  2.  This  chapter  will  emphasize  instead  the 
hypersonic  flight  conditions  unique  to  ballistic  mis¬ 
sile  design  and  performance.  Similarly,  we  shall  not 
emphasize  the  aspect  of  forces  relating  to  lift  (the 
force  component  normal  to  the  flight  velocity  vector) 
since  lift  is  employed  primarily  in  satellite  or  manned 
re-entry  to  minimize  heating  and/or  deceleration  and 
to  facilitate  terminal  guidance.  Furthermore,  wing¬ 


like  surfaces  which  are  usually  used  solely  to  provide 
lift  or  control,  assume  less  importance  in  hypersonic 
flight,  since  their  contribution  to  lift  is  dominated 
by  that  of  the  body  itself.  Thus  the  body  shapes 
which  will  be  of  principle  interest  will  be  spheres, 
cones,  sphere-cones  and  sphere-cone-cylinder-flare 
combinations. 

An  additional  simplification  arises  from  the  fact 
that  for  many  configurations  the  portion  of  the  drag 
(force  component  parallel  to  the  flight  velocity 
vector)  which  stems  from  viscous  skin  friction  is 
comparatively  small  and  may  be  neglected  or  roughly 
estimated.  This  is  not  the  case  for  very  slender 
bodies,  such  as  cones  with  included  angles  on  the 
order  of  10°  or  less,  for  which  skin  friction  drag 
becomes  important. 

The  protion  of  the  drag  not  contributed  by  skin 
friction  is  termed  pressure  drug.  This  drag  is  de¬ 
termined  by  integrating  the  pressures  over  the  sur¬ 
face  and  taking  the  component  of  the  resultant 
parallel  to  the  flight  velocity  vector.  The  body 
pressures,  aside  from  those  over  the  base,  which 
have  a  minor  effect  on  hypersonic  forces  are  estab¬ 
lished  essentially  by  the  ilow  external  to  the  bound¬ 
ary  layer.  In  this  external  region  viscosity  may  be 
neglected.  Subclasses  of  pressure  drag,  termed  wave 
drag  and  base  drag,  need  not  be  distinguished  here. 

When  the  ratio  of  pressure  drag  to  friction  is 
large,  the  body  is  termed  bluff;  the  opposite  type 
of  body  is  termed  streamlined.  Viscous  and  pressure 
force  components  normal  to  the  flight  velocity  vector 
comprise  the  lift.  This  characteristic  is  dominated 
by  the  pressure  forces. 

Within  the  current  state  of  the  art  it  is  quite 
possible  to  make  detailed  theoretical  calculations  of 
the  hypersonic  flow  fields  around  missile  configura¬ 
tions,  even  including  the  effects  of  flow  chemistry. 
However  the  procedures  required  are  highly  complex 
and  involve  numerical  integrations  on  large  digital 
computers. 

I  n  fact,  such  calculation  programs  are  considered 
to  be  “numerical  experiments”  and  are  used  in  touch¬ 
stone  cases  to  check  approximate  theories,  or  in  the 
final  stages  of  design.  Examples  of  such  calculations 
are  given  in  References  3-6.  These  deal  with  bodies 
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at  zero  angle  of  attack,  i.e.,  axially-symmetric  flow. 
The  cases  of  11011-axial  symmetry  which  arise  due 
to  angle  of  attack  are  much  more  difficult  to  obtain, 
and  at  present  are  at  the  limit  of  available  tech¬ 
niques.  Fortunately,  however,  the  simpler  problem 
of  estimating  surface  forces  in  hypersonic  flight  does 
not  require  the  aforementioned  detail  since  it  is 
amenable  to  approximation.  Appropriate  engineering 
approximations  for  the  resultant  forces  and  moments 
form  the  main  subject  of  the  subsequent  paragraphs 
of  this  chapter. 

7-2  STEADY  FRAME  OF  REFERENCE 

Aerodynamic  effects,  such  as  forces  acting  on 
a  body  in  flight,  are  considered  in  terms  of  an  equiv¬ 
alent  system  consisting  of  a  fixed  body  subjected 
to  an  initially  undisturbed  stream  moving  relative 
to  the  body.  This  is  the  situation  actually  produced 
in  wind  tunnel  testing.  The  initial  flow  is  taken  to 
have  a  velocity  vector  u,  equal  in  magnitude  and 
direction,  but  opposite  in  sense  to  the  vehicle’s  flight 
velocity  vector.  The  static  pressure  pf,  density  pf 
and  temperature  Tf  of  the  initial  stream  are  those 
of  the  undisturbed  atmosphere  at  the  flight  altitude. 
In  this  frame  of  reference  the  flow  relative  to  the 

*  Here  subscript  /  denotes  freestream  flight  conditions. 


body  is  steady;  i.e.,  at  each  point  there  are  no  time- 
wise  variations.  Hence  the  determination  of  flow 
effects,  such  as  forces,  is  simplified.  The  usual  varia¬ 
tions  cf  flight  conditions  along  a  trajectory  are  taken 
into  account  by  a  “quasi-steady”  approximation 
which  considers  a  succession  of  instantaneously 
achieved  steady-states  along  the  flight  path.  Special 
conditions,  involving  extreme  curvatures  of  tra¬ 
jectory,  extreme  time-variations  in  flight  condi¬ 
tions,  or  vibrating  surfaces,  for  which  the  quasi¬ 
steady  approach  may  not  be  valid  are  too  limited 
for  consideration  here. 

7-3  QUALITATIVE  STRUCTURE  OF  THE 
HYPERSONIC  FLOW  FIELD 

Use  of  the  approximate  methods  of  force  estima¬ 
tion  to  be  given  here  does  not  require  a  detailed 
knowledge  of  hypersonic  flow  fields,  whose  complex 
structure  can  be  observed  from  the  typical  schematic 
diagram  shown  in  Figure  7-1  (from  Ref.  6).  However, 
several  salient  features  are  cited  to  serve  as  an  in¬ 
troductory  background. 

Shock  waves  are  formed  in  supersonic  flows, 
wherein  local  velocities  exceed  local  sound  speeds, 
due  to  the  coalescence  of  disturbances  which  cannot 
propagate  away  from  the  body  more  rapidly  than 
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Figure  7-1.  Schematic  Diagram  of  Complete  Flow  Field  Around  a  Blunt-Nosed  Body 6 
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the  oncoming  flow  sweeps  them  back.  As  the  flow 
traverses  the  shock  its  velocity  decreases,  whereas 
the  static  pressure,  density  and  temperature  in¬ 
crease.  In  hypersonic  flow  the  bow  shock  lies  very 
close  to  the  body  and  is  roughly  parabolic  in  shape 
downstream  of  the  nose.  The  region  between  the 
shock  and  the  body  is  termed  the  shock  lager.  Its 
thickness  in  the  nose  region  cf  a  blunt  body  in 
hypersonic  flow  is  on  the  order  cf  |  cf  the  nose  radius. 

The  boundary  layer  is  a  region  adjacent  to  the 
body  wherein  significant  shear  stresses  are  created 
by  the  adhesion  of  the  fluid  to  the  surface  (no-slip 
condition).  The  surface  shear  stress  has  the  rough 
magnitude  r„  ~  pw(2ue/8),  where 

a i„  is  the  absolute  viscosity  evaluated  at  the  sur¬ 
face  temperature, 

is  the  local  velocity  outside  the  boundary  layer 
and 

8  is  the  local  thickness  of  the  boundary  layer. 

The  ratio  2 ujh  signifies  the  order  of  magnitude  cf 
the  velocity  gradient  or  rate  of  strain  within  the 
layer. 

The  thickness  6  has  the  following  order  cf 
magnitude  at  a  distance  x  away  from  the  initial 
point  f  the  flow  remains  laminar: 


A  ~  6  +  (?  —  -)a/J 

where 

7  is  the  ratio  cf  specific  heats  with  typical  values 
between  1.2  and  1.4  in  air, 

M  is  the  Mach  number  and  subscript  e  denotes 
values  outside  the  boundary  layer  at  the  sta¬ 
tion  x. 

The  term  in  parentheses  is  a  dimensionless  group 
called  the  Reynolds  number  ( RN ),  or  in  this  instance 
the  “local”  Reynolds  number  because  its  parameters 
are  evaluated  at  the  local  conditions  indicated. 

The  local  density  pe  varies  in  rough  proportion 
to  the  ambient  density  pf,  decreasing  from  about 
10  pf  to  2 pf  around  a  typical  blunt  body.  The  local 
velocity  m,  increases  from  zero  at  the  nose  cf  the 
blunt  body  to  about  0.5  uf  at  the  end.  The  temper¬ 
atures  which  determine  n,  are  more  difficult  to 
evaluate  generally  because  of  their  sensitivity  to 
chemical  composition  ;  however,  they  may  decrease 
from  levels  around  20  to  100  T,  in  the  nose  region 


cf  a  blunt  body  to  10  T  at  the  rear.  Thus  p,ur,/p.r  ~ 
(^PfUf/p.,  and  8/ x  ~  1.5 (pfitfx/nf)~l/2.  At  an  altitude 
cf  150,000  ft,  flight  Reynolds  numbers  per  foot 
(PfUf/pt)  are  on  the  order  cf  105.  Therefore  for 
x  =  10ft  a  typical  thickness  6~  0. 15  ft  is  estimated. 
At  200.000  ft  8  increases  to  about  0.5  ft,  and  at 
100.000ft  it  decreases  to  0.05  ft. 

The  foregoing  rough  estimates  indicate  that  at 
low  altitudes  (e.g.,  below  150.000ft)  the  boundary 
layer  thickness  is  relatively  small  compared  to  char¬ 
acteristic  body  dhnensions  and  to  the  shock  layer 
thickness,  whereas  at  higher  altitudes  the  boundary 
layer  may  become  sufficiently  thick  to  occupy  a 
major  part  of  the  shock  layer. 

f  rom  the  standpoint  of  forces,  a  key  feature 
cf  the  boundary  layer  is  the  fact  that  it  is  virtually 
parallel  to  the  surface,  and  thus  transmits  the  pres¬ 
sure  at  its  outer  edge  essentially  undistorted  to  the 
surface  along  a  surface-normal.  As  a  result,  the  effect 
cf  the  boundary  layer  thickening  on  forces  is  two¬ 
fold:  First,  the  shear  stress  is  decreased  so  that  at 
high  altitudes  its  contribution  is  smaller  than  at 
low  altitudes.  Secondly,  its  displacement  effectcauses 
an  effective  thickening  cf  the  body  with  respect  to 
the  flow  outside  the  boundary  layer.  It  has  been 
shown  both  theoretically  and  experimentally8  that 
this  interaction  has  only  a  modest  effect  on  pressures 
over  bodies  of  revolution.  For  example  on  10°  half 
angle  cones,  pressures  about  20%  higher  than  those 
obtained  with  thin  boundary  layer  are  found.  Al¬ 
though  the  effect  is  much  more  pronounced  on  two- 
dimensional  wing-like  surfaces,  the  effect  cf  viscous- 
inviscid  interactions  on  larger  angle  cones  and  on 
blunt  bodies  of  revolution  is  diminished  still  further. 
Therefore,  we  are  justified  in  neglecting  this  low 
Reynolds  number  phenomenon  for  present  purposes. 
Indeed  at  the  high  altitudes  where  the  interaction 
would  be  observed,  the  overall  aerodynamic  forces 
are  small,  and  have  only  a  minor  influence  on  ballistic 
missile  behavior. 

A  w  ord  may  be  said  about  the  separated  flow  or 
base  region  cf  the  body  surface;  denoted  as  the 
recirculation  region  in  Figure  7-1.  When  the  flow 
outside  the  boundary  layer  creates  regions  of  in¬ 
creasing  pressure  over  the  boundary  layer  (and  body 
surface),  the  boundary  layer  is  usually  forced  to 
separate  from  the  body  and  thus  creates  an  eddying 
wake  region.  On  wings,  this  condition  is  known  as 
stall  and  leads  to  increased  drag  and  decreased  lift. 
On  blunt-based  bodies  cf  revolution,  the  cause  cf 
separation  is  somewhat  more  complex,  being  related 
to  the  ultimate  redeflection  cf  the  flow  beyond  the 
body.  Estimates  cf  the  mean  pressure  in  the  re- 


7-3 


circulation  region,  i.e.,  the  base  pressure,  are  usually 
made  by  semi-empirical  methods.  In  low-speed  flow 
the  contribution  of  the  base  force  to  the  overall 
forces  on  the  body  is  significant  and  requires  accurate 
evaluation.  However,  In  hypersonic  flight  the  base 
force  usually  is  negligible  and  therefore  will  not  be 
considered  here. 

Finally,  it  may  be  remarked  that  the  boundary 
layer  flow  may  be  laminar  or  turbulent.  Usually 
turbulent  boundary  layers  will  not  be  found  at  the 
higher  altitudes  (above  150,000ft).  However,  at  low7 
altitudes  they  are  almost  certain  to  arise.  Estimates 
of  where  the  transition  from  laminar  to  turbulent 
flow  will  occur  are  extremely  difficult  to  make.  For 
this  purpose  one  usually  relies  on  experiments  or 
prior  experience  in  similar  cases  and,  when  in  doubt, 
makes  alternative  estimates  based  on  both  laminar 
and  turbulent  possibilities.  The  main  effect  cf  tur¬ 
bulence  on  the  force  is  to  increase  the  surface  shear 
and  thus  the  skin  friction  drag  contribution.  It  also 
leads  to  boundary  layers  which  are  thicker  than 
laminar  ones. 

Textbooks  dealing  in  detail  with  the  subject 
matter  cf  the  chapter  are  listed  as  Refs.  7,  10  and  11. 

7-4  FORCE  ESTIMATES  BY  THE 

NEWTONIAN  IMPACT  CONCEPT 

The  Newtonian  impact  concept,  previously 
cited  in  Chapter  5  in  connection  with  surface  pres¬ 
sures.  provides  a  very  useful  tool  for  making  initial 
estimates  of  forces  in  hypersonic  flow.  In  this  ap¬ 
proximation,  the  approaching  stream  is  assumed  to 
lose  completely  its  momentum  component  normal 
to  the  surface,  whereas  the  tangential  momentum 
is  unchanged.  Thus  the  pressure  is  increased  on  a 
segment  of  the  surface  that  can  be  "seen”  by  the 
oncoming  stream,  i.e.,  a  segment  whose  surface- 
normal  makes  an  oblique  angle  with  the  oncoming 
stream.  The  surface  segments  “shaded”  from  the 
oncoming  stream  are  assumed  to  make  a  negligible 
contribution  to  the  force  impressed  on  the  body, 
whereas  those  perpendicular  to  the  stream  are  sub¬ 
ject  only  to  ambient  stagnation  pressure  P,.  In 
principle  the  pressures  over  shaded  or  slightly  in¬ 
clined  surfaces  are  not  accurately  represented  by 
these  assumptions,  however,  the  errors  involved  have 
only  a  secondary  effect  on  the  overall  forces  which 
depend  on  the  integrated  pressures. 

The  impact  theory,  which  arose  from  physical 
reasoning,  has  been  studied  on  a  rational  analytic 
basis  in  limiting  hypersonic  cases  for  which  the  con¬ 
ditions  Mf  — h  and  y  — *-  1  are  considered  ( M,  is 


the  flight  Mach  number  and  y  the  ratio  of  specific 
heats). 

The  aforementioned  reasoning  results  in  the 
local  relation : 

Cv  =  =  2  sin2  p  (7-1) 

(*?) 

where  C„  is  termed  a  pressure  coefficient  (dimen¬ 
sionless)  and  P  is  the  angle  of  inclination  cf  the 
surface  with  respect  to  the  direction  of  u,  .Use  of 
this  relation  reduces  all  force  determinations  to 
geometric  operations.  Its  nonlinear  dependence  on 
the  deflection-measure  p  is  characteristic  of  pressure 
and  force  relations  in  the  hypersonic  regime. 

A  correction  to  Eq.  7-1,  which  takes  into  account 
a  decrease  in  the  value  cf  the  absolute  pressure 
p  predicted  by  Eq.  7-1  as  a  result  of  centrifugal 
forces  which  accompany  body  curvature,  will  not 
be  considered  here  for  several  reasons.  First,  the 
pressures  over  spheres  and  cylinders  have  been  found 
to  be  in  good  agreement  with  experiment  F  Eq.  7-1 
is  modified  for  use  in  the  following  form: 

=  sin2  p  (7-2) 

where  C,  is  the  value  corresponding  to  the  stagna¬ 
tion  point  on  the  body,  and  can  be  estimated  quite 
accurately.  Indeed  for  hypersonic  flow  over  blunt 
bodies  which  produce  strong  bow  shocks,  an  ex¬ 
cellent  approximation  is 

CP,^^  2  (7-3) 

Furthermore,  when  centrifugal  forces  are  con¬ 
sidered,  it  can  be  deduced  that  the  value  cf  p  may 
decrease  to  zero  at  a  point  on  the  body  where  this 
value  is  not  measured.  Attempts  have  been  made 
to  explain  this  result  by  postulating  the  presence 
cf  regions  cf  separated  flow  and  by  contending  that 
the  hypersonic  theory  which  produces  analytically 
the  “Newtonian”  pressure  plus  a  centrifugal  force 
correction  is  invalid  near  the  zero-pressure  point. 
Finally,  on  straight-sided  surface,  such  as  cones, 
the  centrifugal  force  is  zero.  Thus  force  estimates 
in  this  section  are  based  on  either  Eq.  7-1  or  Eq.  7-2. 

7-5  TYPICAL  FORCE  COEFFICIENTS 

The  pressure  drag  coefficient  for  a  cone  at  zero 
angle  cf  attack  is 

CD  =  — 4  =  2  sin2  pe  (7-4) 

PfUr 

2 
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where 


D  is  the  drag 

A  is  the  area  of  the  base  cf  the  cone,  and 
ft  is  the  cone  half-angle. 


It  is  convenient  to  express  the  forces  on  bodies 
at  angle  of  attack  a  in  terms  cf  force  components 
Fn  normal  to  the  body  axis  x  and  Fe  along  the  body 
axis.  The  corresponding  force  coefficients  are.  re¬ 
spectively. 


C,  = 


FCA 

Prur 

2 


and  Cm 


FnA 

Prur 

2 


where  A  is  a  reference  area  usually  taken  to  be  the 
maximum  cross-sectional  area  of  the  body. 

II'  M  represents  the  moment  cf  the  forces  on 
a  body  taken  about  an  axis  at  the  foremost  point 
of  the  body,  the  moment  coefficient  may  be  de¬ 
fined  as 


where  ^  =  arc  sin  (tan  /3,,/tan  a)  signifies  the 
meridian  angle  at  which  the  pressure  has  dropped 
to  zero:  the  value  4*  =  rr/2  represents  the  most 
leeward  generator. 

Relations  7-7  lead  to  the  following  expressions 
for  the  slope  of  the  cone's  lift  curve  and  its  center 
cf  pressure,  respectively : 

/ 'dCj\  =  2  cos2  ft  (7-8) 

V  da  / “-° 

o\ 

CP  -  (7-9) 

where  \  is  the  cone  length. 

The  coefficients  for  the  frustum  of  a  cone  of 
minimum  radius  rm  and  maximum  radius  r„  where 
the  reference  area  A  =  irr%  can  be  obtained  from 
the  cone  values  by  simple  subtraction  of  the  force 
contributed  by  the  initial  section  of  a  sharp  cone 
cf  base  radius  rB.  The  results  may  be  expressed  as 
follows : 


C 


M 


M  A  X 
2 

Pruf 

2 


where  X  is  the  body  length.  The  center  of  pressure 
CP  may  be  defined  as 


=  M_  =  \C„ 

t  S'  C.v 


(7-5) 


The  formulas  which  relate  the  lift  and  drag 
coefficients  CL  and  CD  to  CN  and  Cc  are: 

CL  =  C.v  cos  a  —  Cc  sin  a  ^  ^ 

CD  =  CN  sin  a  C,  cosa 

where  a  is  the  angle  of  attack. 

For  cones  at  angle  of  attack  a  there  arise  the 
following  equations: 

(i)  a  <  ft: 

CN  =  cos'  ft  sin  2a 

Cc  =  2  sin'  ft  "I"  sin3  a(l  —  3  sin'  ft);  (7-7a) 


(ii)  a  >  ft 

CN  =  (cos'  ft  sin  2ft(=)|.3(/  +tt/2) 

"I"  cos  ^(cot  a  tan  ft  +  2  tan  a  cot  ft)] 

Cc  =  +  |)(~)[2  sin2  Pc  +  sin2  «(!  “  3  sin2ft)] 

+  (j-)  cos  f  sin  2 a  sin  2/3c  (7-7b) 


(C  A')  f  tub  turn  (C.v) 

and  similarly  for  C„  i.e., 


rj\ 

rl) 


(Cr)trustum  =  (Cc)cone(  1  -  rl/rl). 


(7-10) 


The  center  of  pressure  measured  from  the  front 
of  the  frustum  is  given  by: 


(7-11) 


where  X,  is  the  length  cf  the  frustum. 

The  force  coefficients  of  a  circular  cylinder  are 
obtained  as  special  cases  of  the  cone-frustum  equa¬ 
tions.  equations  in  the  a  >  ft  case,  with  the  addi¬ 
tional  observations  that 


(rn  +  rm)(rB  -  r,) 
2 


2X 

rB 


tan  /3c 


The  above  relations  are  substituted  into  Eq. 
7-10  with  the  use  of  Eqs.  7-7b,  and  the  limit  ft  — *  0, 
imposed.  The  following  result  is  obtained  for  the 
lateral  surface  of  a  circular  cylinder: 


(7-12) 


It  is  of  interest  to  note  that  an  analysis  taking 
into  account  pressure  relief  due  to  centrifugal  effects, 
which  are  neglected  completely  in  this  chapter,  leads 
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to  a  reduction  of  about  10  percent  in  the  value  of 
CN  given  in  Eq.  7-12.  (See  Ref  7,  p.  85.) 

For  a  flat  surface  normal  to  the  flow, 

CD  =  2  (7-13) 

For  a  sphere  or  a  hemisphere  whose  axis  is  parallel 
to  uf, 

Cn  =  1  (7-14) 


For  a  cone-frustum  at  zero  angle  of  attack  including 
the  drag  of  the  front  face,  and  using  A  =  it rB, 


CD  =  2\?v/  +  2  sin2  ft 


■M 


(7-15) 


where 


1  -  H  =  —  tan  /3„ 


It  is  readily  seen  that  for  a  fixed  base  rB  and  body 
length  \f,  the  frustum  of  minimum  drag  is  not  a 
sharp  cone  but  a  frustum  for  which 

tan  2/3,,  =  (7-16) 

A 

For  example,  for  rB/\,  =  1,  the  minimum  drag 
is  obtained  when  /?„  =  31 .7"  as  compared  to  /?„  =  45" 
for  the  sharp  cone.  For  the  optimum  frustum  in  this 
case  CD  =  0.77,  whereas  the  sharp  cone  gives 
CD  =  1,  as  does  the  hemisphere  for  which  rB/\  =  1 
also. 

There  are  many  ramifications  and  refinements 
connected  with  the  optimization  of  body  forces  under 
a  variety  of  constraints  such  as  fixed  volume,  fixed 
surface,  etc.  They  are  not  considered  to  be  within 
the  scope  of  this  chapter. 

Bodies  of  a  variety  of  shapes  may  be  treated 
approximately  as  composites  of  segments  of  cones, 
cone-frustums,  cylinders,  etc.  However,  in  this  proc¬ 
ess  it  is  probably  most  convenient  to  sum  actual 
forces  and  moments  in  order  to  avoid  confusing 
the  various  reference  areas  A  which  are  used  to 
form  the  different  force  and  moment  coefficients. 

For  sharp-nosed  slender  bodies  of  revolution, 
whose  shape  can  be  expressed  in  terms  of  a  power-law 


the  drag  coefficient  may  be  approximated  by 


CD 


2n{2n  —  1)  no  centrifugal  forces; 

"n2(3n  —  l)(2n  —  l)-1,  including  cen¬ 
trifugal  forces  (7-17) 


For  the  case  of  no  centrifugal  forces,  n  = 
gives  minimum  drag  and  Cn/(rB/\f  =  H;  whereas 


including  centrifugal  forces  the  terms  are  n  =  f 
and  CD/ (rb/\)2  =  f . 

When  these  power-law  bodies  are  at  small  angles 
of  attack,  the  coefficients  are  approximately 


C,v  =  K  sin  2 a 


Cc  =  K(0.05+  sin2  a) 


K  =  1  + 


/1\|>  ~  1) 
W  (r  +  1) 


+ 


(7-18) 


In  practical  use,  slender  bodies  usually  will  be 
blunted  somewhat  because  of  thermal  factors.  There¬ 
fore,  it  is  of  interest  to  return  to  reexamine  Eq.  7-15, 
which  gives  the  drag  coefficient  of  a  blunted  cone, 
to  estimate  the  ratio  of  the  blunt-nose  drag  to  that 
contributed  by  the  conical  afterbody  for  small  cone 
semi-angle  /3e.  The  result  is: 


which  indicates  that  the  effect  of  the  blunted  nose 
is  strongly  dominant  for  small  cone  angles. 


7-6  AERODYNAMIC  CHARACTERISTICS  OF 
A  FAMILY  OF  SPHERICAL-NOSED 
BODIES  WITH  CONVERGING 
CONICAL  AFTERBODIES 

The  aerodynamic  characteristics  of  this  family 
of  bodies  have  been  calculated  using  impact  theory 
by  Dickey.12  The  importance  of  this  report  and  the 
brevity  cf  its  presentation  prompted  the  inclusion 
of  it  in  its  entirety.  Another  NASA  report,  by 
Rainey,13  has  a  set  cf  working  charts  to  assist  in 
the  calculation  of  forces  and  pressure  coefficients  on 
arbitrary  bodies  cf  revolution ;  these  calculations  also 
being  based  on  impact  theory. 

Reproduced  below,  in  its  original  form  is  the 
report  Forces  and  Moments  on  Sphere-Cone  Bodies  in 
Newtonian  Flow.'2 


SUMMARY 

The  static  longitudinal  aerodynamic  charac¬ 
teristics  of  a  family  cf  sphere-cone  combinations 
(fineness  ratios  from  1.0  to  6.0)  were  computed  by 
means  cf  Newtonian  impact  theory.  The  effects  of 
angle  cf  attack,  fineness  ratio,  and  center-of-gravity 
location  are  shown.  The  results  indicate  that,  with 
the  center  of  gravity  at  or  near  the  center  of  volume, 
the  spliere-cone  combinations  are  statically  stable 
at  trim  points  that  provide  low  to  moderate  lift-drag 
ratios.  In  general,  the  lift-drag  ratio  increased  with 
increasing  fineness  ratio.  As  an  example,  with  the 
center  cf  gravity  at  the  center  cf  volume,  the  lift- 
drag  ratio  at  trim  was  increased  from  approximately 
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0.05  to  0.56  bv  increasing  the  fineness  ratio  from 

1.2  to  6.0. 


INTRODUCTION 

Consideration  cf  the  heating,  deceleration, 
guidance,  and  stability  problems  encountered  during 
high-speed  entry  into  the  earth’s  atmosphere  indi¬ 
cates  that  it  would  be  desirable  for  an  entry  vehicle 
to  have  a  blunt  forward  face  and  to  be  statically 
stable  at  a  trim  point  which  would  provide  a  mod¬ 
erate  lift-drag  ratio.  It  would  also  be  desirable  for 
the  launch  configuration  to  have  axial  symmetry. 
Preliminary  analysis  indicated  that  a  simple  body 
cf  revolution  consisting  cf  a  sphere  with  a  con¬ 
verging  conical  afterbody  possesses  these  desirable 
characteristics  and,  consequently,  that  such  shapes 
merit  further  investigation.  An  analytical  study 
utilizing  Newtonian  impact  theory  was  therefore 
undertaken  to  gain  further  knowledge  cf  the  static 
aerodynamic  characteristics  cf  sphere-cone  combi¬ 
nations.  The  purpose  of  this  report  is  to  present  the 
results  of  that  study.  The  effects  cf  angle  cf  attack, 
fineness  ratio,  and  center  cf  gravity  on  the  lift, 
drag,  and  pitching-moment  characteristics  cf  a 
family  of  sphere-cone  combinations  having  fineness 
ratios  from  1.0  to  6.0  are  shown. 
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NOTATION 


Axial  force 

A 

Axial-force  coefficient,  — - 
qS 

Local  axial-force  coefficient  per  unit  length, 

dt\ 

d(x/d)  , 

Drag-force  coefficient,  CA  cos  « r  ■"  CN  sin  a 
Lift-force  coefficient,  CN  cos  a <  —  CA  sin  a 


Pitching-moment  coefficient,  — — 

qba 

Local  pitching-moment  coefficient  per  unit 


length, 


dCm 

d{x/d) 


Normal-force  coefficient,  — 

qb 

Local  normal-force  coefficient  per  unit  length, 


_dCx_ 

d(x/d) 

Pressure  coefficient 
Drag  force 

Diameter  of  the  sphere 

Distance  of  the  center  cf  gravity  from  the 

longitudinal  axis 

Fineness  ratio  of  the  sphere-cone  combina- 


tion,  - 
d 

Lift  force 


Lift-drag  ratio 

Length  of  the  sphere-cone  combination 
Pitching  moment  about  the  center  cf  gravity 
Normal  force  (force  perpendicular  to  the 


longitudinal  axis  measured  in  the  a  plane) 
q  Free-stream  dynamic  pressure 

R  Radius  cf  the  sphere 

r  Local  cross-sectional  radius  cf  the  sphere- 

cone  combination 

S  Frontal  area  cf  the  sphere,  tR2 
x  Distance  along  the  longitudinal  axis  from 
the  nose  cf  the  sphere-cone  combination 
xa  Distance  along  the  longitudinal  axis  from 
the  nose  to  the  center  cf  gravity  cf  the 
sphere-cone  combination 

a  Angle  cf  attack  (angle  between  the  longi¬ 
tudinal  axis  and  the  free-stream  velocity 
vector) 

Radial  angle  measured  in  a  plane  normal  to 
the  longitudinal  axis  (See  Fig.  7-1.) 
du  Radial  angle  denoting  the  upper  limit  cf  the 
body  surface  exposed  to  the  air  stream 
(See  Fig.  7-1.) 

5  Local  body  slope  with  respect  to  the  longi¬ 
tudinal  axis  (See  Fig.  7-1.) 
r?  Angle  between  the  free-stream  velocity  vec¬ 
tor  and  the  perpendicular  to  the  local  surface 
of  the  body 

SPHERE-CONE  COMBINATIONS 

The  geometry  of  the  family  of  sphere-cone 
combinations  considered  in  this  report  is  shown  in 
Figure  7-2.  Each  combination  consisted  cf  a  sphere 
with  a  converging  conical  afterbody  tangent  to  the 
sphere.  The  length  of  the  conical  afterbody  relative 
to  the  diameter  of  the  sphere  was  varied  so  that  the 
resultant  fineness  ratios  cf  the  sphere-cone  combi¬ 
nations  ranged  from  1.0  to  6.0. 


METHOD 

In  order  to  obtain  the  total  force  and  moment 
coefficients  of  a  sphere-cone  combination,  it  is  con¬ 
venient  to  determine  first  the  local  force  coefficients 
acting  on  a  transverse  section.  The  general  expres¬ 
sions  for  the  local  axial-force  and  normal-force 
coefficients  per  unit  length  of  an  arbitrary  body  cf 
revolution  are  given  in  Reference  1  as 

4r  f8“ 

CY  =  ^-tan  d  I  Cri  dji  (1) 

and 

4r  f'3” 

C/  - - r,  I  Cv  sin  d  dfi  (2) 

TTlt  J-t/2 

The  foregoing  equations  are  for  that  portion  cf  the 
body  exposed  to  the  air  stream.  In  the  analysis  the 
portion  of  the  body  shielded  from  the  air  stream  is 
assumed  to  have  a  pressure  coefficient  of  zero  and 
therefore  not  to  contribute  to  the  body  forces. 

The  pressure  coefficient,  C„  of  a  surface  exposed 
to  the  air  stream  is  given  by  Newtonian  theory  as 


C„  =  2  cos2  T) 

where,  for  an  arbitrary  body  of  revolution 

cos  n  =  cosa!  sin  6  —  sin  a  cos  S  sin  0 


(3) 

(4) 
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I 


Side  view 


Transverse  section 

Figure  7-2.  Geometry  cf  Sphere-  Cone  Combinations 


Substituting  Equations  3  and  4  into  Equations  1 
and  2  and  integrating  yields 

(*■  +  l) 

•  (2  cos2  a  sin'  5  sin2  a  cos2  5)  "i"  cos  /3U 
•(sin  2a  sin  25  —  sin2  a  cos2  5  sin  j3„)  j  (5) 

and 

C'N  =  vR,  {f  (/3“  i) sin  2a  sin  25  cos  /3“ 


4r 

=  -yjtan  5 

TVIl 


2  cos2  o/  sin2  5  —  2  sin  2a  sin  25  sin  /3U 


"i"  f  sin2  0/  cos2  5(sin2  /3„  "i"  2). 


(6) 


The  local  pitching-moment  coefficient  per  unit  length 
can  be  expressed  in  terms  cf  the  local  geometric 
parameters  and  the  local  axial-  and  normal-force 
coefficients  as 

CJ  =  (j  -  -  2  tan  5(7/  -  \  C/  (7) 

The  total  axial-force,  normal-force,  and  pitch¬ 
ing-moment  coefficients  presented  in  this  report 
were  obtained  by  graphically  integrating  the  corre¬ 
sponding  local  force  and  moment  coefficients  per 
unit  length  (determined  from  Eqs.  5,  6  and  7  over 
the  total  length  cf  the  configuration  as  indicated 
below: 


ft/d 

CA  =  C’A  d(z/d) 

Jo 

(8) 

ft/d 

Cn  —  1  «  d(x/d) 

Jo 

(9) 

ft/d 

cm  =  c:  d(x/d) 

J  0 

(10) 

RESULTS 

The  static  longitudinal  aerodynamic  charac¬ 
teristics  of  a  family  of  sphere-cone  combinations, 
with  the  center  cf  gravity  located  at  the  center  cf 
volume,  are  presented  in  Figures  7-3  to  7-9.  The 
results  show  that  these  blunt-faced  bodies  cf  revolu¬ 
tion  are  statically  stable  at  trim  points  that  provide 
low'  to  moderate  lift-drag  ratios.  In  general,  for  a 
given  angle  of  attack,  an  increase  in  the  fineness 
ratio  results  in  an  increase  in  the  lift  and  drag 
coefficients  and  in  the  lift-drag  ratio. 

The  trimmed  (C,  =  0)  aerodynamic  charac¬ 

teristics  cf  sphere-cone  combinations,  with  center  cf 
gravity  at  the  center  cf  volume,  are  plotted  as  a 
function  cf  the  fineness  ratio  in  Figure  7-10.  The 
results  show  that  increasing  the  fineness  ratio  from 
1.2  to  6.0  decreased  the  trim  angle  from  approxi¬ 
mately  79"  to  31°  and  increased  the  lift-drag  ratio 
from  approximately  0.05  to  0.56.  It  should  be  noted 
that  because  cf  the  rotational  symmetry  of  the 
sphere-cone  combinations,  the  trim  angles  may  be 
considered  as  angles  in  pitch  or  yaw  or  any  vector 
combination  thereof.  In  order  to  direct  the  lift  force 
in  the  desired  direction,  some  auxiliary  method  of 
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Figure  7-3.  Axial-Force  Coefficients  (f  Sphere- Cone 
Combinations 


Figure  7-4.  Drag  Coeflcients  (f  Sphere-Cone  Com 
binations 


Figure  7-5.  Normal-Force  Coefficients  cf  Sphere- 
Cou  e  Combinations 


Figure  7-6.  Lift  Coefficients  cf  Sphere-Cone  Com¬ 
binations 


control  (such  as  reaction  jets)  would  be  required  for 
positioning  and  maintaining  the  configuration  in  the 
proper  attitude  cf  roll,  pitch,  and  yaw.  If,  however, 
the  rotational  symmetry  of  the  configuration  is 
removed  by  locating  the  center  cf  , gravity  off  the 
longitudinal  axis,  the  lift  force  could  be  placed  in 
the  desired  direction  by  controlling  only  the  roll 
attitude  cf  the  vehicle. 

The  effect  cf  center-of-gravity  location  on  the 
pitching-moment  coefficients  cf  a  representative 
sphere-cone  combination  (fineness  ratio  =  2)  is 
shown  in  Figure  7-11  for  the  complete  360°  angle-of- 
attack  range.  The  results  shown  in  Figure  7-ll(A) 
(effect  cf  longitudinal  shift  of  the  center  cf  gravity) 
indicate  that  the  sphere-cone  combination  trims  at 
a  usable  angle  cf  attack  only  when  the  center  cf 
gravity  is  located  very  near  the  center  cf  volume 
(i.e.,  closer  than  about  15  percent  cf  the  body 
length).  For  other  center-of-gravity  locations,  the 
trim  angle  is  either  0"  or  180"  and,  as  a  result,  no 
lift  would  be  produced.  With  the  center  cf  gravity 
located  on  the  longitudinal  axis  and  near  the  center 
cf  volume,  the  configuration  has  two  stable  trim 
points,  one  at  a  positive  angle  of  attach  and  one  at 
a  negative  angle  of  attack.  However,  the  results 
shown  in  Figure  7-1 1(B)  (effect  cf  transverse  shift 
of  the  center  of  gravity)  indicate  that  by  locating 
the  center  cf  gravity  off  the  longitudinal  axis,  the 
sphere-cone  combination  can  be  made  to  have  only 
one  stable  trim  point. 


Figure  7-7.  Pitching-Moment  Coefficients  cf  Sphere- 
Cone  Combinations  with  Center  of  Gravity  at 
the  Center  of  Volume 


binations 


Figure  7-9.  Drag  Polars  of  Sphere-Cone  Com¬ 
binations 
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mxtigm 


(A)  Longitudinal  Shift  of  Center  of  Gravity  (B)  Transverse  Shift  cf  Center  cf  Gravity  from  the 

Center  cf  Volume 

Figure  7-11.  Effect  cf  Center-of-Gravity  Location  on  the  Pitching-Moment 
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7-7  FRICTIONAL  EFFECTS 

A  fundamental  idea  in  the  understanding  of 
fluid  flow  phenomena  is  the  Prandtl  boundary  layer 
concept.  According  to  this  idea  the  behavior  of  air 
flowing  over  solid  surfaces  can  usually  be  studied 
in  two  parts.  One  is  a  thin  “boundary'  layer”  region 
near  the  surface; thin  in  relation  to  distances  meas¬ 
ured  along  the  surface.  Within  this  region,  shear 
forces  caused  by  adhesion  of  the  air  to  the  surface, 
are  stong,  and  the  relative  velocity  drops  off  sharply 
to  zero  at  the  surface.  The  second  part,  outside  the 
boundary  layer,  may  be  considered  frictionless,  its 
character  being  determined  primarily  by  the  dis¬ 
turbances  caused  by  the  surface  contour. 

The  transmission  of  shear  to  adjacent  layers 
within  the  boundary  layer  is  termed  a  “transport 
process”  and  is  attributable  to  the  exchange  of 
momentum  between  molecules  of  different  velocities 
which  intermix  and  collide.  It  is  characterized  by 
the  “coefficient  of  absolute  viscosity”  n-  Other  trans¬ 
port  processes  important  in  the  boundary  layer 
region  are  heat  conduction  and  diffusion.  Conduction 
involves  an  energy  transport  and  exchange  between 
intermixing  molecules  of  different  kinetic  energies, 
i.e.,  temperatures.  Diffusion  involves  a  transport  of 
mass,  momentum  and  energy  due  to  intermixing  of 
various  species  of  particles  in  a  non-homogencous 
mixture  of  different  fluids.  For  example  dissociated 
air  may  behave  as  a  non-homogeneous  mixture  of 
molecules,  atoms,  ions  and  electrons.  Likewise,  in¬ 
jected  or  ablated  coolants  different  from  air  form 
mixtures  with  the  surrounding  air. 

The  transport  processes  are  important  only  in 
boundary  layers,  wakes  and  shock  waves.  Outside 
the  boundary  layer,  their  explicit  effects  are  neg¬ 
lected.  Under  most  practical  circumstances  the  char¬ 
acteristic  Reynolds  number  of  the  flow  is  not  small, 
being  on  the  order  cf  100,000  or  more;  while  the 
Mach  number  is  not  extremely  high.  In  these  cases 
the  deflection  of  the  boundary  layer  streamlines  does 
not  alter  appreciably  the  effective  body  shape,  since 
the  layer  remains  thin.  Therefore,  if  the  boundary 
layer  does  not  separate  from  the  surface,  the  outer 
flow  values  may  be  obtained  by  first  disregarding 
the  presence  of  the  boundary  layer,  thereby  con¬ 
sidering  the  outer  flow  to  slide  without  friction  over 
the  boundary  layer.  The  pressure  and  velocity  dis¬ 


tributions  over  the  surface  are  obtained  in  this  way 
Then  the  boundary  layer  characteristics  can  be  de¬ 
termined  in  terms  of  the  nonviscous  flow  properties. 
A  primary  feature  of  the  boundary  layer  concept 
is  the  result  that  the  pressure  does  not  vary  along 
surface-normals  within  the  boundary  layer.  This 
permits  the  determination  cf  forces  other  than  shear 
forces  by  neglecting  the  presence  of  the  boundary 
layer. 

Exceptions  to  this  procedure  cf  negligible  vis¬ 
cous  displacement  are  cf  3  main  types: 

(1)  Very  slow  “creeping  motions”  with  Reynolds 
numbers  on  the  order  cf  10  or  less.  In  this 
case  a  boundary7  layer  does  not  form.  However, 
this  situation  is  not  cf  practical  importance 
for  missile  flight. 

(2)  Hypersonic  flows  with  simultaneous  high 
Mach  numbers  and  low  Reynolds  number, 
where  the  low7  Reynolds  numbers  stem  from 
low  deiisities  and  proximity  to  sharp  leading 
edges  and  tips.  This  leads  to  pronounced 
boundary7  layer  growth  arid  displacement  of 
the  outer  flow.  After  first  neglecting  this 
interaction,  the  displacement  effect  can  be 
estimated  and  the  changed  inviscid  properties 
calculated  by  a  second  iteration.  It  is  im¬ 
portant  to  note  that  the  constancy  cf  the 
pressure  along  surface  normals  still  holds  in 
this  case.  A  rough  criterion  for  neglecting  this 
interaction  effect  is  that  the  quantity  M'JrR~U2 
be  much  smaller  than  one.  This  criterion  is 
based  on  the  pressure  change  induced  by  the 
viscous  flow  deflection.  This  change  is  pro¬ 
portional  to  M ,  times  the  induced  flow  de¬ 
flection.  The  flow  deflection  being  proportional 
to  the  boundary  layer  growth-rate  which 
varies  like  M2R~1/X. 

(3)  Under  conditions  cf  rising  pressure,  the  bound¬ 
ary  layer  may  lose  its  energy  and  become 
separated  from  the  surface  causing  a  wake 
to  form  in  the  separated  region,  and  causing 
a  major  change  in  the  inviscid  flow  pattern, 
especially  in  the  wake  region.  An  estimate  cf 
the  separation  point  can  be  made  in  a  first 
approximation  by  assuming  that  separation 
does  not  occur  and  thus  deriving  the  in¬ 
viscid  and  boundary -layer  properties.  These 
will  indicate  the  incipient  separation  point  by 
producing  a  point  of  zero  shear  stress.  The 
shape  and  flow  characteristics  cf  the  separated 
flow  pattern  are  not  readily  amenable  to 
purely  theoretical  analysis  but  can  be  treated 
empirically. 
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CHAPTER  8 

AERODYNAMIC  TESTING 


8-1  INTRODUCTION 

The  experimental  simulation  of  high  speed 
thermal  conditions  encountered  in  hypersonic  bal¬ 
listic  flight  requires  the  use  of  a  variety  of  com¬ 
plementary  devices.  Whereas  the  measurement  of 
forces,  heat  transfer  and  surface  friction  have  dom¬ 
inated  aerodynamic  testing  in  the  past,  additional 
factors  of  current  importance  are  flow  details  which 
bear  on  the  chemical  and  electrical  properties  of  the 
high  energy  air  surrounding  the  vehicle  and  in  its 
wake.  The  latter  properties  are  significant  in  relation 
to  communications,  detection  and  discrimination 
problems,  and  for  some  purposes  require  relatively 
new  testing  techniques. 

The  need  for  tests  on  a  basis  other  than  full- 
scale  is  clear  from  a  financial,  design  and  scientific 
standpoint.  The  economy  of  model  tests  has  been 
well  proven  although  accepted  design  practice  re¬ 
quires  the  acquisition  of  all  available  theoretical  and 
experimental  data  on  a  vehicle  hr  parameter  ranges 
as  close  to  full-scale  as  possible.  Moreover,  detailed 
testing  of  a  type  only  possible  with  models,  leads 
to  improved  performance  which,  in  the  full  scale 
ballistic  missiles,  may  involve  faster  re-entry  with 
increased  payloads,  and  improved  stability  and 
launching  behavior. 

Finally,  the  scientific  aspect  of  the  information 
gathered  in  experiments  gives  rise  to  new  ideas  and 
technical  advancement. 

The  variety  and  severity  of  ballistic  flight  con¬ 
ditions,  as  well  as  the  aforementioned  multiple 
objectives  of  testing,  make  it  clear  that  no  single 
facility  will  suffice  for  all  purposes.  Rather  a  set 
of  coordinated  experimental  facilities  is  required  to 
provide  maximum  flexibility  and  research  capability7. 

The  experimental  aerodynamics  connected  with 
ballistic  vehicles  may  be  divided  into  three  categories : 

(1)  The  first  is  essentially  of  a  fluid-mechanical 
nature,  and  involves  the  pressure  distribution  and 
frictional  resistance  of  the  vehicle  which  give  rise 
to  the  aerodynamic  forces.  These  aspects  may  be 
considered  to  be  influenced  indirectly  by  the  thermo¬ 
dynamics  of  the  air  at  high  speed  and  high  temper¬ 
ature.  In  most  ballistic  missiles,  the  friction  drag  is 


a  relatively  small  fraction  of  the  total  drag,  being 
on  the  order  cf  0.1  or  0.01  of  the  total.  Hence,  it 
does  not  require  very  accurate  evaluation,  and  is  not 
given  prominence  in  testing  programs.  The  remaining 
forces  require  rather  accurate  determination  but  do 
not  necessarily  require  the  simulation  of  the  full- 
scale  thermodynamic  effects. 

(2)  The  second  involves  the  heat  transfer  to 
the  vehicle.  Even  this  type  of  information,  which 
at  first  appears  to  be  so  intimately  connected  with 
the  thermodynamic  behavior  of  the  air,  niay  be 
largely  uncoupled  from  the  explicit  thermodynamic 
and  chemical  effects  and  thus  may  be  dealt  with 
by  the  use  of  various  types  cf  partial  simulation. 

Also  in  this  category  may  be  placed  the  study 
cf  cooling  methods,  and  the  behavior  of  materials 
such  as  heat  sinks  and  ablators  under  various  types 
cf  thermal  loading. 

(3)  The  third  category  involves  the  explicit 
investigation  of  the  chemical  state  of  the  gas  sheath 
around  the  vehicle  and  in  the  wake.  Moreover,  it 
concerns  the  chemical  behavior  of  the  flow  in  various 
parts  of  the  test  facilities  which  utilize  high  energy' 
gases  in  their  operations,  i.c.,  the  chemical  nature 
cf  the  flow  in  shock-tubes  and  in  hypersonic  nozzles 
supplied  by  high  energy'  air  sources  is  cf  importance. 

Fortunately  the  problem  cf  dealing  with  chem¬ 
ical  effects  in  flows  around  vehicles  can  be  simplified 
in  some  cases  by  the  expedient  of  uncoupling  the 
fluid-mechanical  behavior  from  the  chemical  be¬ 
havior.  For  example,  in  many  hypersonic  flows  the 
general  flow-field  configuration,  including  pressure 
distributions,  are  not  affected  strongly  by  the  details 
cf  the  chemical  processes.  Thus,  the  shock  pattern, 
streamline  shapes  and  pressures  can  be  determined 
to  a  first  order,  theoretically  or  experimentally,  by 
taking  into  account  simplified  chemical  and  thermo¬ 
dynamic  effects.  On  the  other  hand  the  chemistry 
may  be  strongly  influenced  by  the  pressure  distribu¬ 
tion  in  the  field  and  in  particular  along  streamlines. 
Thus,  as  a  second  step,  the  pressure  variation  along 
streamlines  may  be  prescribed  and  the  chemistry 
determined  by  one-dimensional  analysis  along  the 
Streamlines.  Experimentally  this  means  that  simple 
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nozzle  flows  may  be  used  to  simulate  the  chemical 
behavior  of  the  flow  along  the  streamline.  It  is  true 
that  this  procedure  does  not  provide  a  full  chemical 
configuration;  however,  it  is  an  important  step  in 
the  simulation  of  a  very  complicated  and  difficult 
flow  property.  Other  aspects  of  this  problem  will  be 
discussed  later. 

8-2  TYPES  OF  TEST  FACILITIES 

Listed  below  are  various  types  cf  aerodynamics 
test  facilities  with  brief  statements  on  the  features 
of  each  type.  Most  require  heat-supply  systems  for 
their  operation.  Sources  of  heat  for  these  facilities 
are  outlined  separately  in  Paragraph  8-3.  A  more 
detailed  description  of  several  of  these  testing  sys¬ 
tems  is  given  in  Paragraph  8-4.  Considered  are: 

(1)  Free  flight-testing  of  prototypes  or  models45 

(2)  Ballistic  firing  ranges.2-8 

(3)  Wind  tunnels,  continuous  or  intermittent,  with 
preheating  systems  to  prevent  test  section 
condensation  of  the  expanded  gas  or  to  achieve 
increased  stagnation  temperatures. 10-17 

(4)  Shock  tunnels.18-22 

(5)  Electrical  discharge  methods  to  generate  a 
stream  of  high  temperature  fluid  (plasma  jets 
and  “hot-shot”)10’23-28 

(6)  Low-density  tunnels.27-29 

Free  flight-testing  usually  is  the  most  expensive 
method.  Unmanned  rocket  vehicles  are  expensive 
and  are  usually  considered  to  be  expendable.  Their 
use  presents  problems  in  the  collection  of  data  which 
must  be  recoverable  or  accurately  telemetered  to 
the  ground.  Furthermore,  accurate  knowledge  and 
control  of  the  flight  plan  and  altitude  cf  the  vehicle 
is  required  but  is  often  difficult  to  acquire  in  practice. 
This  method,  however,  can  supply  test  conditions 
close  to  those  in  actual  flight. 

A  “ballistic  range”  is  generally  a  long  tube 
which  may  be  pressurized  or  filled  with  gases  other 
than  air  and  through  which  models  are  fired  from 
guns.  It  is  a  relatively  economical  device  and  can 
supply  excellent  material  for  visual  study  by  various 
photographic  means.  Electrical,  radiative  and  chem¬ 
ical  effects  cf  the  flow  around  the  body  and  in  the 
wake  can  be  made  with  the  use  of  microwaves, 
photosensors,  and  spectrometers,  respectively.  At 
present,  miniature  recording  and  telemetering  de¬ 
vices  are  being  developed  to  permit  certain  detailed 
measurements  to  be  made.  Size  is  still  a  restricting 
factor,  however,  and  in  this  sense  range  testing  has 
some  of  the  disadvantages  of  flight-testing. 

The  firing  of  a  projectile  upstream  in  a  wind 


tunnel  augments  the  relative  velocity  and  hence  the 
stagnation  temperature.3,9  Certain  aspects  of  actual 
trajectories  can  also  be  duplicated  by  this  means. 

In  “wind  tunnels,”  models  are  usually  held 
fixed  with  respect  to  the  ground.  This  facilitates 
the  making  of  controlled  and  accurate  measure¬ 
ments.  Wind  tunnels  are  usually  found  to  be  econom¬ 
ical  when  the  cost  is  compared  to  the  amount  of 
data  obtained.  The  intermittent  or  blow-down  type 
of  tunnel  is  less  expensive  than  its  continuously 
operating  counterpart.  Duration  of  the  test  run  can 
vary  from  several  seconds  to  minutes  depending 
upon  the  test  conditions  and  the  reservoir  capacity. 
There  are  many  areas  in  which  the  desired  informa¬ 
tion  can  be  obtained  in  a  relatively  short  test  run, 
thus  justifying  this  type  of  tunnel.  In  certain  applica¬ 
tions,  continuous  operation  may  be  desirable,  but 
the  cost  and  power  requirements  may  be  prohibitive. 

In  wind  tunnels  using  air,  Mach  numbers  of 
4  or  5  can  be  achieved  without  preheating  the  air 
supply  to  avoid  condensation  in  the  test  section.10,30 
Otherwise  some  form  of  preheating  is  required.  A 
number  cf  these  preheating  systems  are  discussed 
separately  in  the  next  paragraph. 

High  Mach  number  wind  tunnels  using  certain 
gases  other  than  air,  particularly  helium,  do  not 
require  preheating  to  prevent  condensation  phenom¬ 
ena.  However,  with  these  devices  it  is  necessary  to 
relate  the  flow  phenomena  occurring  in  helium  to 
those  in  air. 

“Shock  tubes”  have  been  used  to  produce  rela¬ 
tively  low  speed  flows,  on  the  order  of  Mach  2, 
at  temperatures  as  high  as  15,000°R.  Flow  durations 
on  the  order  of  1  msec  are  typical  and  the  entire 
measurement  problem  under  these  circumstances  is 
a  difficult  one.  To  overcome  these  difficulties  it  has 
been  proposed  that  a  series  cf  shock  tubes  can  be 
used  successively  to  supply  their  air  to  a  collection 
chamber  which  could  act  as  a  settling  or  supply 
chamber  for  a  supersonic  nozzle.19,33  This  device 
may  modify  the  shock  tubes  disadvantage  cf  short 
running  times. 

Mach  numbers  above  2  may  be  obtained  by 
expanding  the  high  energy  gas  sample  provided  by 
the  shock  tube  through  a  nozzle  attached  to  the 
end  cf  the  shock  tube.  The  augmented  shock  tube 
is  termed  a  “shock  tunnel.” 

It  is  important  to  note  that  “shock  tubes,” 
which  form  the  driving  portions  cf  shock  tunnels, 
also  have  an  extremely  important  use  as  tools  for 
the  study  cf  the  high-temperature  properties  of  gas 
samples,  i.e.,  as  tools  for  physical-chemical  studies 
cf  gas  properties.  The  shocked  gas  samples  may  be 
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examined  for  composition,  reaction  rates  and  radia¬ 
tive  properties  under  high  energy  conditions.  Such 
information  relative  to  the  properties  of  air  and 
possible  gaseous  contaminants  clearly  is  important 
for  hypersonic  flow  studies. 

Electrical  discharges  from  arcs  have  been  used 
to  generate  streams  of  high-temperature  fluid.  The 
temperatures  are  reported  to  be  in  the  order  of 
15,000°R,  but  the  properties  of  the  so-called  “plasma” 
fluid  are  not  well-known.  This  method  appears  to 
be  useful  for  materials  testing.  Gas  stabilized  plasma 
jets  and  fluid  transpiration  arcs  have  been  developed 
which  supply  cleaner  high  energy  streams  at  low 
density. 

Low  density  phenomena  of  importance  at  high 
altitudes  are  studied  by  means  of  special  wind 
tunnels  whose  dominant  components  are  the  vacuum 
exhaust  systems. 

8-3  TYPES  OF  HEATERS 

Among  the  various  types  of  heating  systems 
available  for  aerodynamic  test  facilities,  the  follow¬ 
ing  are  perhaps  the  most  widely  used: 

(1)  Conventional  tubular  or  electric  heaters  for 
air.13 

(2)  Ceramic  storage  heaters  heated  electrically  or 
with  combustion  products.34 

(3)  Compression  heaters."'1'1 

Heating  devices  such  as  induction  heaters  and 
radiation  heaters  are  not  considered  here  inasmuch 
as  fluid  is  not  involved. 

Relatively  large  tubular  or  electric  heaters  are 
currently  available  for  raising  wind  tunnel  stagna¬ 
tion  temperatures.  These  are  now  mainly  in  the 
1500°R  range,  but  an  increase  to  the  2500°R  range 
is  quite  likely  in  the  near  future. 

A  ceramic  heater  intermittenly  operated  and 
electrically  powered  has  provided  air  at  3000°R  and 
600  psi  with  flow  rates  up  to  12  lb/sec.  This  tem- 
perture  can  be  raised  to  4000°R  or  5000°R  if  com¬ 
bustion  firing  is  used. 

By  operating  these  heaters  with  gases  other 
than  air,  the  heat  transfer  rates  can  be  increased 
for  the  same  conditions  of  supply  temperature  and 
pressure.  The  use  of  gases  other  than  air  has  been 
investigated  from  the  standpoint  of  power-saving. 
Some  wind  tunnels  have  been  operated  directly  with 
the  products  of  chemical  reactions  such  as  the  usual 
combustion  exhaust  gases  or  special  reaction 
products. 

Adiabatic  compression  heaters  may  also  be  used 
to  achieve  high  supply  temperatures  and  pressures. 


Such  a  heater  is  exemplified  by  the  light  gas  free- 
piston  gun  tunnel  described  in  Paragraph  8-4.6. 

8-4  DESCRIPTION  OF  SELECTED  TESTING 
FACILITIES 

In  the  following  description  of  test  facilities, 
emphasis  is  placed  on  the  blow-down  wind  tunnel 
with  combustion-fired  ceramic  storage  heater  for 
fluid  mechanical  information  in  all  ranges  cf  rare¬ 
faction  and  on  the  shock  tunnel  for  real-gas  in¬ 
formation. 

8-4.1.  Blowdown  Hypersonic  Wind  Tunnel  with 
Storage  Heater10  34 

Perhaps  the  most  generally  useful  and  versatile 
single  device  for  fluid  mechanical  information  and 
testing  is  the  intermittent  hypersonic  blow-down 
wind  tunnel,  using  axially  symmetric  nozzles  with 
air  supplied  through  a  ceramic  storage  heater 
(Figure  8-1). 34  This  facility  is  a  ground  test  device 
and  has  the  evident  advantages  of  fixed  models  and 
instrumentation, 

Up  to  the  present  tune,  the  maximum  stagna¬ 
tion  temperatures  achieved  in  conjunction  with  con¬ 
sistent  operation  have  been  in  the  3000°R  range. 
It  appears  possible  to  increase  this  maximum  capabil¬ 
ity  of  ceramic  storage  heaters  to  4000°R  or  5000°R 
if  the  elevated  temperature  characteristics  of  the 
ceramic  (e.g.,  zirconium  oxide)  would  be  improved. 
Improved  performance  can  also  be  achieved  F  the 
scarcer  refractories  such  as  beryllium  oxide  could  be 
produced  in  pebble  form  a  t  a  nonprohibitive  cost. 

Even  supply  temperatures  as  high  as  5000°R, 
however,  are  not  sufficient  to  produce  actual  physical 
chemical  effects,  such  as  molecular  dissociation  of 
air  which  is  encountered  in  hypersonic  flight.  The 
elevated  temperatures  and  the  corresponding  en¬ 
thalpies  in  this  range  have  had  the  following  main 
purposes : 

(a)  to  prevent  component  condensation  at  Mach 
numbers  above  4  ; 

(b)  to  provide  temperature  or  enthalpy  differences 
sufficient  to  yield  well-defined  heat  transfer 
measurements  ; 

(c)  to  establish  a  temperature  or  enthalpy  ratio, 
particularly  a  ratio  cf  stagnation-to-surface 
enthalpy,  which  is  an  important  similarity 
parameter  for  determining  boundary  layer  and 
heat  transfer  characteristics; 

(d)  to  establish  an  enthalpy  flux  of  sufficient 
magnitude  and  duration  to  permit  meaningful 
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mass  transfer  coolmgtests  and  tests  of  thermal 
effects  in  structures; 

(e)  to  permit  the  study  cf  interactions  between 
gas  molecules  and  solid  surfaces,  for  example, 
free  molecular  accommodation  coefficients  at 
high  temperatures. 

Mach  number  limits  above  which  air  component 
condensation  is  possible  under  equilibrium  conditions 
at  a  given  stagnation  temperature  T.  are  shown  in 


Figure  8-2. 10  In  this  figure,  the  abscissa  is  the  static 
temperature  T,  corresponding  to  a  given  Mach 
number  M .  It  is  important  to  remember  that  whereas 
condensation  is  a  rate  process,  the  above  limits  are 
based  upon  an  infinitely  rapid  condensation  rate, 
i.e.,  equilibrium  condensation,  and  therefore  rep¬ 
resent  conservative  estimates.  These  limiting  values 
could  therefore  be  exceeded  with  currently  available 
supply  temperatures,  particularly  at  low  density 
levels. 


Figure  8-1.  Schematic  Cross-Sectional  View  cf  Convection  Heater  and  Inlet  Sec¬ 
tion,  Polytechnic  Institute  cf  Brooklyn  Aerodynamics  Laboratory 34 
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Figure  8-2.  Limit  cf  Mach  Number  for  Component  Saturation  Lines 11 


With  the  blow-down  system,  stagnation  pres¬ 
sures  as  high  as  3000  psia  are  considered  quite 
feasible.  There  appears  to  be  no  limit  in  principle 
on  the  extent  to  which  the  stagnation  pressure  can 
be  decreased  for  low  density  testing;  the  capability 
of  the  vacuum  system  being  the  determining  factor 
in  this  connection. 

Supply  pressures  p,  required  for  equal  Reynolds 
numbers  in  the  test  section  and  in  flight  are  shown 
in  Figures  8-3(A)  and  8-3(B)'°  for  T,  =  3000°R  and 
5000°R  and  for  the  following  ratios  of  prototype 
to  model  characteristic  lengths;  L,JLm  =0.1  and  1.0. 
Altitudes  cf  100,000and  200.000  feet  have  been  con¬ 
sidered.  ,Itis  seen  that  Reynolds  numbers  at  altitudes 
around  200,000  feet  and  above  may  be  readily 
achieved  with  reasonable  stagnation  pressures. 

Mass  flows  achievable  in  a  blow-down  hyper¬ 
sonic  facility  may  be  as  high  as  100  lbs/sec.  In 
general,  much  lower  mass  flows  are  required  for  low' 
density  hypersonic  operation  than  for  high  density 
operation.  For  example,  with  p.  —  100  psia,  Ts  = 
3000°R  and  a  mass  flow  cf  2.5  lbs/scc,  a  Mach  14 
tunnel  8.8  feet  in  diameter  can  be  operated.  This 
yields  full  scale  flight  Reynolds  number  simulation 
at  200.000  feet  when  L„/Lm  =  1  (see  Figure  8-3 (A)). 
It  should  be  noted,  however,  that  large  test  section 


dimensions  result  in  long  nozzles — for  the  case  cited 
above,  the  nozzle  length  would  be  about  100  feet 
long. 

Running  times  for  facilities  of  this  type  can 
vary  front  several  seconds  to  several  hours  depending 
upon  the  particular  type  of  lest  or  nozzle  size.  For 
low  density,  low  mass  flow  operation,  runs  may  be 
virtually  continuous,  being  limited  primarily  by  the 
air  supply  and  exhaust  system  capabilities. 

Worthy  of  mention  in  connection  with  the  blow¬ 
down  tunnel  is  the  shrouded  model  technique  which 
has  been  developed  primarily  for  high  Reynolds 
number  testing.  Sketches  cf  the  system  as  applied 
in  forcbody  studies,  and  in  base  and  wake  flow 
studies,  are  presented  in  Figures  8-4  and  8-5.  For 
forehody  flows,  the  shroud  forms  a  channel  around 
the  body  and  thus  impresses  on  it  a  prescribed 
pressure  distribution  which  may  correspond  to  that 
cf  actual  flight,  if  desired.  Boundary  layer  and  heat 
transfer  characteristics  are  thus  produced  at  stagna¬ 
tion  pressures  which  exist  in  the  settling  chamber 
ahead  cf  the  shroud.  Large  models  may,  therefore, 
be  utilized  with  relatively  small  amounts  cf  fluid. 

In  connection  with  afterbodies  and  wakes,  the 
shroud  is  utilized  to  duplicate  a  streamtube  in  the 
inviscid  region  where  shape  may  be  deduced  from 
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Figure  8-4.  Schematic  Diagram  cf  Shrouded 
Mode2  Technique  for  Forebodies  (Reference  31 
cf  Chapter  6) 

calculations  or  small  body  tests.  Proper  flow  condi¬ 
tions  nrior  to  the  afterbody  are  regulated  within 
the  settling  chamber  of  the  shroud. 

Gases  other  than  air  may  be  employed  as  a 
test  media  in  the  blow-down  system  but  they  would 
not  >e  used  in  conjunction  with  a  ceramic  pebble 
ben  heater.  A  simple  gas  such  as  helium  may  be 
used  to  avoid  condensation  or  to  produce  high  heat 
flec-os. 

The  direct  use  of  other  combustion  reactions 
m<  produce  gas  mixtures  resembling  air  has  received 
nine  attention.  It  appears  possible  to  produce  rela¬ 
tively  large  mass  flows  of  test  gas  which  can  be 
supplied  continuously  and  conveniently.  Their  supply 
temperatures,  however,  are  in  the  range  3000°- 
■5000°R  and  do  not  produce  the  dissociation  phe¬ 
nomena  of  hypersonic  flight.  Up  to  the  present  time, 
combustion  products  have  been  used  primarily  as 
a  source  of  high  temperature  gas  for  structural  test¬ 
ing  rather  than  as  a  medium  for  the  study  cf  flow 
field  properties.  A  schematic  diagram  of  a  rig 
utilizing  a  heater  of  this  type  recently  put  into 
operation  at  the  Polytechnic  Institute  of  Brooklyn 
is  shown  in  Figure  8-6. 

It  is  important  to  note  that  under  certain  con¬ 
ditions  there  may  not  be  the  necessity  to  test  at 
extremely  high  Mach  numbers.  For  a  body  enveloped 


in  a  shock  which  is  every  where  strong,  the  flow 
between  the  body  and  the  shock  depends  only  on 
the  free  stream  speed  and  density  and  is  independent 
of  all  other  free-stream  conditions,  including  the 
free-stream  Mach  number.  This  ‘  Mach  number 
independence”  principle  includes  all  real  gas,  trans¬ 
port,  relaxation,  and  rarefied  gas  effects.  Thus  it 
is  seen  that  for  these  studies  around  bodies  enveloped 
in  strong  shocks,  the  full-scale  Mach  number  need 
not  be  duplicated.  It  is  emphasized,  however,  that 
the  principle  applies  to  flows  around  the  same  body 
and  that  further  reasoning  is  involved  in  comparing 
two  bodies  of  different  scale. 

The  Mach  number  independence  principle  can 
also  be  applied  to  hypersonic  boundary  layers.  It  has 
been  demonstrated  that  F  the  fluid  composition, 
pressure  distribution,  surface  temperature  distribu¬ 
tion,  and  the  virtually  constant  external  speed  are 
given,  the  flow  in  major  portions  of  the  boundary 
layer  is  independent  of  the  external  Mach  number. 

An  interesting  application  of  the  blow-down 
wind  tunnel  with  axially  symmetric  nozzle  arises 
in  connection  with  the  firing  cf  a  free  flight  model 
upstream  through  the  tunnel.3,10  The  relative  veloc¬ 
ity  between  the  model  and  stream  creates  stagnation 
enthalpies  and  Reynolds  numbers  close  to  those 
encountered  in  free  flight.  This  method  can  be 
applied  to  duplicate  some  features  of  the  time- 
dependent  changes  occurring  during  the  course  of  a 
trajectory,  e.g.,  the  heat  flux  history. 

It  can  be  concluded  that  there  appears  to  be 
very  few  high  speed  phenomena  whose  dominant 
characteristics  are  fluid  mechanical  (including  heat 
transfer)  rather  than  physio-chemical  which  cannot 
be  investigated  or  tested  for  by  means  of  the  blow¬ 
down  tunnel.  The  similarity  parameters  or  model 
laws  are  quite  well  understood  for  these  types  of 
flows,  and  are  the  same  for  both  high  and  low 
densities,  e.g.,  Mach  number,  Reynolds  number, 
Nusselt  number  and  other  similarity  parameters. 


Figure  8-5.  Schematic  Diagram  of  Shrouded  Model  Technique  for  Afterbodies  and 
Wakes  (Reference  31  cf  Chapter  6) 
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Figure  8-6.  Schematic  Diagram  d  Combustion  Type  Wind  Tunnel  with  Shroud  Rig,  Polytechnic 
Institute  d  Brooklyn  Aerodynamics  Laboratory” 


This  permits  test  data  to  be  readily  interpreted  in 
terms  of  free  flight  conditions. 

In  the  full  use  of  this  facility,  ingenuity  is  a 
factor  not  to  be  disregarded.  A  set  of  test  configura¬ 
tions  employing  high  energy  gas  sources  is  shown 
in  Figure  8-7. 27 

8-4.2.  Shock  Tunnel 

At  air  temperatures  in  the  4000°  to  5000°R 
range  and  higher,  the  effects  of  molecular  dissocia¬ 
tion,  ionization  and  radiation  become  important. 
The  blow-down  system,  just  discussed  is  incapable 
of  developing  supply  energies  sufficiently  high  to 
produce  these  “real  gas”  effects.  In  order  to  obtain 
air  energies  of  proper  order  in  ground  test  facilities, 
it  has  been  found  feasible  to  use  systems  which 
produce  high  energy  bursts  for  very  short  periods 
of  time — on  the  order  of  milliseconds. 

The  most  highly  developed  and  versatile  of 
these  systems  at  present  is  the  shock-tube  wind 
tunnel,  sometimes  simply  referred  to  as  the  shock 
tunnel.  (Figure  8-8  (D) ) . 

Basically  the  shock  tunnel  consists  of  a  shock 
tube  which  provides  a  small  amount  of  high-en¬ 
thalpy  high  pressure  gas,  which  is  then  expanded 
to  high  velocity  through  a  nozzle.  The  shock  tube 
itself  consists  of  a  duct  separated  by  a  diaphragm 
into  regions  of  high  pressure  (driver  gas)  and  low 
pressure  (driven  or  test  gas).  When  the  diaphragm 
is  ruptured,  artificially  or  by  a  further  increase  in 
driver  pressure,  the  driver  gas  expands  into  the  low 
pressure  section,  sending  before  it  a  shock  which 
compresses  the  driven  gas.  The  interface  between 


the  driver  and  driven  gases  follows  the  shock,  and 
in  practice  the  usable  test  gas  is  contained  betwe  n 
the  shock  and  the  interface.  As  the  shock  passes 
through  the  driven  gas,  it  sets  the  compressed  gas 
in  motion  at  a  Mach  number  which  varies  roughly 
between  1.8  and  3.0  depending  upon  the  real  gas 
effects.  A  model  may  be  placed  in  this  high  enthalpy, 
high  pressure,  low  supersonic  Mach  number  gas  in 
one  method  of  using  the  shock  tube,  this  being 
referred  to  as  the  “straight  through”  method 
(Figure  8-8(A)). 

I  n  a  second  method,  the  shock  may  be  permitted 
to  impinge  upon  the  downstream  face  of  the  shock 
tube  and  then  be  reflected  back,  further  compressing 
the  gas  in  its  wake.  A  nozzle  may  be  placed  in  the 
downstream  face  of  the  shock  tube  for  further  expan¬ 
sion  cf  the  small  amount  cf  test  gas  in  both  the 
straight-through  and  reflected  shock  methods.  In 
still  another  arrangement,  the  shock  tube  may  be 
left  open  at  the  downstream  end  so  that  the  test 
gas  is  expelled  in  a  free  jet  from  the  open  end. 

One  cf  the  most  significant  difficulties  in  the 
use  cf  the  shock  tube  is  the  attenuation  of  the  shock 
as  it  travels  down  the  tube.  This  attenuation  is  due 
to  the  viscous  effects  in  a  tube  of  high  length-to- 
diameter  ratio,  the  latter  being  required  to  achieve 
practical  durations  cf  run  in  the  test  gas. 

In  an  important  development,  this  problem  has 
been  greatly  diminished  by  the  use  cf  the  “tailored- 
interface”  modification  cf  the  reflected  shock  method 
(Figure  8-8(D)).  This  permits  the  reflected  shock 
to  pass  through  the  driver  interface  without  further 
reflection,  and  thus  allows  all  the  air  processed  by 
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the  reflected  shock  to  be  used  as  supply  gas  for 
au  expansion  nozzle.  The  increased  utilization  cf 
such  procedures  has  allowed  significant  reductions 
in  the  required  shock  tube  lengths  and  accompanying 
reduced  effects  of  attenuation.  Testing  times  on  the 
order  of  10  times  that  obtainable  with  the  simple 
reflected  shock  are  possible  by  this  means. 

The  many  technological  and  engineering  factors 
involved  in  the  effective  use  of  shock  tunnels  are 
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discussed  in  detail  in  the  literature.  Significant  ad¬ 
vances  in  the  art  have  evolved  since  its  inception 
Measurements  of  heat  transfer  and  pressure  dis¬ 
tributions  over  models ;  visual  observations  of  waves, 
and  some  wake  studies  have  been  made.  Reliable 
measurements  in  the  wake  or  base  region  of  bodies 
are  most  difficult  to  achieve  in  the  shock  tunnel 
because  of  the  relatively  long  times  required  to  pro¬ 
duce  steady  flow  in  these  regions.  Force  balances 
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Figure  8-7.  Test  Configurations  for  Use  with  High  Energy'  Gas  Sources' 
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for  test  runs  of  5  to  10  milliseconds  duration  appear 
to  be  feasible.  In  addition  several  types  of  testing 
configurations  can  be  utilized  in  connection  with  a 
shock  tube  gas  generator.  (Figure  8-7),  hi  which  a 
wide  range  of  high  altitude  test  conditions  are 
attainable. 

An  advantage  of  the  shock  tunnel  lies  in  the 
fact  that  the  test  gas  is  virtually  free  cf  impurities. 


which  may  not  be  the  case  in  systems  such  as  the 
plasma  jet  and  hot-shot.  This  is  because  the  test  gas 
is  not  contaniinated  by  the  driver,  although  particles 
from  the  ruptured  diaphragm  may  pass  down-stream 
and  strike  a  model  mechanically.  Nevertheless,  the 
gas  remains  chemically  pure,  and  diaphragm  par¬ 
ticles  can  be  removed  mechanically  by  means  of  a 
turning  passage  prior  to  the  nozzle. 
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(A)  Schematic  Diagram  of  Shock  Tube 


(B)X—T  Diagram  Showing  the  Progress  of  the 
Shock  Wave  and  the  Expansion  Wave  Following 
the  Diaphragm  Burst  (The  gases  which  were 
originally  separated  by  the  diaphragm  are  sepa¬ 
rated  by  the  contact  discontinuity  (interface).) 


(C) Pressure  Distribution  at  Several  Typical  Times 
in  the  Shock  Tube  (Time  /0  and  t^from  (B).) 
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(D)  Tailored  Interface 
Figure  8-8.  Nomenclature  for  Shock  Tube" 


Ill  the  hot-shot  system,"  on  the  other  hand, 
the  gas  in  the  stagnation  chamber  is  heated  and 
compressed  directly  by  means  of  an  arc  discharge 
immersed  within  the  gas.  This  gas  is  actually  the 
test  gas,  which  is  expanded  through  the  nozzle.  It, 
therefore,  has  the  possibility  cf  being  contaminated 
by  electrode  material  or  of  possessing  extraneous 
ionized  components.  The  same  is  true  of  the  plasma 
jets,  although  these  jets  my  operate  for  time  periods 
on  the  order  cf  seconds.  Moreover,  the  hot-shot 
may  produce  erosion  of  the  nozzle  throat  which  also 
may  increase  the  flow  contamination. 

An  advanced  and  complex  type  of  gas  heating 
system  developed  at  the  Cornell  Aeronautical  Lab¬ 
oratory,  Inc.  consists  cf  a  rotating  bank  cf  shock 
tubes  which  discharge  into  a  collector  to  provide 
a  steady  flow  cf  shock-heated  air  for  an  interval 
on  the  order  of  tens  of  seconds.19'33  This  device  is 
termed  the  "wave  superheater."  Its  aim  is  to  over¬ 
come  the  disadvantage  of  short  running  times  in 
the  shock  tube  and  to  provide  real-gas  test  media 
for  mass  transfer  studies,  tests  cf  thermal  effects  in 
structures  and  materials,  pressure  distribution  and 
force  studies,  and  other  conditions  in  which  extended 
testing  times  are  required  or  are  convenient. 

8-4.3.  Combined  Wind  Tunnel  and  Shock  Tunnel 

A  new  variation  of  test  facilities  is  obtained  by 
the  use  of  the  shock  tunnel  combined  with  a  blow¬ 
down  hypersonic  wind  tunnel  or,  alternatively,  by 
the  combination  of  two  shock  tunnels,  in  a  manner 
under  development  at  the  Polytechnic  Institute  of 
Brooklyn  Aerodynamics  Laboratory.22  The  objective 
is  to  add  energy  to  a  hypersonic  stream  (the  primary 
flow)  which  has  already  been  expanded  to  low  static 
temperature.  This  is  accomplished  by  sending  a 
shock  downstream  along  the  primary  flow.  The 
shock  is  provided  by  a  secondary'  shock  tube  whose 
downstream  end  is  attached  obliquely  to  the  test 
section  of  the  primary  nozzle  through  which  its  shock 
discharges  into  the  primary  stream.  It  is  important 
to  note  that  the  stagnation  temperature  cf  the 
primary  flow  may  be  kept  sufficiently  low  (e.g.,  less 
than  3000°R)  to  avoid  the  introduction  cf  real  gas 
effects.  Thus  the  stagnation  enthalpy  and  limiting 
velocity  of  the  basic  flow  are  increased  by  the  second¬ 
ary  shock,  while  the  resultant  test  gas  is  not  signifi- 
icantly  influenced  by  chemical  history,  since  its 
static  temperature  is  always  maintained  at  a  rela¬ 
tively  low  level. 

This  method  is  of  importance  since  it  appears 
that  under  low  density  conditions  it  will  not  be 
possible  to  maintain  thermodynamic  equilibrium 


hypersonic  nozzle  flows  which  expand  their  test 
media  from  a  high  energy  source  in  the  real  gas 
range,  regardless  of  whether  the  source  is  a  shock 
tube,  hot-shot,  or  plasma  jet.  It  is  believed  that 
the  nonequilibrium  state  in  the  nozzle  will  not  be 
of  extreme  importance  for  fluid  mechanical  or  heat 
transfer  effects,  provided  that  the  actual  state  and 
velocity  conditions  in  the  nonequilibrium  flow  in 
the  test  section  are  known  and  used  to  explain  the 
data.  For  chemical  and  physical  problems  in  the 
shock  layer  and  wake,  however,  it  inay  be  important 
to  approach  the  precise  thermodynamic  state  of  a 
prototype  in  a  model  test.  For  this  purpose,  a  method 
cf  creating  flows  simulating  the  equilibrium  condi¬ 
tions  of  flight  in  a  ground  test  facility  becomes 
extremely  useful.  This  problem  may  also  be  partially 
overcome  by  the  use  of  high  velocity  guns  and 
firing  ranges,  which  form  useful  adjuncts  to  the 
fixed  model  techniques. 

Another  use  of  this  method  lies  in  the  fact  that 
the  stagnation  enthalpy  level  of  a  basic  flow  can 
be  increased  by  the  action  of  the  superimposed 
shock.  In  principle,  the  effect  of  several  shocks  may 
be  superimposed  in  this  way  to  provide  a  flow  with 
extremely  high  stagnation  enthalpy  or  limiting  ve¬ 
locity.  The  difficulty  lies  in  the  effective  running  time 
which  is  diminished  with  each  superimposed  shock. 

8-4.4.  Hot-Shot  Tunnel 

In  this  system,  electrical  energy  stored  in 
capacitor  bank  (which  may  have  a  capacity  of  a 
million  joules  or  more)  is  released  in  an  arc  chamber, 
and,  through  a  constant-volume  process,  produces 
a  high  pressure,  high  energy  air  mass.  For  example, 
the  power  supply  in  use  at  AEDC"  consists  cf  twin 
47-ton  fly  wheels  driven  by  a  pair  cf  lOOOhp  motors, 
which  operate  a  1.28  X  108  joule  induction  coil. 
A  plastic  diaphragm  upstream  of  the  nozzle  throat 
is  ruptured  and  expansion  into  the  nozzle  takes  place. 

Among  the  advantages  of  this  type  cf  test 
facility  is  the  relatively  long  duration  cf  the  test 
run  compared  to  that  cf  the  shock  tunnel.  This 
permits  sufficient  time  for  the  measurement  of  heat 
transfer,  pressures,  forces,  and  moments.  The  hot- 
shot  tunnel  also  provides  a  wide  operating  range 
and  convenience  of  operation. 

A  critical  feature,  previously  mentioned,  is  the 
impurity  of  the  test  gas  caused  by  the  arc  in  the 
explosion  chamber  and  by  throat  erosion.  However, 
substantial  improvements  in  the  purity  of  the  test 
gas  have  resulted  from  recent  work.  In  this  method, 
expansion  cf  the  test  gas  from  a  constant-volume 
stagnation  chamber  causes  an  unsteadiness  in  the 
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stagnation  conditions.  This  unsteadiness  can  be 
diminished  by  increasing  the  chamber  volume  for 
a  given  length  of  run  but  at  the  expense  of  an 
mcrease  in  the  power  required  for  operation. 

8-4.5.  Plasma-Jet 

The  arc-powered  jet  is  a  device  in  which  a  high 
intensity  arc  formed  between  an  anode  and  cathode, 
stabilized  by  a  gas  which  is  passed  over,  around 
or  through  the  arc  (Figure  8-9),  results  in  a  high 
enthalpy  ionized  gas  jet  consisting  of  the  heated 
stabilizer  plus  residue  from  the  electrodes.  Some 
modern  arc  techniques  utilize  wall-cooling  methods 
for  stabilization;  these  do  not  require  the  addition 
of  a  foreign  fluid  to  the  jet. 

At  present,  the  power  input  requirements  are 
high  for  relatively  small  mass  flows.  For  example, 
in  one  jet,  700  kw  of  AC  power  are  required  to 
produce  0.080  lb/sec  of  "air"  at  a  stagnation  tem¬ 
perature  of  about  11,000°R  and  a  stagnation  pres¬ 
sure  of  150  psia.  A  2000  kw  version  is  cited  as 
yielding  0.30  Ib/sec  of  test  fluid  at  a  stagnation 
pressure  of  2000  psia  and  temperature  of  13,000°R. 
Running  times  on  the  order  of  seconds  or  even 
minutes  may  be  achieved. 

The  plasma  jet  provides  an  interesting  high 
enthalpy  ionized  medium  for  testing  materials  and 
structures.  It  is  particularly  useful  for  ablating  sys¬ 
tems  which  are  primarily  sensitive  to  the  enthalpy 
absorbed  per  pound  of  material  rather  than  to  the 
fluid  mechanical  and  other  detailed  parameters  of 
the  jet.  One  of  the  difficulties  in  the  operation  of 
plasmas  for  aerodynamic  testing  has  been  minimizing 
contamination  of  the  stabilizing  gas  which  forms 


the  test  medium.  This  contamination  stems  primarily 
from  the  anode,  usually  carbon,  rather  than  the 
cathode,  which  may  be  water-cooled  copper.  De¬ 
creasing  the  power  input,  and  the  arc  chamber 
pressure  (stagnation  pressure)  serves  to  diminish 
anode  material  loss.  In  addition,  anode  cooling 
methods  have  been  utilized  to  advantage  to  deter 
material  loss.  Thus  clean  plasma  jets  become  more 
practical  in  aerodynamic  testing  under  low  density 
conditions. 

Theplasmajethasthe  sameproblem,with  regard 
to  non-equilibrium  conditions  within  the  hypersonic 
expansion  nozzle,  as  has  been  mentioned  in  connec¬ 
tion  with  the  shock  and  hot-shot  tunnels.  Under 
low  density  conditions  the  nozzle  flows  may  be 
expected  to  be  frozen.  Moreover,  a  problem  of  un¬ 
steadiness  (small  scale)  exists  under  many  condi¬ 
tions.  The  use  of  settling  chambers  between  the 
arc  and  expansion  nozzle,  to  minimize  non-uni¬ 
formity'  and  unsteadiness,  results  in  substantial 
heat  loss. 

Plasmas  have  also  provided  a  useful  medium 
for  the  study  of  the  behavior  of  ionized  media  and 
for  magnetohydrodynamic  investigation. 

8-4.6.  Light  Gas  Free-Piston  Gun  Tunnel 

An  interesting  device  for  generating  a  short 
burst  of  high  temperature,  high  pressure  air  is  a 
compression  heater  superficially  resembling  a  shock 
tunnel.  In  this  device,  termed  a  gun  tunnel,14  a  free 
piston  is  interposed  between  the  driver  and  driven 
gas,  as  shown  in  Fig.  8-10.  High  thermal  energy 
on  the  downstream  end  of  the  piston  is  achieved 
if  a  light  gas  such  as  helium  or  hydrogen  is  used 
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Figure  8-9.  Plasma  Generator 26 
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Figure  8-10.  Schematic  Diagram  cf  Light  Gas-Free  Piston  Gun  Tunnel 14 


in  the  breech  as  a  driver,  since  the  high  temperature 
results  from  multiple  shock  reflections  which  take 
place  between  the  moving  piston  arid  the  down¬ 
stream  end  of  the  barrel.  The  temperatures  attain¬ 
able  are  higher  than  those  corresponding  to  isentropic 
compression.  In  principle,  temperatures  as  high  as 
10,000°K  can  be  achieved  after  3  shock  reflections 
with  initial  pressure  ratios  on  the  order  of  600  across 
the  piston.  Running  times  on  the  order  of  hundreds 
of  milliseconds  can  be  obtained.  However,  the  erosion 
of  the  nozzle  throat  restricts  the  running  time  during 
high  stagnation  pressure  runs. 

This  device  appears  to  be  quite  attractive  for 
hypersonic  testing,  at  least  on  the  pilot  scale  for 
which  information  is  available.  Factors  such  as 
reliability  and  steadiness  will  influence  final  evalua¬ 
tions  cf  this  device. 

8-4.7.  Firing  Ranges 

A  wide  variety  of  techniques  for  studying 
models  fired  in  free  flight  has  been  accumulated  in 
ballistic  studies  over  the  years. 

The  possibilities  for  flight  simulation  are  ex¬ 
cellent  in  firing  ranges.  However,  a  disadvantage 
is  the  difficulty  in  obtaining  data  on  the  flow  near 
the  model,  or  effects  on  the  model  without  the  use 
cf  fixed  instrumentation.  In  wake  studies  this  factor 
is  not  as  important. 

Past  observation  methods  have  been  primarily 
photographic.  More  recently,  small  telemetering 
devices  have  also  been  developed  for  making 
rudimentary  heat  transfer  and  pressure  measure¬ 


ments  on  the  model.  These  devices  must  withstand 
accelerations  as  high  as  200,000  g  in  some  cases. 
A  ample  resistance,  inductance,  capacitance  circuit 
developed  at  the  U.S.  Naval  Ordnance  Laboratory 
responds  to  an  intermittent  radar  pulse  by  reemitting 
a  decaying  signal.  The  circuit  resistance,  a  measure 
of  temperature  or  pressure,  can  be  estimated  from 
the  decay  rate  of  the  signal.  This  system  is  inexpen¬ 
sive,  rugged,  and  can  be  applied  to  multiple  channels. 

Projectile  studies  have  become  even  more  im¬ 
portant  with  the  advent  of  aerophysical  problems 
connected  with  the  communication  through  and 
detection  and  diagnosis  of  high  energy  flows  in  the 
shock  layer  and  in  the  wake  cf  hypersonic  vehicles. 
These  problems  have  been  studies  with  the  use  cf 
optical  techniques  such  as  the  schliercn  and  shadow¬ 
graph  systems,  electro-magnctic  wave  diagnostics 
(usually  in  the  microwave  range),  molecular  beams, 
spectroscopy,  radiation  sensors  and  other  methods 
of  experimental  physics.  Wake  measurements  can 
be  made  at  points  fixed  in  laboratory  space  over 
a  period  of  time  after  the  model  has  passed. 

The  modern  development  of  light  gas  guns2-6'9 
has  led  to  projectile  velocities  on  the  order  cf 
30,0(X)  fps  with  certain  types  cf  models.  Problems 
arise  in  connection  with  the  launching  of  light 
delicate  models  without  injury  to  the  model  and 
in  connection  with  the  stabilization  cf  slender  pro¬ 
jectiles.  A  typical  modern  flight  range  for  purposes 
of  aerophysical  measurements  may  consist  cf  a  gun 
section  on  the  order  cf  50  ft  in  length,  a  measure¬ 
ment  section  about  50  ft  long  and  an  enlarged 
terminal  section  for  radar  measurements. 
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